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SECTION  I 

STIFFNESS  OF  UNIDIRECTIONAL  COMPOSITES 


SCOPE 

The  stiffness  of  unidirectional  composites,  like  any  other  structural 
material,  can  be  defined  by  appropriate  stress-strain  relations.  We  will 
show  that  the  coefficients  or  material  constants  of  these  relations  can  be 
packaged  in  a set  of  engineering  constants,  compliance  components,  or 
modulus  components.  The  components  of  any  one  set  are  directly  express- 
ible in  terms  of  the  components  of  the  other  sets.  Each  set,  however, 
possesses  unique  characteristics  that  make  it  suitable  for  specific  usage. 

We  intend  to  demonstrate  that  the  stiffness  of  unidirectional  composites  is 
governed  by  the  same  stress -strain  relation  that  is  valid  for  ordinary  mate- 
rials. The  number  of  independent  constants  are  four  for  composites  and  two 
for  ordinary  materials.  We  will  show  the  stiffness  data  for  several  typical 
unidirectional  composites.  We  can  then  compute  stress  from  strain,  or 
strain  from  stress.  We  intend  to  show  that  composites  are  conceptually  and 
operationally  as  simple  as  ordinary  materials.  Two  key  relation*  are: 
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NOMENCLATURE 


Young's  modulus  for  isotropic  materials 
Longitudinal  modulus 

Transverse  modulus 

Shear  modulus  for  isotropic  materials 
Longitudinal-transverse  shear  modulus 

Dimensionless  multiplying  constant  = £l  - 

Modulus  components;  i,  j = 1,  Z,  6 


Compliance  components;  i, j = 1,  Z,  6 


u 

v 

w 

ri 

si 

v 

V 

LT 

v 

TL 
Sub  1 
Sub  2 
Sub  6 


= Displacement  along  x-axis 
= Displacement  along  y-axis 
= Stored  elastic  energy 
= Stress  components;  i,  j = 1,  Z,  6 

= Strain  components;  i,  j = 1,  Z,  6 

= Poisson's  ratio  for  isotropic  materials 
= Major  Poisson's  ratio 

= Minor  Poisson's  ratio 


= Normal  component  along  1 -axis 
= Normal  component  along  2-axis 
= Shear  component  in  1 -2  plane 
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STRESS 

Stress  is  a measure  of  internal  forces  within  a body.  This  together 
with  strain  are  the  key  variables  for  the  determination  of  stiffness  and 
strength  of  a material.  The  mechanisms  of  deformation  and  failure  are  also 
interpreted  in  terms  of  the  state  of  stress  and  strain.  They  are  the  funda- 
mental variables  for  the  mechanical  behavior  of  materials  similar  to  tem- 
perature and  heat  flux  for  heat  conduction,  and  pressure,  volume,  and  tem- 
perature for  gas. 

There  is  no  direct  measurement  for  stress.  Instead,  stress  is  inferred 
or  derived  from  the  following: 

• Applied  forces  using  stress  analysis. 

• Measured  displacements  also  using  stress  analysis. 

• Measured  strains  using  stress-strain  relations. 

When  we  talk  about  stress  we  usually  mean  the  average  stress  over  some 
physical  dimension.  This  is  similar  to  population  measured  over  a city, 
county,  or  state.  In  our  study  of  composites  we  deal  with  three  levels  of 
average  stress: 

• Micromechanical  stress  0~  is  that  calculation  based  on  distinct 
phases  of  fiber,  matrix,  and  in  some  cases  the  interface  and 
voids. 

e Ply  stress  7 is  that  calculation  based  on  assumed  homogeneity 
with  each  ply  or  ply  assembly  in  which  the  fiber  and  matrix 
are  smeared. 

• Laminate  stress  or  stress  resultant  N/h  is  the  average  of  ply 
stresses  across  the  thickness  of  a laminate. 

In  Fig.  1 we  show  two  levels  of  this  idealization  of  average  stresses.  On  the 
micromechanical  level  in  (a),  the  fiber  and  matrix  stresses  vary  from 
point  to  point  within  each  constituent  phase.  The  average  of  these  stresses 


is  the  ply  stress  0” . In  a laminate  or  on  the  macromechanical  level,  each 
ply  or  ply  assembly  has  its  own  ply  stress.  The  average  of  several  ply 
stresses  is  the  laminate  stress  or  stress  resultant  N/h. 


Fig.  1.  Schematic  relations  between  local  and  average  stresses: 

(a)  Micromechanical  level  where  stresses  in  fiber  and 
matrix  are  recognized. 

(b)  Macromechanical  level  where  stresses  in  plies  and 
ply  assemblies  are  recognized. 

We  will  use  contracted  notation  in  this  book.  Single  subscripts  for 
stress  and  strain,  and  double  subscripts  for  compliance  and  modulus  will 
be  followed.  The  conversion  from  the  conventional  or  tensorial  notation 
to  the  contracted  notation  is  shown  in  Table  1. 


Table  1.  STRFSS  COMPONENTS  IN  CONTRACTED  NOTATION 


CONVENTIONAL  OR 

TENSORIAL  NOTATION 

CONTRACTED 

NOTATION 

<T 

(T 

cr 

<T 

0" 

X 

XX 

1 

11 

1 

<T 

<T 

01 

a;. 

01 

y 

yy 

2 

22 

2 

<r 

*y 

<r 

xy 

°I2 

r6 

The  use  of  subscript  6 for  the  shear  stress  component  is  derived  from  the 
6 components  in  3-dimensional  stress.  Although  subscript  3 has  occasionally 
been  used  for  this  shear  component,  it  is  a source  of  confusion  since  3 can 
be  used  for  the  3rd  normal  stress  component  in  3-dimensional  problems. 
Subscript  6 is  used  to  avoid  this  confusion. 

The  state  of  stress  in  a ply  or  ply  assembly  is  predominantly  plane 
stress.  The  nonzero  components  of  plane  stress  are  those  listed  in  Table  1. 
The  remaining  three  components  are  of  secondary  and  local  nature,  and  will 
not  be  treated  in  this  book.  It  is  convenient  to  represent  the  state  of  plane 
stress  in  a 3-dimensional  stress-space  where  the  three  orthogonal  axes 
correspond  to  the  three  stress  components.  The  stress-space  is  shown  in 
Fig.  2.  Here  each  applied  stress,  represented  by  three  stress  components, 
can  be  readily  portrayed  as  a vector  in  a 3-dimensional  space.  The  unit 
vector  which  signifies  the  direction  of  the  applied  stress  are  represented  by 
the  conventional  notation  of 

(i#  y 

where  the  components  of  the  unit  vectors  are  directional  cosines.  All  three 
unit  vectors  are  shown  in  Fig.  2.  Typical  unit  vectors  for  simple  states  of 
stress  will  be  shown  in  the  following  table: 


Table  2.  UNIT  VECTORS  FOR  SIMPLE  STRESS  STATES 


Type  of  Stress 

Unit  Vector 

Longitudinal  tension 

* 

o 

V 

o 

Longitudinal  compression 

1 

N* 

% 

o 

w 

© 

Transverse  tension 

(0,  1,  0) 

Transverse  compression 

(0,  -1,  0) 

Positive  longitudinal  shear 

(0,  0,  1) 

Negative  longitudinal  shear 

(0,  0,  -1) 

Figure  No 


space.  Unit  stress  vectors  compression.  The 

are  also  shown  as  arrows.  respective  unit  vectors 


are  (±1,  0,  0). 

I 


Fig.  4.  Transverse  uniaxial 
stress  in  tension  and 
compression.  The 
respective  unit  vectors 
are  (0,  *1,  0)  . 


Fig.  5.  Positive  and  negative 
longitudinal  shears. 
The  respective 
unit  vectors  are 
(0,  0,  ±1) . 
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2.  STRAIN 

Strain  is  a purely  geometric  quantity;  no  material  property  or  conser- 
vation principles  are  involved.  First  we  will  show  the  geometric  relations 
between  strain  components  and  displacements.  Then  we  can  establish  the 
stress-strain  relation.  The  constants  in  this  relation  govern  the  stiffness 
of  composites  as  well  as  ordinary  materials. 

Relative  displacements  in  a plane  will  induce  2-dimensional  strain. 

If  the  displacements  do  not  vary  from  point  to  point  within  a material,  there 
will  only  be  rigid  body  motion  and  no  strain.  Thus,  strain  is  simply  the 
spatial  variation  of  the  displacements. 

Let  Au  = Relative  displacement  along  x-axis 
Av  = Relative  displacement  along  y-axis 

From  Fig.  7 we  can  define: 

_ lim  Au 
1 Ax»0  Ax 

_ lim  Av 
e2  “Ay^O  Ax 


Fig.  7.  Normal  strain  and  displacement  relations 

The  partial  differentiation  is  used  because  the  displacements  are 
functions  of  both  x and  y coordinates.  Strain,  like  stress,  is  a local  prop- 
erty. In  general  it  varies  from  point  to  point  in  a material.  Only  in  special 


du 

6x 


&v 

by 


8 


i 


! 


cases  is  the  state  of  strain  or  stress  uniform;  we  call  this  homogeneous  strain 
or  stress.  This  special  case  is  pertinent  to  testing  for  property  determination 
where  we  deliberately  try  to  create  a simple,  homogeneous  strain  or  stress. 


Note  that  the  normal  strain  components  are  associated  with  changes  in 
the  lengths  of  an  infinitesimal  element.  The  rectangular  element  before 
deformation  remain  rectangular  although  the  lengths  of  sides  may  change. 
There  is  no  distortion  produced  by  the  normal  strain  components.  Distortion 
is  measured  by  the  change  of  angles.  The  original  rectangular  element  would 
be  distorted  into  a parallelogram.  Geometrically  this  is  equivalent  to  stretch- 
ing one  diagonal  and  compressing  the  other.  This  combined  action  will  pro- 
duce distortion  which  is  measured  by  shear  strain.  Fig.  8 shows  the  com- 
bined action  produced  by  the  same  two  incremental  displacements  that  pro- 
duced the  normal  strain  components  in  Eq.  1.  The  desired  shear  strain  is: 

66  = a + b (2) 


Fig.  8.  The  strain-displacement  relation  for  shear  strain.  The 
arrows  show  the  direction  of  stretching  and  compressing 
the  diagonals. 
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The  resulting  strain  displacement  relation  is 


i 


(4) 


This  is  the  engineering  shear  strain  which  is  twice  the  tensorial  strain. 
Engineering  shear  strain  is  used  because  it  measures  the  total  change  in 
angle,  or  the  total  angle  of  twist  in  the  case  of  a rod  under  torsion.  This 
factor  of  2 is  often  a source  of  confusion.  When  in  doubt,  the  strain- 
displacement  equation,  in  Eq.  4,  is  the  best  place  for  reference. 


As  with  stress,  contracted  notation  will  be  used  for  strain  compo- 
nents. The  conversion  table  between  the  components  of  the  conventional  or 
tensorial  strain  and  the  contracted  strain  will  be  listed  in  Table  3. 


Table  3.  STRAIN  COMPONENTS  IN  CONTRACTED  NOTATION 


CONVENTIONAL  OR 

TENSORIAL  NOTATION 

CONTRACTED 

NOTATION 

6 

X 

6 

XX 

G1 

en 

G1 

ey 

6 

yy 

G2 

G22 

G2 

26 

xy 

2€r 

xy 

2e!2 

2G12 

G6 

Strain  vectors  can  also  be  portrayed  in  strain- space.  Because  of  the 
coupling  between  the  normal  strain  components  known  as  the  Poisson' s effect, 
the  response  to  a uniaxial  stress  creates  a biaxial  strain  state.  For 
example,  for  ordinary  material  an  extension  is  coupled  with  a lateral 
contraction  if  the  applied  uniaxial  stress  is  tensile.  In  Fig.  9 we  will 
show  the  unit  strain  vectors  as  the  result  of  uniaxial  longitudinal  and 
transverse  tensile  stresses  in  (a)  and  (b),  respectively:  If  the  applied 
stress  is  compressive,  the  direction  of  all  the  stress  and  strain  unit 
vectors  shall  be  reversed. 
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(°)  (b) 


Fig.  9.  Unit  strain  vectors  resulting  from  uniaxial  stresses 

(a)  Biaxial  strain  (1,  -v,  0)  resulting  from  uniaxial 
stress  (1,  0,  0). 

(b)  Biaxial  strain  (- v , 1,  0)  resulting  from  uniaxial 
stress  (0,  1,  0) 

3.  STRESS-STRAIN  RELATIONS 

We  will  limit  the  composites  of  this  report  to  the  linearly  elastic 
materials.  The  response  of  materials  under  stress  or  strain  follows 
a straight  line  up  to  failure  in  stress  or  strain  space.  Assumed  linearity 
moreover  permits  us  to  use  superposition  which  is  a very  powerful  tool. 

For  example,  the  net  result  of  combining  two  states  of  stress  is  precisely 
the  sum  of  the  two  states- -no  more  and  no  less.  The  sequence  of  the  stress 
application  is  immaterial.  We  can  assemble  or  disect  components  of  stress 
and  strain  in  whatever  pattern  we  choose  without  affecting  the  result. 
Combined  stresses  are  the  sum  of  simple  stresses. 

Secondly,  elasticity  means  full  reversibility.  We  can  load  and  unload 
and  reload  a material  without  incurring  any  hysteresis.  Elasticity  also 
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means  that  the  material's  response  is  instantaneous.  There  is  no  time 
lag,  no  time  or  rate  dependency. 

Experimentally  observed  behavior  of  composites  follows  closer  to 
linear  elasticity  than  nearly  all  metals  and  non -reinforced  plastics.  The 
assumed  linear  elasticity  for  composites  appears  to  be  reasonable.  If  we 
are  to  improve  upon  this  simplification,  such  as  the  incorporation  of  non- 
linear elasticity,  plasticity  and  viscoelasticity,  the  degree  of  complexity  is 
increased  significantly  beyond  the  scope  of  this  report. 

For  unidirectional  composites,  the  stress- strain  relations  can  be 
derived  by  the  superposition  principle.  We  must  recognize  that  two  orthog- 
onal planes  of  symmetry  exist  for  unidirectional  composites:  one  plane  is 
parallel  to  the  fibers,  and  the  other,  transverse  to  the  fibers.  Symmetry 
exists  when  the  structure  of  the  material  on  one  side  of  the  plane  is  the 
mirror  image  of  the  structure  on  the  other  side.  The  two  orthogonal  planes 
are  shown  in  Fig.  10,  where  the  l-axis  is  along  the  longitudinal  direction  of 
the  fiber  while  the  2-axis  is  in  the  transverse  direction.  When  the  reference 
axes  1-2  coincide  with  the  material  symmetry  axes,  we  call  this  the  on-axis 
orientation.  The  stress-strain  relation  in  this  section  is  limited  to  this 
special  case,  which  is  shown  in  Fig.  10.  The  off-axis  orientation  will  be 
discussed  in  later  sections. 


Fig.  10.  Two  orthotropic  planes  of  symmetry  of  unidirectional 
composites.  Axes  1-2  coincide  with  the  longitudinal 
and  transverse  directions.  This  material  symmetry 
is  called  orthotropic  and  on-axis. 
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The  on-axis  stress-strain  relation  can  be  derived  by  superpositioning 
the  results  of  the  following  simple  tests: 
a.  Uniaxial  Longitudinal  Test 

The  applied  uniaxial  stress  and  the  resulting  biaxial 
strain  were  shown  in  Fig.  9a.  The  stress-strain  curves 
for  this  test  are  shown  in  Fig.  11,  from  which  we  can 
establish  the  following  stress-strain  relations: 


e = - — cr 

1 E 1 

li 


(5) 


e2  -- 


LT 


where  E = Longitudinal  Modulus 


LT 


= Major  Poisson' s Ratio  = - 


Fig.  11.  Uniaxial  longitudinal  tensile  test 

b.  Uniaxial  Transverse  Test 

The  applied  uniaxial  stress  and  the  resulting  biaxial  strain 
were  shown  in  Fig.  9b.  The  stress-strain  curves  for  this 
test  were  shown  in  Fig.  12,  from  which  the  following  stress- 
strain  relations  can  be  established: 
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where  E„  = Transverse  Modulus 
T 


v - Minor  Poisson's  Ratio  = — 

TL  6 


Fig.  12.  Uniaxial  transverse  tensile  test 
c.  Longitudinal  Shear  Test 

We  apply  another  simple  state  of  stress,  the  pure  shear, 
to  our  unidirectional  composite.  This  is  shown  in  Fig.  13. 
The  resulting  stress-strain  relation  is: 

« 1 


(7) 


where  G_  _ * Longitudinal  Shear  Modulus 
Xj  X 
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Fig.  13.  Longitudinal  shear  test 

By  applying  the  principle  of  superposition,  we  can  sum  up  the  con- 
tribution of  each  stress  component  in  Eq.  5,  6,  and  7 to  the  resulting 
strain  components.  The  final  stress-strain  relation  for  our  unidirectional 
composite  is: 


This  is  the  on-axis  stress-strain  relation  of  a unidirectional  compos- 
ite; i.  e.,  the  material  is  in  its  orthotropic  symmetry  orientation. 

These  simultaneous  equations  can  be  repackaged  in  a matrix  multiplication 
table,  wherein  each  row  in  the  table  is  equal  to  the  sum  of  products  from  each 
column  and  its  column  heading.  This  rule  should  be  self-evident  if  we  com- 
pare the  first  of  Eq.  8 with  the  first  row  of  Table  4.  This  and  all  subse- 
quent tables  will  be  drawn  in  italics  when  matrix  multiplication  is  in  force. 
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Table  4.  ON-AXIS  STRESS-STRAIN  RELATION  FOR 
UNIDIRECTIONAL  COMPOSITES  IN  TERMS 
OF  ENGINEERING  CONSTANTS 


OJ  CTg 

€. 

1 -Ell  0 

El  Et 

Ell  J-  0 

el  Et 

0 0 ± 

6 

Glt 

All  the  material  constants  of  the  stress- strain  relation  shown  in  this  table 
are  called  engineering  constants.  They  are  the  familiar  constants  used  for 
ordinary  materials  with  subscripts  added  to  denote  the  directionality  of 
properties.  Many  design  formulas  for  structural  elements  are  written  in 
terms  of  engineering  constants.  Thus  the  use  of  engineering  constants  will 
often  facilitate  the  use  of  composites  for  structural  applications.  This  con- 
cession to  the  state-of-the-art  design  methodology,  however,  can  lead  to 
an  unnecessarily  complicated  design  procedure.  In  fact,  engineering  con- 
stants for  composites  are  clumsy,  and  should  be  replaced  by  the  components 
of  compliance  and  modulus.  A change  of  notation  from  engineering  constants 
in  Table  4 to  components  of  compliance  in  Table  5 can  be  done  by  direct 
substitution. 


16 


* * 

- *+■  • 


>.*  / k \ 


* 


The  relations  between  these  two 

sets  of  elastic  constants  are: 
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From  Eq.  8 we  can  solve  for  stress  in  terms  of  strain  for  which 

we  have  the  following  equations: 
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where  m = 


f1  " *LT  ‘tl] 
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In  order  to  eliminate  the  clumsiness  of  engineering  constants  in  this 
stress- strain  relation,  we  will  introduce  components  of  modulus  in  Table  6, 


Table  6.  ON- AXIS  STRESS-STRAIN  RELATION  FOR  UNI- 
DIRECTIONAL COMPOSITES  IN  TERMS  OF  MODULUS 


On 

o,t 

0 

Ogt 

0M 

0 

0 

0 

Oee 

The  following  relations  exist  between  engineering  constants  and  the  compo- 
nents of  modulus : 


where  m 


«-» 

a 

°22 

el  ■ 

m 

•U  _ * — — 

T m 

is  s 

D 

pj 

00 

„ - 2li. 

LT 

°22 

YL  Qn 

glt 

°66 

L °12 

°21  1 

I"1 

1 - 


Q Q 
11  22 


(13) 


We  have  seen  three  sets  o£  material  constants,  any  of  which  can  com- 
pletely describe  the  stiffness  of  on-axis  unidirectional  composites.  The 
characteristics  of  each  set  is  summarized  in  the  following: 


• Modulus  is  used  to  calculate  the  stress  from  strain.  This  is 
the  basic  set  needed  for  the  stiffness  of  multidirectional 
laminates. 

• Compliance  is  used  to  calculate  the  strain  from  stress.  This 
is  the  set  needed  for  the  calculation  of  engineering  constants. 

• Engineering  constants  are  the  carryover  from  the  ordinary 
materials.  Old  designers  feel  more  comfortable  working  with 
the  engineering  constants. 

As  stated  earlier,  from  one  set  of  constants  we  can  readily  find  the  other 
sets.  They  are  all  equivalent.  There  is  a direct  relationship  between  the 
modulus  and  compliance.  One  is  the  inverse  of  the  other.  We  will  discuss 
the  process  of  inversion  later. 


4.  SYMMETRY  OF  COMPLIANCE  AND  MODULUS 

We  wish  to  show  that  the  coupling  components  of  compliance  and  those 
of  modulus  are  equal,  or  in  the  terminology  of  matrix  algebra,  that  the 
compliance  and  modulus  matrices  are  symmetric.  Since  the  only  coupling 
that  we  have  seen  thus  far  is  the  Poisson  coupling,  the  symmetry  condition 
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states  that  the  coupling  components  are  equal,  as  follows: 


(14) 


We  can  demonstrate  the  validity  of  these  equalities  from  the  stored  elastic 
energy  in  a body  subjected  to  stress  and  strain.  Let  the  stored  energy  at 
a point  in  the  orthotropic  body  be 


W ■ T [°i  61  + + °6  e6. 


From  this  definition  we  see 
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(16) 


Substituting  the  stress- strain  relation  in  terms  of  compliance  from  Table  5 


into  Eq.  15, 

W ■ 1 [S„  °i2  + 5 (s.2  + s2l)  + S22  °-22  + S66  <"}  <17> 

We  will  obtain  the  stress-strain  relation  by  differentiation  of  this  energy 
term  in  accordance  with  Eq.  16, 


(18) 

62  * i[s,2  + s2i]ai  +s22'r2 

Matching  the  like  constants  between  this  set  and  those  in  Table  5,  the  only 
condition  that  satisfies  both  sets  is 


6i  ■ sn°i  +2  [su  + s2i] 


S12  S21 


(19) 


By  substituting  the  modulus  relations  in  Table  6 into  Eq.  15  we  can  also 
•how  that 


(20) 
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The  last  two  equations  state  the  symmetry  conditions  of  the  Poisson  coupling. 
A similar  symmetry  condition  can  be  applied  to  engineering  constants. 

From  Eq.  9,  for  example,  we  have 


v _ E - v E 
LT  T TL  L 


(21) 


With  these  symmetry  conditions,  the  number  of  independent  constants  for 
the  on-axis  orthotropic  unidirectional  composite  are  reduced  by  one,  from 
five  to  four  in  Tables  4 to  6.  If  additional  symmetry  conditions  exist,  the 
number  of  constants  can  be  further  reduced.  Specifically,  two  such  cases 
exist: 


e Square  Symmetric  Materials 

If  the  longitudinal  and  transverse  properties  are  equal;  i.  e. , 


Q11  “ °22’  S11  " S22'  and  EL  = ET 


(22) 


we  have  a square  symmetric  material.  But  because  of  the 
additional  relation  in  Eq.  22,  the  number  of  independent  constants 
are  three,  one  less  than  the  orthotropic  material.  A cross-ply 
laminate  is  a square  symmetric  material  in  the  plane  of  the 
laminate. 


• Isotropic  Materials 

We  know  that  isotropic  materials  have  only  two  independent 
constants  because  there  is  another  relation  among  the  three 
remaining  constants;  i.  e.. 


2Q66  = Q11  ‘ °12 


S66  = 2(SU'S12) 


2 (1+*) 


(23) 
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This  relation  is  derived  from  the  equivalence  between  the  state 
of  pure  shear  and  that  of  equal  tension- compression.  This 
equivalence  is  only  valid  for  isotropic  materials.  The  deriva- 
tions of  these  relationships  will  be  discussed  later. 

In  summary,  the  stress- strain  relations  which  govern  the  stiffness  of 
all  materials  have  the  identical  form  for  unidirectional  composites  as  for 
ordinary  materials.  No  additional  terms  or  more  complex  relationships 
exist.  The  only  difference  is  the  number  of  independent  constants;  four  for 
composites  versus  two  for  ordinary  materials.  But  there  are  no  conceptual 
and  operational  barriers  that  would  make  composites  intrinsically  difficult 
to  work  with.  In  fact,  once  we  understand  composites,  we  automatically 
will  understand  ordinary  materials.  Ordinary  materials  can  be  treated  as 
special  cases  of  composites. 

5.  STIFFNESS  DATA  FOR  TYPICAL  UNIDIRECTIONAL  COMPOSITES 
a.  Measured  engineering  constants  for  a number  of  unidirectional 
composites  are  listed  in  Table  7.  The  fiber  volume  fraction  and  specific 
gravity  are  also  included. 


Table  7.  ENGINEERING  CONSTANTS,  FIBER  VOLUME 
AND  SPECIFIC  GRAVITY  OF  TYPICAL 
UNIDIRECTIONAL  COMPOSITES 


glt  vf 
GPa  % 


Sp. 

Grav, 


b.  The  compliance  components  for  the  same  composites  in  Table  7 
are  listed  in  Table  8.  These  components  are  computed  from  Table  7 using 


the  formulas  in  Eq.  9. 

Table  8.  COMPLIANCE  COMPONENTS  OF  TYPICAL 
UNIDIRECTIONAL  COMPOSITES  (TPa)'1 


Type 

sn 

S22 

S12 

S66 

T300/5208 

5.525 

97.09 

-1.  547 

139.  5 

B(4)5505 

4.902 

54.05 

-1.128 

172.  7 

AS/3501 

7.246 

111.6 

-2.174 

140.  8 

Scotchply 

/1002 

25.91 

120.  9 

-6. 744 

241.5 

Kevlar  49 
/Epoxy 

13.16 

181.8 

-4.474 

434.8 

24 


I 


c.  The  modulus  components  for  the  same  composites  are  listed  in 
Table  9.  These  components  are  calculated  using  the  formulas  in  Eq.  12. 


Table  9.  MODULUS  COMPONENTS  OF  TYPICAL 
UNIDIRECTIONAL  COMPOSITES  (GPa) 


T 


Type 

Mat'l 

m 

Q11 

°22 

°12 

°66 

T300/5208 

Gr/Ep 

1.  0045 

181.8 

10.34 

2.  897 

7.17 

B(4)/5505 

B/EP 

1.0048 

205.0 

18.  58 

4.  275 

5.79 

AS/3501 

Gr/Ep 

1.0059 

138.8 

9.013 

•2.  704 

7.1 

Scotchply 

/1002 

Gl/Ep 

1.0147 

39.16 

8.  392 

2.182 

4.14 

Kevlar  49 
/Epoxy 

Aramid 

/Ep 

1.0084 

76.64 

5.  546 

1.886 

2.3 
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a.  Find  Strain  from  Stress 
Given  stress  vector:  (400,60,15)  MPa 

For  compliance  of  T300/5208  from  Table  8: 


S..  = 5.  525  (TPa) 


Using  stress- strain  relation  in  terms  of  compliance,  such  as  that 
in  Table  5: 

Gj  = (5.  525  x 400  - 1.  547  x 60)  x 10-6 
= 2.  117  x 10"3 


If  a different  material  is  used,  we  only  need  to  replace  the  compliance 
components  in  Eq.  25  with  different  data.  If  the  new  material  is  Scotchply 
1002,  we  can  get  the  compliance  from  Table  8. 

S,  =25.91  (TPa)'1 


With  the  same  applied  stress  as  Eq.  24,  the  resulting  strain  is 
€,  =(25.91  x 400  - 6.744  x 60)  x 10-6  =9.959  x 10-3 


Since  the  glass  composite  is  less  stiff  than  the  graphite  composite, 
the  strain  produced  by  the  same  applied  stress  is  expected  to  be  larger  in  the 
glass  composite.  If  we  compare  the  strain  components  by  components  between 
Eq.  26  and  28,  the  strain  in  the  glass  composite  is  larger  in  two  components, 
and  smaller  in  one.  The  moral  of  the  story  is  that  biaxial  stress  and  strain 
states  are  complex.  Disciplined,  analytic  approach  is  straightfoward  and  is 
definitely  preferred  over  guesswork, 
b.  Find  Stress  from  Strain 

This  process  is  the  inverse  of  the  previous  process.  If  we  are 
given  the  strain  in  Eq.  26  and  apply  it  to  a T 300/ 5208  composite,  the  resulting 
stress  must  be  calculated  by  using 

• Stress- strain  relation  in  terms  of  modulus,  such  as  that 
in  Table  6,  and 

• Modulus  components  in  Table  9, 

Qn  = 181.8  GPa 

Q22  = 10.34  » 

Q12=2.897 

Q66=7-17 

The  resulting  stress  is: 

0^  = 181.8  x 2.  117  + 2.897  x 5.206  =400  MPa 
(T2  =2.897  x 2.  117  + 10.34  x 5.  206  =60  MPa  j (30) 

(T6  =7.  17  x 2.092  = 15  MPa  ) 


Note  that  the  original  stress  of  Eq.  24  has  been  recovered.  If  our  composite 
is  Scotchply  1002,  we  should  use  the  modulus  components  listed  in  Table  9 
for  this  material: 


Qn  =39.16  GPa 
Q22  =8.  392 
Q12=2.182 

Q66=4*14 

The  resulting  stress  from  the  applied  strain  in  Eq.  28  is: 
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0^  = 39. 16  X 9.  959  + 2.  182  x 4.  556  = 400  MPa 

(Tz  = 2. 182  x 9.  959  + 8.  392  x 4.  556  = 60  MPa 

(IT  = 4.14  x 3.623  = 15  MPa 
o 

Note  again  that  the  original  stress  of  Eq.  24  has  been  recovered. 


SECTION  II 

TRANSFORMATION  OF  STRESS  AND  STRAIN 
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SCOPE 

The  change  of  stiffness  of  unidirectional  composites  as  a function  of 
ply  orientation  is  a unique  feature  of  composites.  These  orientational  var- 
iations of  stress  and  strain  are  the  fundamental  underlying  issues  which 
must  be  understood.  We  will  derive  the  relations  that  govern  these  varia- 
tions; namely,  the  transformation  equations.  We  will  show  three  formula- 
tions for  the  transformation;  viz.,  the  conventional  power  functions,  the 
double  angle  functions,  and  the  invariant  functions.  Each  formulation  has 
its  unique  characteristic  and  is  useful  for  special  purposes.  All  three  formu- 
lations are  equivalent  and  will  yield  the  same  answer.  Two  key  relations 
for  stress  and  strain  transformations  in  double  angle  functions  are: 


PRINCIPAL  NOMENCLATURE 


first  order  invariant  of  stress  or  strain,  depending  on  the 

subscript. 

cos  0 

sin  6 

average  of  normal  stress  or  strain  components;  equal  to  first 
invariant  I. 

one  half  difference  between  normal  stress  or  strain  components, 
another  symbol  for  shear  stress  or  one  half  shear  strain;  this 
notation  is  used  in  conjunction  with  p and  q above, 
a second  order  invariant  of  stress  or  strain;  it  is  the  radius 
of  the  Mohr  circle, 
coordinate  axes 

new  or  transformed  coordinate  axes 
displacements  along  x and  y axes, 
displacements  along  x*  and  y'  axes. 

stress  components  in  material  symmetry  axes,  i = 1,2,6. 

1 

stress  components  in  1 ' — 2*  coordinate  system,  i-1,2,  6. 
principal  stress  components. 

strain  components  in  material  symmetry  axes,  i=l,2,  6. 

strain  components  in  1 • - 2*  coordinate  system,  1=1,2,  6. 

angle  of  ply  orientation;  counterclockwise 
rotation  is  positive. 

phase  angle  for  stress  or  strain  transformation  in  the  invariant 
formulation;  it  is  the  orientation  of  the  principal  axes, 
positive  and  negative  or  inverse  transformation  of  stress  or 
strain,  depending  on  the  subscript.  The  sign  corresponds  to 
that  of  the  ply  orientation. 
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1.  BACKGROUND 

Up  to  this  point,  we  have  only  dealt  with  the  stiffness  of  unidirectional 
composites  in  their  material  symmetry  axes,  as  shown  in  Fig.  14.  In  this 
reference  coordinate  system,  we  call  this  type  of  symmetry  orthotropic; 
see  Fig.  10  for  graphic  illustration.  General  orthotropic  configuration 
occurs  when  the  ply  orientation  is  different  from  0 or  90.  This  is  shown  in 
Fig.  15.  We  also  call  the  latter  configuration  the  off-axis  as  distinguished 
from  the  on-axis  in  Fig.  14. 


Fig.  14.  Material  symmetry  axes  of  Fig.  15.  Off-axis  or  generally 

orthotropic  configura- 
tion of  a unidirectional 
composite.  Counterclock- 
wise rotation  of  the  ply- 
orientation  is  positive; 
clockwise  rotation,  neg- 
ative. The  sign  of  angle 
of  rotation  is  critical. 

In  order  to  evaluate  the  properties  of  an  off-axis  unidirectional  com- 
posite; i.  e. , the  properties  in  l'-2'  coordinates  in  Fig.  15,  one  of  the  direct 
methods  is  to  transform  the  stress  from  the  off-axis  1 * - 2*  coordinates  to 
the  on-axis  1-2  coordinates.  Thus,  we  recover  the  orthotropic  symmetry 


a unidirectional  composite. 
The  1-axis  is  along  the  fiber 
and  is  in  the  longitudinal  dir- 
ection. This  on-axis  config- 
uration is  called  orthotropic. 


in  this  new,  transformed  coordinates,  and  can  now  apply  the  on-axis  stress- 
strain  relations  developed  in  Section  1,  such  as  Table  5.  We  can  determine  the 
induced  strain  from  the  applied  stress;  both  stress  and  strain  are  in  the  1-2 
system.  If  we  want  to  know  what  the  induced  strain  is  in  the  original  off- 
axis  1*  — Z*  system,  we  need  only  to  apply  the  inverse  transformation  to  the 
induced  strain  in  the  on-axis  1-2  system. 

The  sequence  of  operations  just  described  can  be  illustrated  in  Fig.  16. 

It  can  be  used  to  derive  the  off-axis  compliance  of  a unidirectional  com- 
posite. 

• The  originally  applied  stress  to  an  off-axis  composite  is  shown 
in  (a),  expressed  as  CT/. 

• If  we  apply  stress  transformation  to  the  components  of  (a)  in  1' -2' 
system,  we  will  get  (b),  the  same  state  of  stress  but  expressed  in  the 
different  components  of  the  on-axis  1-2  system;  i.  e.,  CT,  . 

This  positive  stress  transformation  is  designated  as  T^  . 

• Since  we  know  the  relations  between  stress  and  strain  from  Table 

5,  we  can  determine  the  induced  strain  in  the  on-axis  1-2  system,  €. 
which  is  shown  in  (c). 

• We  can  then  apply  inverse  strain  transformation  to  get  the  strain 
components  in  the  off-axis  l'-2'  system  from  the  on-axis  1-2 
system;  i.  e. , from  (c)  to  (d)  in  Fig.  16.  This  negative  strain 
transformation  is  designated  as  T . Then  we  have  the  induced 
strain  in  (d)  as  the  result  of  applied  stress  in  (a),  both  of  which 
are  in  the  off-axis  l'-2'  system. 
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Off-Axis  On -Axis  On -Axis  Off-Axis 

Stress  Stress  Strain  Strain 


Fig.  16.  Determination  of  the  off-axis  compliance: 

From  (a)  to  (b):  use  positive  stress  transformation 

From  (b)  to  (c):  use  the  on-axis  stress-strain  relation  in 
compliance. 

From  (c)  to  (d):  use  negative  strain  transformation 

We  can  go  from  (a)  to  (d)  directly  if  we  know  the  off-axis 
stress-strain  relation.  This  can  be  derived  by  merging  the  three  steps 
in  this  figure  into  one.  If  imposed  strain  is  given  in  Fig.  16(a)  instead  of 
stress,  the  induced  off-axis  stress  can  be  determined  by  a very  analogous 
method.  The  off-axis  modulus  of  a unidirectional  composite  can  be  derived 
from  this  process.  The  sequence  of  operations  is  illustrated  in  Fig.  17. 

• From  off-axis  strain  to  on-axis  strain,  use  positive  strain  trans- 
formation. This  is  the  operation  from  (a)  to  (b). 

• From  on-axis  strain  to  on-axis  stress,  use  the  on -axis  stress- 
strain  relation  in  modulus,  as  in  Table  6.  This  is  the  operation 
from  (b)  to  (c). 

• From  on-axis  stress  to  off-axis  stress,  use  negative  stress 
transformation.  This  operation  is  from  (c)  to  (d). 
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Off-  Axis 
Strain 


On  - Axis 
Strain 


On  - Axis 
Stress 


Off-Axis 

Stress 


Fig.  17.  Determination  of  the  off-axis  modulus: 

From  (a)  to  (b):  use  positive  strain  transformation. 

From  (b)  to  (c):  use  the  on-axis  stress-strain  relations  in 
modulus. 

From  (c)  to  (d):  use  negative  or  inverse  stress  transformation, 
Alternatively,  we  can  go  from  (a)  to  (d)  directly  if  we  know 
the  off-axis  modulus. 

The  scope  of  this  section  is  to  show  stress  and  strain  transformation. 
The  formulas  of  the  transformation  are  simple  and  easy  to  use,  but  the 
most  critical  part  of  the  operation  is  the  sign  convention.  As  we  have  re- 
peatedly mentioned,  signs  are  critical  for  the  study  of  composites.  Such 
emphasis  is  not  called  for  in  the  case  of  ordinary  materials  because  their 
behavior  is  often  insensitive  to  signs  and  directions. 

The  notations  associated  with  coordinate  transformation  are  arbitrary. 
The  components  of  the  original  versus  the  transformed,  the  old  versus  the 
new,  the  1 * — Z*  versus  the  1-2  systems,  or  the  primed  versus  the  unprimed 
are  based  on  a matter  of  judgment,  and  may  certainly  vary  from  author  to 
author  and  from  situation  to  situation.  Only  the  definition  of  the  on- axis 
and  the  off-axis  are  normally  fixed.  The  key  issue  is  the  initial  definition 
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such  as  that  shown  in  Fig.  15,  where  the  coordinates  and  the  ply  orientation 
are  illustrated.  This  definition  holds  only  for  the  present  situation  where 
we  are  concerned  with  unidirectional  composites.  When  we  study  the  lami- 
nated plates  a different  definition  will  be  used;  e.  g. , the  primed  and  unprimed 
coordinates  may  be  changed.  One  definition  that  we  will  keep,  however,  is 
that  the  counterclockwise  rotation  be  noted  as  positive. 

2.  TRANSFORMAT  ION  OF  STRESS 

Now  we  would  like  to  derive  the  relations  between  two  sets  of  stress 
components;  one  set  expressed  in  the  1 * - Z*  system,  and  the  other  in  1-2 

I 

system.  The  latter  is  rotated  from  the  former  by  a positive  angle  as  shown 
in  Fig.  15  and  repeated  in  Fig.  18(a).  In  Fig.  18(b)  & (c),  the  two  sets  of 
stress  components,  one  primed,  the  other  unprimed,  are  also  shown.  The 
arrows  indicate  the  direction  of  the  positive  component. 


Fig.  18.  Stress  transformation:  changes  in  stress  components  due  to 
coordinate  rotation  or  transformation. 

(a)  Relation  between  primed  and  unprimed  systems. 
Counterclockwise  rotation  is  positive. 

(b)  Primed, off-axis  stress  components. 

(c)  Unprimed, on-axis  stress  components. 

All  arrows  for  the  components  are  pointing  in  a positive 
direction. 
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The  transformation  of  stress  can  be  derived  from  the  balance  of 
forces.  Consider  a free-body  diagram  shown  in  Fig.  19(a)  which  is  a wedge 
sliced  across  fibers  in  a typical  infinitesimal  unit  area  like  that  in  Fig.  18(b). 
The  sides  of  this  wedge  have  the  following  lengths  relative  to  unity  hypotenuse; 
also  shown  in  Fig.  19(a):  1 


m = cos  8 
n * sin  6 


(33) 


The  forces  exerted  on  the  sides  of  this  triangular  free  body,  shown  only 
schematically  in  Fig.  19(b),  are  the  products  of  the  stress  components  mul- 
tiplied by  the  appropriate  lengths  of  the  sides. 

• Balance  of  horizontal  (along  the  l'-axis)  forces  yields: 


m(J^  - nO^  = mO"j 1 + nd^1 

(34) 

• Balance  of  vertical  (along  the  2' -axis)  forces  yields: 

n(TJ  + md^  = n(T^'  + mO^' 

(35) 

Keeping  in  mind  that  the  original,  l'-2'  components  are  given,  and  we  are  ' 
looking  for  the  new,  1-2  components.  We  can  find  the  unknown  components 

by  solving  simultaneously  Eq.  34  and  35,  as  follows: 

(Tj  = mVj’  + n20£’  + 2mn0^' 

(36) 

= -mnCT'  + mn<r;'  + (m2  - n2)  0^' 

(37) 
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(b) 


Fig.  19*  Free-body  diagram  for  the  balance  of  stress  components. 

The  components  of  on-axis,  1-2  coordinates  can  be  expressed 
in  terms  of  those  of  the  off-axis,  i'-2'  coordinates.  All 
stress  components  shown  are  positive. 

We  can  now  repeat  the  process  by  taking  a different  wedge  by  slicing 
parallel  to  fibers  in  the  unit  area  as  in  Fig.  18(b).  On  this  plane  the  normal 
stress  component  is  0^.  The  free  body  diagram  for  this  wedge  is  shown 
in  Fig.  20,  from  which  the  following  relations  can  be  established: 

• Balance  of  horizontal  forces  yields: 

n(T2  + m(T6  = n0^ ' - mO^'  (38) 

• Balance  of  vertical  forces  yields: 

m0"2  - nO^  = mO^'  - mO^'  (39) 
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Fig.  20.  Free-body  diagram  for  balance  of  stress  components. 


This  is  the  same  as  Fig.  19  except  the  new  plane  is  sliced 
along  the  fibers. 

If  we  solve  the  last  two  equations  simultaneously,  we  get 

CT  = n2(T'  + mV'  - 2mnQ-'  (40) 

2 12  o 

0^  = -mn(Tj'  + mnO^'  + (m2  - n )(J^'  (*l) 

Note  that  the  shear  stress  expressed  in  Eq.  41  is  the  same  as  that  in 

Eq.  37,  as  it  should  be.  Thus,  the  three  equations  for  stress  transformation 

are  Eqs.  36,  40  and  41.  These  equations  can  be  packaged  in  a matrix 

multiplication  table  as  we  did  in  Table  4,  as  follows: 

Table  10.  STRESS  TRANSFORMATION  EQUATIONS 
IN  POWER  FUNCTIONS 


The  transformation  equations  above  are  expressed  in  terms  of  second  power 
of  sines  and  cosines.  We  can  rewrite  these  equations  using  double  angle 
trigonometric  identities  as  follows: 


2 2 11 
m = cos  0 = — +—  cos29 

« L» 

2 2 11 
n = sin  9 = — - — cos29 

2mn  = sin29 
2 2 

m -n  = cos20 


When  we  substitute  these  identities  into  the  equations  in  Table  10,  we  get 
o;  = ^ ((T1'+  (T2')  + J (0^ ' - (r2')  cos20  + (T6'sin2e 

K,  = 2 ((ri'+  ^ ~ 2 (°l  ‘ °2>  cos20  ' (r6'8in2e  (43) 

1 I | | | 

= - — (CTj  - (T^  ) sin20  + (J^  cos20  ' 


Introducing  a notation  commonly  used  in  photoelasticity,  where 
p = |((r1+cr2),  q = |<<r1-(r2),  r = o"6 

or  p'  = ilTj'f  0-2'),  „•  = r = ^ 


(44) 


we  can  now  express  the  stress  transformation  equations  in  terms  of  double 
angles  and  the  notation  in  Eq.  44.  This  new  formulation  is  shown  in  a ma- 
trix multiplication  table  as  follows: 
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Table  11.  STRESS  TRANSFORMATION  IN  DOUBLE 
ANGLE  FUNCTIONS  - I 


cos20 
-cos20 
-sin 20 


sin20 

-sin20 

cos20 


0 is  positive  in  counter-clockwise 
rotation 

There  ia  an  alternative  arrangement  for  Table  11  where  the  column  headings 
and  the  trigonometric  functions  are  interchanged.  This  arrangement  is  use- 
ful for  certain  ply  orientation  such  as  45  degrees,  in  which  case  the  column 
with  the  cosine  function  vanishes. 

Table  12.  STRESS  TRANSFORMATION  IN  DOUBLE 
ANGLE  FUNCTIONS  - II 


sin20 


Either  table  can  be  used.  The  common  feature  is  the  first  column,  where 
the  influence  of  the  angle  of  rotation  does  not  exist.  The  constant  p'  is 
called  an  invariant  of  this  coordinate  transformation.  If  we  add  the  first 
two  rows  of  the  tables  above,  we  get 


2P  = o-j  + (r2  = 2P-  * (Tj1  + <r2 


(45) 


Thus  the  sum  of  the  two  normal  stress  components  remain  constant,  inde- 
pendent of  the  angle  of  rotation  or  ply  orientation.  We  call  this  invariant 
the  first-order  invariant  for  stress  transformation;  i.  e.. 


I = p = p'  (46) 

There  is  a second-order  invariant  that  we  can  show  as  follows: 

From  Eq.  44 

q2  + r2  =i((T-  (T,)2  + 0? 


From  Table  12 


4 '~l  ~2'  ' ”6 

2 2 2 2 

= q'  cos  20  + 2q'r' sin  28  cos  20  + r'  sin  20 

2 2 2 2 

+ r'  cos  20  - 2q,r'sin  20  cos  28  + qf  sin  20 


2 2 
= q'  + r' 


(47) 


This  is  another  invariant  because  the  quantity  remains  the  same  for  any 
value  of  angle  or  ply  orientation.  We  label  this  second-order  invariant  as 


_2  2 2 2 ,2 

R = q + r = q'  +r' 


(48) 


where  R is  the  radius  of  the  Mohr  circle  for  stress  transformation.  The 
geometric  relationship  of  Eq.  48  is  shown  in  Fig.  21.  Also  shown  in  this 
figure  are  the  phase  angle  0q  and  the  following  relations: 

q'  = Rcos  20 


r*  = Rsin  20 


28  .,a„-14..in-I4=co.-14 

o q'  R R 


(49) 
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Fig.  21.  Geometric  relatione  of  second-order  stress  invariant  +R 

and  Mohr  circle.  This  location  of  the  center  in  the  conven- 
tional stress  coordinates  is  specified  by  the  first  invariant  I . 
Thus  the  three  stress  components  that  characterize  the  state  of  plane  stress 
can  be  represented  by  at  least  three  sets  of  variables  for  a given  coordinate 
system;  i.  e. , for  a given  angle  of  orientation,  say,  0. 

• First  set:  (T. , 01,  (T.. 

ICO 

• Second  set:  p,  q,  r,  in  accordance  with  Eq.  44. 


• Third  set:  I,  R,  9 , in  accordance  with  Eq.  46,  48  and  49. 

o 

Similarly,  the  stress  transformation  can  be  formulated  in  terms  of  each  of 
the  sets  above.  We  have  done  the  first  two;  the  third  set  can  be  used  to  de- 
rive the  transformation  equations  by  substituting  Eq.  46,  48  and  49  into  the 
appropriate  column  headings  in  Table  11,  the  transformation  in  terms  of 

p'»  q'» 


* Ig_+  Rjy/cosZG^  cos26  + Rg*in20^  sin20 
= V+  Racos2<0  - 0o) 

85  V Vos2(e  - V 

= -Rgjcos20o  sin20  + R^sin20o  cos20 


(50) 


* -R(Jsin2(0  - ©e) 
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the  signs  are  correctly  built  in.  The  classical  power  functions  formulation 
in  Table  10  appears  most  frequently  in  existing  textbooks.  This  formulation 
is  most  convenient  to  use  when  one  or  two  of  the  stress  components  are  zero; 
e.g. , if  0^ ’ = <T6'  = 0. 

3.  NUMERICAL  EXAMPLES  OF  STRESS  TRANSFORMATION 
a.  Problem  #1:  Given  stress  in  1 * - 2*  system 

0-.'  = (2,4,6) 

Find  stress  components  in  1-2  system  for  0 = 45  degrees.  The 
two  reference  coordinates  are  shown  in  Fig.  22.  (Same  as  Fig.  18. ) 


Fig.  22.  Stress  transformation:  changes  in  stress  components  due  to 
coordinate  transformation. 

Solution: 

(1)  From  the  power  function  transformation  in  Table  10. 
cos  9 = sin  9 =V"2~ 


<T  = 2 (Z  + 4 + 2x6)  = 9 

CT2  = \ (2  + 4 - 2x6)  = -3 

<r6  = i (-2  + 4 + 0)  = 1 


(51) 
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(2)  From  the  double  angle  function  transformation  in  Table  12: 


cos  28  = 0 , sin  20  = 1 

p«  * \ (2  + 4)  = 3 | 

| (52) 

q'  = I (2  - 4)  = -1  | 

r'  = 6 

Then, 

Or  = 3+6  = 9 

<T2  = 3 - 6 = -3  j (53) 

°i  =1  ' 

(3)  From  the  invariant  function  transformation  in  Table  13: 

I = p*  =3 

R = l2  + 62  =V37 

0 = ^ tan  * -6  = -40.269  degree 

or  = -^sin”*  = 40.269  degree  (54) 

2 VtT 

or  = 4 cos ~ ^ - 49.731  degree 

2 V37 
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The  correct  answer  for  the  phase  angle  is  that  which  recovers  the 
original  given  stress  when  0 is  zero  in  stress  transformation  in  invariant 
functions  in  Table  13.  Thus,  using  this  table  and  the  invariants  in  Eq.  54, 
the  four  possible  states  of  stress,  each  corresponding  to  a phase  angle,  are 
shown  in  the  equivalent  Mohr  circle  in  Fig.  23.  Only  one  such  state  of 
stress  is  shown  in  Fig.  21. 


Fig.  23.  Four  possible  states  of  stress.  For  a given  Mohr  circle, 
the  sign  of  the  stress  components  depend  on  the  signs  of 
q'  and  r'  in  Eq.  52  (primes  are  omitted  in  this  figure). 

Each  state  of  stress  is  associated  with  a phase  angle.  The 
four  phase  angles  shown  in  this  figure  are  related  in  Eq.  56. 
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Note  that  the  last  phase  angle  is  the  correct  one  for  this  present 
example.  This  difficulty  of  the  invariant  formulation  is  caused  by  the 
multivalued  trignometric  functions.  Note  relationships  between  phase 
angles;  e.  g.. 


e(1) 

+ 

e(2)  = 90 

o 

e<3> 

+ 

e(4)  = -9C 

o 

! (56) 

e(1) 

_ 

-e(4) 

o 

° i 

e<2> 

o 

= 

-e<3> 

o 

(2) 

Using  the  correct  phase  angle  0^  and  the  given  angle  of  rotation  at  45 
degrees,  we  can  now  determine  the  transformed  stress  components  from 
Table  13: 

0^  = 3 +V37cos  2(45-  49.731)  = 9 

CT  = 3 -V37" cos  2(45-49.731)  = -3 
c* 

0*6  = -V^Tsin  2(45-49.731)  = 1 


Note  that  all  three  formulations  yield  the  same  answers,  as  we  can  see 
in  Eq.  51,  53,  and  57.  In  spite  of  the  multivalued  phase  angles  for  the  in- 
variant formulation  of  the  transformation,  the  phase  angle  has  one  impor- 
tant feature.  When  the  angle  of  rotation  6 is  equal  to  the  phase  angle  6^, 
from  Table  13  we  have: 

0^  = 3 +V37'  = 9.082 

<72  = 3 - VTf  * -3.082  (58) 

°6  S ° 


[• 

K ' 
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This  orientation  is  called  the  principal  direction.  For  in  this  state  of  stress, 
the  shear  stress  is  zero,  and  the  normal  stress  components  reach  maximum 
and  minimum  values.  They  can  be  determined  immediately  from  the  two 


invariants : 


T.  = <T  = I + R 
I max 


CTt  = (T  . = I - R 

II  mm 


where  stress  (T  and  CT  are  the  two  principal  stress  components.  See  Fig.  24. 


Fig.  24-  Principal  stress  components.  They  are  the  maximum  or 
minimum  values  in  the  Mohr  circle. 

Thus,  the  principal  direction  is  derived  from: 

0 - 9 =0  (60) 

o 

There  is  another  important  orientation;  i.  e. , 45  degrees  from  the  principal 
direction,  or  when 


At  this  angle,  we  have 


0^  = 1 = 3 


ff2  = I -3 


0 , - -R  « -vSf  = -6.082  I 

Here,  both  normal  stress  components  are  equal  to  the  first  invariant;  the 
shear  stress  component  reaches  its  minimum  value.  The  latter  would  have 
been  the  maximum  shear  stress  if  -45  is  used  in  Eq.  61. 

b.  Problem  #2:  Given  the  same  stress  as  Problems  #1,  i.  e. , 

0.'  = (2,4,6)  (63) 

But  the  direction  of  rotation  is  reversed;  i.  e. , 0 * -45;  see  Fig.  25.  This 
is  an  inverse  transformation  of  Problem  #1. 


1 

(b) 


Fig.  25.  Inverse  stress  transformation.  This  is  identical  to  that  in  Fig.  22 
except  the  angle  of  rotation  is  negative. 

Solution: 

(1)  From  the  power  function  transformation  in  Table  10: 
cos  0 = VZT  sin  0 = -y/i 


°1  = 2 (2  + 4 - 2x6)  = -3 


<^2  = g-  (2  + 4 + 2x6)  = 9 


°6  = 2 t2  " 4 + 0)  = -1 


, - v:  -— ;;  ' 


(2)  From  the  double  angle  transformation  in  Table  12: 


cos  20  = 0,  sin  20  = -1 


As  in  Problem  No.  1, 

P'  = 3,  q'  = -1,  r'  = 6 (65) 

Then  (Jj=3-6=-3 

r2  = 3 * 4 = ’ 

Note  that  when  the  angle  of  orientation  changes  sign,  the  normal  com- 
ponents are  interchanged,  and  the  shear  component  changes  sign.  Al- 
though we  did  not  show  the  transformation  by  the  invariant  formulation, 

the  quantities  in  Eq.  54  (I  , R , 0 ) remain  the  same. 

o 

c.  Problem  #3.  Given  a state  of  stress  in  Eq.  66;  l.  e. , 

0J  = (-3,9, -1)  (67) 

Find  the  stress  components  for  0 = +45.  This  is  the  inverse  transforma- 
tion of  Problem  #2;  see  Fig.  26.  This  transformation  is  the  same  as 
Problem  #1  except  the  original  state  of  stress  is  different. 


Fig.  26*  Inverse  stress  transformation  of  that  illustrated  in  Fig.  2 5. 
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Solution:  We  will  use  the  double-angle  formulation  here.  When  0 is 

+45  degrees 

cos  20  = 0,  sin  20  = +1 
From  Eq.  44 

p'  = I (-3  + 9)  = 3 = invariant  \ 

q'  = \ (-3  - 9)  = -6  (68) 

r'  = -1  ) 

From  Table  12 

<T  = 3 - l = 2 | 

<r2  = 3 + 1 = 4 j (69) 

*6  * 6 ' 

Note  that  we  have  recovered  the  original  stress  given  in  Problem  # 2. 

As  a final  emphasis  on  the  importance  of  the  sign  of  angles.  Fig.  27 

shows  the  consequence  of  a sign  error.  A positive  transformation  from  the 
1*  - 2'  axes  will  result  in  the  material  symmetry  axes,  designated 
1*  - 2+  axes.  A sign  error  will  result  in  the  1 - 2 axes  which  are  20 

orientation  away  from  the  correct  answer. 


Fig.  27.  Positive  and  negative  angles  of  rotation.  Guesswork  is  not 
good  enough  for  composites.  Keep  track  of  the  signs  ! 


4.  transformation  of  strain 

Strain  transformation  is  as  important  as  stress  transformation.  An 
identical  figure  to  Fig.  18  can  be  drawn  for  the  strain  components.  This  is 
done  in  Fig.  28. 


(a) 


Fig.  28  . Strain  transformation.  Changes  in  strain  components  due  to 
coordinate  rotation  or  transformation. 

(a)  Relation  between  primed  and  unprimed  systems. 
Counterclockwise  rotation  is  positive. 

(b)  Primed,  off -axis  strain  components. 

(c)  Unprimed, on- axis  strain  components. 

All  arrows  for  the  components  are  pointing  in  a positive 
direction. 

Like  strain  itself,  strain  transformation  is  purely  geometric  and  involves 
no  material  property  or  balance  of  forces.  Using  the  notation  shown  in 
Fig.  28,  the  off-axis  orientation  is  the  1* -2*  system,  and  the  on-axis,  the 
1-2  system.  We  will  now  derive  the  strain  transformation  relations  from 
the  strain- displacement  relations  shown  in  Eq.  1 and  4 and  repeated  as 
follows . 
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Since  both  displacements  u and  v and  coordinates  x and  y (in  place  of 


1 and  2)  are  vectors,  and  are  directionally  dependent  quantities,  we  only 
need  to  find  the  relationship  between  the  primed  and  the  unprimed  com- 
ponents of  a vector. 

First,  from  analytic  geometry,  we  can  derive  the  following  rela- 
tions between  the  primed  and  unprimed  coordinate  systems  as  shown  in 
Fig.  29  (a)  and  (b),  respectively. 


x = 

mx' 

+ 

ny1 

y = 

-nx' 

+ 

my' 

conversely. 

X1  = 

mx 

- 

ny 

y*  = 

nx 

+ 

my 

where,  as  before. 

m = 

C08 

6. 

n = 

53 


. . fe  w.  I 


*< 

•v 


/ 


(b) 


Fig.  29.  Coordinate  systems  between  the  primed  and  unprimed  axes. 

(a)  To  go  from  primed  to  unprimed,  0 is  positive. 

(b)  To  go  from  unprimed  to  primed,  9 is  negative. 

From  Eq.  72,  we  can  get  the  following  by  partial  differentiation: 


bx  &x’  by' 

sr  =m»^  = -n'-sr  = 


(73) 


The  relations  between  displacements  in  the  primed  and  unprimed  coordinates 
are  identical  to  those  in  Eq.  71  and  72  because  all  quantities  are  vectors.  We 
can  simply  write  the  following  by  replacing  x,  y,  x’,  y1,  by  u,  v,  u1 , v1 , 
respectively: 


u = mu'  + nv' 
v = -nv'  + mv' 

conversely,  u'  = mu  - nv 

v'  = nu  -f  mv 


(74) 

(75) 


Now  we  are  ready  to  derive  the  strain  transformation  equations. 
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From  Eq.  70 


by  chain  differentiation 

bu'  bx]  bu  by* 
3x‘  bx  by1  Sx 


from  Eq.  73  and  74 


€ 


1 


where 


6. 

i 


r 

bu' 

bv1 

r 

bu' 

bv'  "I 

r 

+ 

nSP 

m + j 

m 7y  i 

+ n 

Sy'J 

2 

m 

2 « 
n €2 

+ mn€ 

« 

6 

(76) 

bu' 

bv' 

bu' 

bv1 

bx' 

by' 

: fe6  = 

3?  + 

3? 

Here  primes  are  added  to  all  quantities  in  Eq.  70.  We  can  do  this  because 
the  relationship  is  invariant;  i.  e. , the  relationship  does  not  change  from 
coordinates  to  coordinates,  and  is  valid  for  all  coordinate  systems.  Note 
that  the  strain  transformation  in  Eq.  76  is  very  similar  to  the  stress  trans- 
formation in  Eq.  36  except  the  factor  2 is  missing  in  the  shear  term. 

This  difference  comes  about  from  the  use  of  engineering  shear  strain  as 
shown  in  Table  3 and  Eq.  4. 


By  an  identical  process  as  that  used  in  the  derivation  of  Eq.  76,  we  can 
show 

€2  = n26J  + m26^  - mnfi^  (77) 

6,  = -2mn6*  + 2mn6'  + (m2-n2)€!  (78) 

6 1 £ D 


This  is  summarized  in  a matrix  multiplication  table  as  follows: 


We  can  express  the  transformation  relations  in  terms  of  double  angle  and 
invariant  functions  like  we  did  for  the  stress  transformation.  Comparable 
to  Tables  11  and  12  for  stress  transformation,  we  have  for  strain  trans- 
formation in  double  angle  functions  shown  in  Tables  15  and  16,  respectively. 


Table  15.  STRAIN  TRANSFORMATION  IN 

DOUBLE  ANGLE  FUNCTIONS  - I 


Table  16.  STRAIN  TRANSFORMATION  IN 

DOUBLE  ANGLE  FUNCTIONS  - II 


where  p'  = I (€J  + €^),  q«  = j (€'  - €p, 

Note  that  the  definition  of  r'  is  different  from  that  for  the  stress 
transformation  in  Eq.  44.  The  use  of  engineering  shear  strain  is  respon- 
sible for  the  difference.  / 

! 

The  invariant  function  comparable  to  Table  13  for  the  stress  trans 
formation  can  be  derived  in  a similar  fashion  and  the  results  are  listed  in 
a matrix  multiplication  table  as  follows: 


Table  17.  STRAIN  TRANSFORMATION  IN 
INVARIANT  FUNCTIONS 


where:  I - I_  = 1(6.*  + 6 ') 

c 2 1 6 

R - V^'2  + r'2  =Vi(e,'  - s2')2  + i«6'2 

26  ~ cos  1 = sin  1 ^ = tan  1 — 

o R R q 

The  advantages  and  disadvantages  of  each  formulation  for  the  strain  trans- 
formation are  similar  to  those  for  the  stress  transformation.  The  double 
angle  formulation  appears  to  provide  the  best  compromise  and  is  recom- 
mended for  general  usage.  This  will  be  our  choice  for  the  balance  of  this 
book. 
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5.  NUMERICAL  EXAMPLES  OF  STRAIN  TRANSFORMATION 


a.  Problem  I.  Given  a state  of  strain  in  I*  - 2*  system 

= (2,4,6) 

Find  transformed  strain  for  6 equal  to  45  degrees.  See  Fig.  30. 
Solution:  From  Table  16,  cos  26  = 0,  sin  20  = 1 

p'  = (2  + 4)/2  = 3 

q'  = (2  - 4)/2  = -1 

r'  = 6/2  = 3 

Then  G^  = 3 + 3 = 6 

*2  = 3 - 3 = 0 

G6  = -2  x(-l)  = 2 


(81) 


(82) 


Note  that  the  transformed  strain  is  quite  different  from  the  trans- 
formed 8 tress  of  (9,  -3,1)  from  Eq.  53.  The  factor  of  2 in  the 
engineering  shear  strain  is  responsible  for  this  dramatic  difference. 

b.  Problem  2.  Try  inverse  transformation;  for  strain  in  Eq.  82 

€!  = (6,0,2) 

Find  strain  at  -45  degree  rotation. 


Solution:  From  Table  16,  cos  26  = 0,  sin  26  = -1 

p'  = (6  + 0)/2  = 3 

q'  = (6  - 0)/2  = 3 (83) 

r'  =2/2  = 1 
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Fig.  30.  Strain  transformation.  To  go  from  (b)  to  (c)  is  a positive 
transformation  when  the  angle  of  rotation  is  positive,  as 
in  Problem  #1.  The  negative  or  inverse  transformation 
in  Problem  #2  is  in  effect  when  the  angle  is  negative. 
(This  figure  is  the  same  as  Fig.  28.  ) 


Then 


Note  that  the  original  strain  is  recovered. 


c.  Problem  3:  Strain  rosettes  record  normal  strain  components  in 
several  directions.  Because  there  are  three  strain  components  in  a given 
2-dimensional  state  of  stress,  we  need  a minimum  of  three  independent 
normal  strain  readings  in  order  to  define  the  strain  state.  Thus,  any  three  or 
more  orientations  of  the  strain  rosette  can  completely  determine  the  desired 
strain  state.  Four- element  rosettes  spaced  at  45  degrees  apart  are  often 
used.  This  arrangement  provides  a redundant  element  for  crosschecking, 
and  insuring  against  a defective  or  malfunctioning  element.  The  strain 
transformation  equations  of  Table  16  are  needed  for  analyzing  the  recorded 
normal  strain  components. 


Solution:  Using  superscript  to  denote  rosette  orientation, 
we  have  from  the  first  row  of  Table  16. 

0 

€ = p + q 

_90 


= p - q 

45 

6 = p + r 

p’45 

e = p - r 


(85) 


There  are  four  equations  for  three  unknowns,  p,  q,  r,  such  that 


p = (€°  + €9°)/2  = (€45  + €"45)/2 


(S' 


(6 


45 


€9°)/2 

6*45)/2 


(86) 


Once  we  know  these  values,  we  can  place  them  in  appropriate 
columns  in  Table  15  from  which  we  can  calculate  the  strain  com- 
ponents for  any  orientation.  Any  one  of  the  four  elements  can 
fail  without  affecting  the  uniqueness  of  the  solution  in  Eq.  86. 

If  for  example  the  -45  degree  element  failed,  the  relationship 
in  the  last  equation  must  be  altered  as  follows: 


r 


(87) 
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SECTION  III 

OFF-AXIS  STIFFNESS  OF  UNIDIRECTIONAL  COMPOSITES 


SCOPE 


The  stiffness  of  unidirectional  composites  with  off-axis  ply  orientation 
are  important  because  composite  laminates  normally  consist  of  many  off- 
axis  plies  in  addition  to  on-axis  plies.  We  must  know  how  to  determine  the 
contribution  to  the  laminate  stiffness  by  each  ply  or  ply  assembly.  We 
will  need  the  transformation  of  modulus  and  compliance  to  determine  the 
off-axis  stiffness.  These  transformation  relations  can  be  formulated  in 
terms  of  power  functions,  the  multiple  angle  functions  and  the  invariants. 
Specific  examples  of  a graphite- epoxy  composite  are  used  to  enhance  under- 
standing of  the  off-axis  stiffness  of  unidirectional  composites.  The  key 
relation  of  this  and  subsequent  sections  is  the  modulus  transformation 
equation.  The  multiple  angle  formulation  of  the  transformation  is: 
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PRINCIPAL  NOMENCLATURE 

E.  , E_  = Young's  modulus  along  and  transverse  to  an  off-axis  direction, 

X U 

respectively. 

I , I = Linear  or  first  order  invariants  of  modulus  or  compliance. 

1 « 

m = cos  0 

n = sin  0 

Q..  = Components  of  modulus;  i,  j = 1,2,6. 

R ,R  = Quadratic  or  second  order  invariants  of  modulus  or  compliance. 
1 £ 

S..  = Components  of  compliance;  i,  j = 1,2,6. 

U = Linear  combinations  of  modulus  and  compliance  in  multiple- 

i 

angle  formulation.  Same  notation  but  different  definition  is 
used  for  modulus  and  compliance. 

07  = Components  of  stress;  i = 1,2,6. 

O’.'  = Transformed  components  of  stress;  i = 1,2,6. 

G.  = Components  of  strain;  i = 1,2,6. 

€.'  = Transformed  components  of  strain;  i = 1,2,6. 

v = Major  Poisson's  ratio  along  an  off-axis  orientation. 

1 £ 

A = Determinant  of  modulus  or  compliance. 


1 . OFF  - AXIS  MODUL  US 

As  we  have  shown  in  Fig.  17  and  repeated  here  in  Fig.  31,  the  off- 
axis  modulus  can  be  determined  in  three  steps:  the  off-axis  to  on-axis  strain 
transformation,  the  on-axis  stress-strain  relations,  and  the  inverse  or  the 
on-axis  to  off-axis  stress  transformation.  This  process  was  initiated  by  a 
given  strain  in  Fig.  31(a)  and  led  us  eventually  to  the  induced  stress  in  Fig. 
31(d).  The  off-axis  compliance  can  be  similarly  derived  in  three  steps, 
as  shown  in  Fig.  16.  The  purpose  here  is  to  derive  the  off-axis  modulus  and 
the  off-axis  stress- strain  relation  for  this  arbitrary  angle  of  orientation.  Then 
we  can  go  directly  from  (a)  to  (d)  in  Fig.  31  in  one  step. 


Off-  Axis  On  - Axis  On  - Axis  Off-Axis 

Strain  Strain  Stress  Stress 


Fig.  31*  Determination  of  the  off-axis  modulus: 

From  (a)  to  (b):  use  positive  strain  transformation. 

From  (b)  to  (c):  use  the  on-axis  stress- strain  relations  in 
modulus. 

From  (c)  to  (d):  use  negative  or  inverse  stress  transformation. 
We  caui  go  from  (a)  to  (d)  directly  if  we  merge  these  three 
steps  into  one.  This  is  the  same  as  Fig.  17. 
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We  will  follow  these  steps  in  Fig.  31. 

• To  go  from  (a)  to  (b),  we  need  the  strain  transformation  listed 
in  Table  14,  repeated  here  as  follows: 

6,  = m^6.'  + n26  ' + mn6  ' 
114o  . 

G = n26  ' + m26,'  - mn6  ' l (88) 

2 12  6 I 

2 2 / 

G , = -2mn6  ' + 2mn6  ' + (m  -n  )€,  1 
6 12  6 

• To  go  from  (b)  to  (c)  in  Fig.  31,  we  need  the  on -axis,  orthotropic 
stress -strain  relation  in  modulus  in  Table  6,  which  when  combined 
with  the  results  in  Eq.  88  produces: 

0*1  = Qj  1(m2G1’  + n262'  + mnG^) 

+ Q12(n261'  + m262'  - mn©6') 

= (m2Qn  + + (n2QH  + 

+ (mnQ^  - mnQl2)€6'  (89) 

Similarly,  = (m2Q12  f n Q22^l'  + ^12  + m ^22^2* 

+ (mnQj2  - mnQ22)6^'  (90) 

(T6  = + 2™Q6662'  + (m2-n2)Q66e6‘  (91) 

• To  go  from  (c)  to  (d)  in  Fig.  31,  we  need  inverse  or  negative 

stress  transformation  as  listed  in  Table  10.  The  angle  of  rota- 
tion is  now  negative,  however,  and  the  primed  and  unprlmed  are 
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interchanged.  The  unprimed  is  now  the  old  (before  transformation),  and 
the  primed  is  now  the  new  (after  transformation). 

(Tj*  = m2^  + n2<J^  - 2mn(T6 

m2  [(m  Q11  + n Q12)€l'  + (n  Q11  + m Q12)€2' 

+ (mnQu  -mnQ12)66'j 

+ -2  [""2qu  + ”2q22,61’  + <“2q12  + m2Q22,e2' 

+ (mnQ12  - mnQ22)66'j 

-2mn  [-2n»Q66ei'  + ZmnQ^'  + (m2-n2,Q ^6^ 

[m  Qll  + “ Q22  + 2m  n Q12  + 4m  n Q66]  *1* 

+ [m2n2Qn  + m2n2Q22  + (m4+n4)Q12  - 4m2n2Q66]  62' 


r 3 3 33  33* 

+ Im  nQjj  - mn  Q22  + (mn  -m  n)Q12  + 2(mn  -m  n)Q^ 

66l 

<7  - Qll'6l'+Q12'e2'tQ16,e6' 

(92) 

Similarly,  0 ’■  = + Q22'  G2'  t 0^6^ 

(93) 

V ■ Q6l’Sl'  ♦ V S2'  + 

1 

(94) 

This  is  the  off-axis  stress -strain  relation  that  directly  relates  the  given 
•train  in  Fig.  31  (a)  to  the  resulting  stress  in  31(d),  and  redrawn  in  Fig.  32. 
This  relation  can  also  be  arranged  in  a matrix  multiplication  table  as  follows: 
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Of f -A  xis  Off  - A xis 

Strain  Stress 


Fig.  32.  The  off-axis  stress- strain  relations  in  modulus.  We  have 

merged  the  three  steps  in  Fig.  31  into  one.  The  relationship 
remains  valid  if  all  primes  are  omitted  in  this  figure.  This 
new  unprimed  applies  to  the  coordinates  of  Fig.  33. 

Table  18.  OFF-AXIS  STRESS-STRAIN  RELATION  FOR 


UNIDIRECTIONAL  COMPOSITES  IN  TERMS 
OF  MODULUS 


Since  the  choice  of  the  primed  and  unprimed  notation  is  arbitrary,  this 
table  remains  valid  if  all  primes  are  eliminated.  Then  the  relation  applies 
to  the  unprimed  coordinates  shown  in  Fig.  33. 


2 


Fig.  33.  A new,  unprimed  coordinate  system.  Unprimed  coordinates 
are  now  the  off-axis  configuration.  This  is  different  from 
Fig.  15.  This  coordinate  system  is  used  for  multidirectional 
laminates. 

The  major  difference  between  the  on-axis  stress -strain  relation  in  Table  6 
and  the  off-axis  relation  in  Table  18  lies  in  the  additional  components  in  the 
modulus.  These  subscripts  16  and  26  components  are  shear  coupling  terms 
that  relate  the  shear  strain  to  normal  stress,  or  normal  strain  to  shear 
stress.  Such  coupling  does  not  exist  in  ordinary  materials,  or  in  unidirect- 
ional composites  in  their  on-axis  orientation.  Geometric  illustration  of  the 
shear  coupling  effects  will  be  done  when  we  develop  the  generally  orthotropic 
compliance.  Symmetry  of  these  components  can  also  be  demonstrated  in  a 
manner  similar  to  that  used  for  the  component  in  Section  1.  The  stored 

energy  in  Eq.  17  must  contain  interaction  terms  of  (T  O",  and  (T,  <J ' , or  their 

16  £.  o 

equivalent  strain  components  € 6 and  €,6, . 

Id  2 o 

The  relationship  between  the  modulus  components  of  the  on-axis  and 
the  off-axis  orientations  can  be  summarized  in  Table  19,  where  matrix 
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multiplication  is  implied.  These  relations  result  from  the  derivation  of  Eq.  92 
Sind  what  was  omitted  in  Eq.  93  and  94.  These  relations  are  limited  to  trans- 
formation from  the  on-axis,  orthotropic  orientation  where  shear  coupling 
components  are  zero.  Note  that  and  no*  aPPear  a8  column 

headings  in  this  table.  * 

Table  19-  TRANSFORMATION  OF  MODULUS  FROM  ON -AXIS 

UNIDIRECTIONAL  COMPOSITES  IN  POWER  FUNCTIONS 


On 

o„ 

0/2 

Oee 

On 

m4 

n4 

2m*n* 

4m2n* 

o'„ 

n 4 

m4 

2m*n* 

4m*n* 

o;* 

m2n2 

m2n2 

m4-hn 4 

- 4m* n* 

O'ss 

m2n 2 

— 2n2 
m n 

-2m2n2 

m2-n2 

Ole 

m3n 

-mn3 

mn3-m3n 

2(mn  3 - m3n) 

Ole 

mn3 

-m3n 

m3n-  mn3 

2(m3n-  mn3) 

m - cos  0, 

n = sin  0 

Note  that  all  the  sums  of  exponents  of  the  trigonometric  functions  in  this 
table  are  in  the  fourth  power  which  are  by  definition  characteristic  of  the 
4th  rank  tensor.  The  stress  transformation  equations  are  governed  by  the 
2nd  power  functions,  as  we  have  shown  in  Table  10,  and  belong  to  the  2nd 
rank  tensor.  The  strain  transformation  equations  in  Table  14  are  also 
governed  by  2nd  power  functions,  but  are  different  from  those  for  stress 
because  engineering  shear  strain  is  used  which  is  twice  the  tensorial 
shear  strain. 

The  critical  issue  again  is  the  sign  convention.  The  angle  used  in  this 
table  is  the  ply  orientation.  Because  of  its  importance.  Fig.  15  and  33  are 

♦If  the  transformation  is  from  one  off-axis  orientation  to  another,  additional 
columns  for  Qjg  and  must  be  present. 
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shown  here  again  for  emphasis.  Although  the  primed  versus  unprimed 
coordinates  may  be  interchanged  from  situation  to  situation,  it  is  imperative 
that  the  sign  convention  is  clearly  understood.  For  unidirectional  composites, 
the  on-axis,  orthotropic,  and  material  symmetry  axes  coincide.  We  normally 
use  the  unprimed  axes  to  denote  this  configuration.  The  off-axis,  generally 
orthotropic  configuration  refers  to  ply  orientations  other  than  0 or  90  degrees. 
We  normally  use  the  primed  axes  for  the  off-axis  situation.  This  is  shown  in 
Fig.  34(a).  But  for  multidirectional  laminates,  there  can  be  many  ply  orienta- 
tions. So  a new  system  of  coordinates  is  more  convenient.  The  1-2  axes  are 
some  reference  coordinates  for  the  laminate.  Each  ply  orientation  0.  can  be 
designated  by  1.-2.  axes.  The  angle  used  in  Table  19  is  that  shown  in  Fig.  34(a). 
This  sign  convention  is  not  used  universally.  Some  authors  define  ply  orienta- 
tion opposite  to  that  shown  in  Fig.  34(a).  Then  their  transformation  relations 
will  be  different  from  those  shown  in  Table  19.  In  particular,  any  term  that 

3 

has  tiie  odd  power  of  sines  (n  and  n ) must  change  its  sign.  Only  the  shear 
coupling  terms  and  Q^are  affected. 


v< 


Fig.  34.  Positive  ply  orientation  is  shown.  The  notation  for  uni- 
directional composites  normally  follows  that  in  (a);  that  for 
multidirectional  composites,  in  (b)  where  0.  is  the  orientation 
of  the  i-th  ply  or  ply  assembly.  Primes  have  been  deleted. 


We  can  further  develop  a multiple -angle  formulation  for  the  modulus 
transformation  in  place  of  the  power  functions  in  Table  19.  This  process 


,>\*t  y >.v  V*  • 


I 


^ 

can  be  done  directly  by  substituting  the  following  trigonometric  Identities 
into  Table  19: 

m^  = ~ (3  + 4cos29  + cos46) 

O 

m^n  = ^-(2sin20  + sin40) 

2 2 1 

m n = ~ (1  - cos49) 

O 

3 1 

mn  = — (2sin26  - sin40) 

n*  = ”"(3  - 4cos20  + cos49) 

O 

| 

We  will  now  show  three  examples  of  the  substitution  of  values  from  Eq.  95 
into  Table  19: 

Q'n  = mQjj  + n Q22  + 2m  n (QJ2  + 20^) 

[ 

= ~(3  + 4cos29  + cos49)  + |-(3  - 4cos20  + cos40)Q22 
+ i(l  -co.46)  <Q12  + 2066> 

- £<*»„  4 3Q22  4 2°l2  4 4°66>  + 2 (Q11  ' Q22)co'29 

+ ?(QlltQ22-  2QU  - 4Q66)co849 
= Uj  + U2cos20  + U3cos40  (96) 

q'i2  - m2"2  <0,1  + q22  - 4tw  4 <m4+"4>Q12 

■ ?(*  - ‘<”49>  ,Qn  +q22  ’ “it1  4 8 4 2“*49>°i2 

■ 8 ,Qn  4 °22  4 6QU  - ‘V  * 8 <Q11  4 °22  ■ 2Q12  • 4Q66,CO,4e 


= U4  - U3cos40  (97) 


(95) 


• W- 


From  the  transformation  equations  in  Table  20,  we  can  show  the  off-axis 
combinations  listed  in  Eq.  99  are: 

Ui->Q'll  + 3Q22  + 2Q'l2  + 4Q66) 

From  Table  20: 

■ !<6U1  4 2U4  4 4U5> 

From  Eq.  99: 

■ 64  [(18QU  4 I8Q22  + l2Ql2  4 24Q66>  + (2Q11  4 2Q22  4 

12Q12  - 8Q66>  + (4Q1.  + 4Q22  - 8Q12  + 16Q66)J 

‘ &24°ll  4 24Q22  4 ,6Q12  4 32Q66J 

■it3Qll43Q22  + 2Ql2  + 4Q66] 

= Uj  (100) 


Similarly 


U'  * U cos20 

W b 


U'j  = U^cos40 


u4 


U5=  U5 


(101) 


When  the  off-axis  and  on-axis  combinations  are  equal,  they  are  by  definition 

invariant.  This  is  analogous  to  the  stress  invariants  in  Eq.  46.  Note  that 

U , U , and  U.  are  first-order  or  linear  invariants,  of  which  two  are  independent 
12  3 

because  we  can  show  from  Eq.  99  that 

2Ug  * - U4  (102) 

We  have  shown  that  stress  and  strain  possess  a second-order  or  quadratic 
invariant  each;  i.  e. , the  radius  of  Mohr  circle.  The  modulus  also  possesses 
second-order  invariants.  They  can  be  derived  as  follows.  We  will  first 
define  two  additional  linear  combinations  for  the  modulus. 
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■N 


(103) 


U6  * i,Q'l6  + Q26>  ■ 7U2,in2e 

U7  ‘ T(Q'i6  - Qk>  - V1”49 


We  can  now  derive  two  second-order  invariants  as  follows: 

2 2 '2 
Rj  = + 4Ug 

*=  U2(co8220  + sin220) 


= U! 


or 


R1  = ±U2 


(104) 


Similarly, 


_2  TTt2  _t12 
r2  . U'5  + w7 

= U2(co8240  + sin240) 

- u2 


or 


R2=±U3 


(105) 


where  R and  R are  invariants  and  U and  U are  not,  the  values  of  U 
12  Z 3 Z 

and  stated  in  Eq.  99  are  based  on  on-axis  orientation  for  which  the  shear 
coupling  terms  are  zero.  The  R' s are  radii  of  the  equivalent  Mohr  circle 
for  the  stress.  They  must  always  be  positive.  The  U.  can  be  positive  or 
negative.  In  fact,  would  be  negative  if  the  longitudinal  and  transverse 
directions  are  interchanged.  The  correct  sign  in  Eq.  104  and  105  must 
be  picked  to  make  R^  and  R^  positive. 


We  can  derive  the  transformation  equations  for  modulus  in  terms  of 

invariants  like  we  did  for  stress  and  strain  transformations  and  shown  in 

Tables  13  and  17.  If  we  limit  ourselves  to  transformed  modulus  from  the 

material  symmetry  or  orthotropic  axes,  the  results  will  be  identical  to 

those  shown  in  Table  20;  in  which  case,  positive  values  of  U and  U are 

2 3 

assigned  to  R^  and  R^»  respectively.  This  is  done  in  Table  21,  where  matrix 
multiplication  is  implied. 
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Table  21 . TRANSFORMED  MODULUS  FROM  ON-AXIS  UNIDIRECTIONAL 
COMPOSITES  IN  INVARIANT  FUNCTIONS 


1 

R, 

R* 

0), 

u, 

cos20 

cos40 

0** 

u, 

- cos20 

cos40 

o\* 

u4 

- cos40 

Oee 

U5 

- cos40 

O', 6 

Y sin 20 

sin40 

0*6 

• j sin20 

- si n40 

This  table  is  valid  for  the  I-axis  to  be  pointed  along  the  fiber  orientation, 
like  that  in  Fig.  14  or  34(a). 


2.  EXAMPLES  OF  OFF-AXIS  MODULUS 

We  will  show  in  this  section  the  transformed  modulus  for  a graphite- 
epoxy  composite.  The  particular  material  system  that  we  are  going  to  use 
is  the  Union  Carbide  and  Toray  T300  filament  and  Narmco  5208  resin,  or 
T300/5208  for  short.  The  modulus  data  for  this  material  were  listed  in  Table  9. 
We  can  immediately  calculate  the  transformed  modulus  by  substituting  the 
data  into  the  transformation  equations  in  Tabic  19.  Numerical  data  for  the 
transformed  modulus  are  listed  in  Table  22  and  plotted  in  Fig.  35.  All  six 
transformed  components  of  the  modulus  are  shown.  The  angle  of  ply 
orientation  is  also  shown,  where  counterclockwise  direction  is  positive. 
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Table  22  . TRANSFORMED  MODULUS  OF  T300/5208 
UNIDIRECTIONAL  COMPOSITES  (GPa) 


e 

Q' 

11 

Q' 

22 

Hz 

Q66 

Qi6 

°26 

m 

181.8 

10.3 

2.90 

7.17 

0 

0 

15 

160.4 

11.9 

12.75 

17.05 

38.50 

4.36 

30 

109.3 

23.6 

32.46 

36.78 

54.19 

20.05 

45 

56.6 

56.6 

42.  32 

46.59 

42.87 

42.87 

60 

23.6 

109.3 

32.46 

36.78 

20.05 

54.19 

75 

11.9 

160.4 

12.  75 

17.05 

4.  36 

38.  50 

90 

10.3 

181.8 

2.90 

7.17 

0 

0 

An  alternative  method  of  arriving  at  the  same  transformed  modulus 
is  the  use  of  the  multiple -angle  or  the  invariant  formulation  in  Table  20 
or  21,  respectively.  We  must  first  determine  the  values  of  the  U^  using  Eq. 
99  and  data  in  Table  9.  The  results  of  several  unidirectional  composites 
including  T300/5208  are  listed  in  Table  23.  A typical  calculation  is  listed 
as  follows: 

From  Eq.  99  U,  - |oQ]  , + 30^  + 2Q[2  + 4Q66> 

From  Table  9 U = -|(3  x 18 1 . 8 + 3 x 10.  34  + 2 x 2.  897  + 4 x 7.  17) 

= 76.  37  GPa  (106) 


Table  23.  LINEAR  COMBINATIONS  OF  MODULUS  FOR 
TRANSFORMATION  OF  MODULUS  (GPa) 


i 1 

U1 

U2 

U3 

U4 

U5 

1 

T300/5208 

76.37 

85.73 

19.71 

22.61 

26.88 

B(4)/5505 

87.80 

93.21 

23.98 

28.26 

29.77 

AS/ 3 501 

59.66 

73.90 

14.25 

16.95 

21.35 

Scotchply 

/1002 

20.47 

15.4 

3.  33 

5.  53 

7.47 

Kevlar  49 
/Epoxy 

32.44 

35.55 

8.65 

10.54 

10.95 

With  the  values  listed  in  this  table  and  the  equations  in  Table  20  or  21,  we 
can  arrive  at  the  same  transformed  modulus  listed  in  Table  22  and  shown 
in  Fig.  3b. 

A typical  calculation  by  the  multiple- angle  transformation  is  listed 
as  follows: 

From  Table  21:  Q'^  = U^  + ^^820  + U^cos40 

From  Table  23  for  T300/5208: 

Q'  - 76.  37  + 85.73cos20+  19.71cos40 

when  0=  45  degrees,  cos20  = 0,  cos40  = -1 

Q'n  ■=  76.37  - 19.71 
= 56.66 

This  result  agrees  with  that  shown  in  Table  22. 


i 
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I 


(107) 


The  calculation  of  transformed  modulus  using  the  invariant  formulation 
as  shown  in  fable  21  will  be  identical  to  the  multiple  angle  because 


This  identity  is  true  only  for  the  modulus  transformation.  As  we  shall  see 
later,  the  transformed  compliance  calls  for  negative  signs  in  Eq.  108. 

Then  the  multiple -angle  formulation  is  not  identical  to  the  invariant 
formulation  because  of  this  sign  change. 

The  following  general  remarks  can  be  made  about  the  transformed 
modulus: 

a.  Mirror  image  exists  between  Q'^  and  and  Q'^  ar*d  Q^* 

This  can  be  shown  by  substituting  9+90  into  appropriate  equations 
in  Tables  20  or  21,  and  seen  in  Fig.  35. 


Q'  (9  + 90)  = U + U cos2(0  + 90)  + U cos4(0  + 90) 

1 i X C>  J 

= U.  - U cos29  + U„cos49 
12  3 


(0) 


(109) 


These  two  components  can  be  superimposed  by  a displacement 
of  90  degrees  along  the  0-axis.  We  can  also  show 


Q'26(0)  = -Q,l6(9  + 90)  (HO) 

These  two  components  can  be  superimposed  by  a displacement 
and  a rotation. 

b.  The  frequency  and  amplitude  of  Q'^  and  are  the  same; 
i.  e. , 40  and  U^,  respectively.  (See  Table  20  and  Fig.  36) 

The  two  transformed  components  are  vertically  displaced  by  the 
amount  of  U,  • U,. 
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Transformed  modulus  as  functions  of  Uj . The  relations 
in  Table  20  are  shown  graphically.  Each  transformed 
component  contains  an  invariant  and/or  cyclic  terms. 


1 


c.  The  first  four  transformed  components  in  Table  20  are  governed 
by  cosine  functions.  They  are  even  functions  and  not  sensitive 
to  the  sign  of  ply  orientation.  So  an  error  in  the  sign  will  not 
affect  these  components.  But  the  last  two  components  in  Table  20 
are  the  shear  coupling  components  and  are  governed  by  sine  or 
odd  functions.  A sign  error  will  lead  to  a real  error. 

d.  Because  of  the  symmetry  relations  between  these  transformed 
components,  we  only  need  to  show  three  curves  instead  of  six 

in  Fig.  35;  either  the  three  curves  on  the  left  of  this  figure  or  the 
three  on  the  right.  In  Table  22,  and  can  deleted. 

e.  The  Q’  or  Q*  component  is  made  up  of  three  terms;  one 
constant  or  invariant  term  and  two  cyclic  terms  with  angular 
displacements  2 and  4 times  that  of  the  coordinate  rotat  on;  i.  e. , 
qi  = u + U cos20  + U cos40.  Since  the  cyclic  terms  do  not 
contribute  to  the  area  under  this  curve,  this  area  is  simply  equal 
to  the  value  of  Up  Thus,  this  invariant  represents  the  Young's 
modulus  or  extensional  stiffness  potential  of  this  unidirectional 
composite.  The  cyclic  terms  contribute  to  the  directional  changes 
or  anisotropy  without  any  effect  on  the  total  stiffness  potential. 

The  increase  in  the  stiffness  in  one  direction  must  be  made  up  by 
a decrease  in  some  other  direction  while  the  total  area  under  the 
transformed  modulus  remains  constant  and  invariant.  This  is  like 
an  incompressible  material  that  can  undergo  a shape  change 
without  any  volume  change.  Fig.  36  shows  the  sum  of  Up 

f.  The  Q'  is  the  Poisson's  component  that  has  an  invariant  U . 

Xb  “ 

The  cyclic  term  has  an  amplitude  of  U^.  Similarly,  the 

is  the  shear  component  with  an  invariant  Ug#  and  cyclic 

amplitude  U^»  These  two  curves  can  be  superimposed  over  one 

another  by  a vertical  displacement: 


U5  - U4  ’ ' Q.2  (11,> 


The  last  step  comes  from  Eq.  99.  Any  linear  combinations 
of  invariants  are  invariants.  We,  therefore,  have  a new  but  not 
independent  invariants  in  Eq.  111. 

g.  The  shear  coupling  terms  Q'^  and  Q'  have  no  invariant  associated 
with  the  transformation.  They  are  not  independent  in  the  sense  that 
they  are  derivable  from  Q'^  and  by  differentiation, 
respectively.  From  Table  20  and  21, 


From  Eq.  112(a) 

Q'l6  = w^en  ® = 0 ant*  90,  or 

when  U + 4U  cos  20  = 0,  or 

J 

cos20  = -U2/4U3  (113) 

For  T300/5208  from  Table  22: 


U2/4u3  = 85.73/4  x 19.71  = 1.087  (114) 

There  is,  therefore,  no  solution  for  Eq.  114  or  it  means  that 
this  shear  coupling  does  not  go  to  zero  other  than  at  0 and  90 


rearrange: 

2 

cos  20  + 


2 1 
—co s 20  - — = 0 


solving  for  0: 


cos29  = — 


For  T300/5208, 


cos20  = .485,  -1. 029 

9 = 30.4,  no  solution 


(116) 


At  this  angle,  it  is  a point  of  inflection  in  the  Q'j  1 curve  and  a 

maximum  in  the  Q'  Similarly,  at  59.6  degrees,  Q has  the 

inflection,  and  Q'  , the  maximum.  These  relations  are  shown 
26 

in  Fig.  37. 

Points  of  inflection  of  Q'  can  be  found  from  letting 

io 


= 8U  sin20  + 64U  sin40  = 0 
v3  Z j 


(117) 


8sin20(U2  + 16U3cos20)  = 0 

sin20  = 0 or  0 = 0,  90  for  all  components 
U2 

cos20  = -TZ77-  or  0 = 52.  88  for  T300/5208 
!6U 


(118) 


The  power  functions  formulation  for  the  modulus  transformation 
can  be  used  to  demonstrate  the  dominance  of  the  longitudinal 
properties  of  unidirectional  composites.  Since  the  first  column 
of  Table  19  is  many  times  higher  than  the  other  components;  i.  e. 
181  GPa  versus  10,  3 and  7 for  T300/5208  based  on  the  data  in 
Table  9,  we  can  show  the  contribution  of  the  first  column  or 
that  of  Q'n*  *n  38.  The  dashed  lines  are  those  transformed 

modulus  based  on  the  first  column  of  Table  19;  the  solid  lines 
are  the  complete  solutions,  as  those  in  Fig.  35.  We  can 
see  that  for  our  highly  anisotropic  unidirectional  composite,  this 
approximation  is  fairly  close  to  the  exact. 
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Fig.  37.  Relationships  between  transformed  modulus.  Special 

relationships  are  expressed  in  Eq.  112  to  116.  Note  the 
point  of  inflection  in  Q 1 ^ is  the  point  of  maximum  value 

in  ^16*  At  ° and  de8rees»  aye  zero,  and  the  slopes 
of  Qjj  are  also  zero.  The  points  of  inflection  for  Q'^are 
also  shown. 


By  the  same  token,  we  can  approximate  the  values  of  Us  in  Eq.  99  as 


The  approximate  transformation  equations  can  be  simplified  from 
Table  20  as  follows: 


(3  + 4cos26  + cos40)Q  /8  = m Q 


This  and  the  other  components  of  transformed  modulus  will  be 
the  same  as  the  first  column  of  Table  19.  Thus,  both  power 
functions  and  multiple- angle  functions  for  the  modulus  trans- 
formation reduce  to  the  same  limiting  case  when  only  the 
is  present. 


— — ■■  - ■ __ 


Analogous  to  the  approach  for  the  off-axis  modulus,  we  can  derive  the 
off-axis  compliance  following  the  sequence  of  Fig.  39,  which  is  a repeat  of 


(C)~ 

On-  Axis 
Strain 


* (d) 
Off  -Axis 
Strain 


Fig.  39.  Derivation  of  off-axis  compliance. 

From  (a)  to  (b):  Positive  stress  transformation. 

From  (b)  to  (c):  Stress-strain  relations  in  compliance. 

From  (c)  to  (d);  Negative  strain  transformation. 

We  can  go  directly  from  (a)  to  (d)  by  merging  the  three  steps 
into  one. 


Since  the  derivation  of  the  compliance  transformation  is  analogous  to 
that  of  the  modulus  transformation  from  Eq.  88  to  94,  we  will  write  only  the 
first  line  of  each  equation  for  this  derivation. 


• To  go  from  (a)  to  (b)  in  Fig.  39,  we  know  the  equations  for  stress 
transformation  from  Table  10: 


To  go  from  (b)  to  (c)  in  Fig.  39  we  need  the  on-axis  stress- strain 
relation  in  compliance  from  Table  5. 


6^  = + n2(J^'  + 2mn<^') 

+ Sl2(n2(J^'  + m2(T2'  - 2TOn(T6')  (122) 

• To  go  from  (c)  to  (d)  in  Fig.  19  we  need  the  negative  strain 
transformation  in  Table  14  where  the  sine  functions  now  have 
negative  signs. 

2 2 

€',  = m 6,  + n - mnC, 

1 14  0 

= [m4Sn  + n4S22  + 2m2n2  + xnVs^J  0^' 

+ [m2n2  (Sn  + S22)  + (m4  + n4)S12  - mVs^]  ^ 

+ ^2m  nS^  - 2mn3522  + (mn3  - mn  ) (2S^ ^ ) 1 0^' 

- s’ii<Y  + + s'16<r6'  023) 

Note  that  shear  coupling  terms  appear  in  this  off-axis  unidirectional 
composite,  in  an  analogous  fashion  as  the  off-axis  modulus  Eq.  92  to  94. 

The  off-axis  stress- strain  relation  in  terms  of  compliance  is  presented  in 
a matrix  multiplication  table  in  Table  24  similar  to  the  off-axis  stress- 
strain  relation  in  terms  of  the  modulus  in  Table  18. 


Table  24.  OFF-AXIS  STRESS-STRAIN  RELATION  FOR  UNI- 
DIRECTIONAL COMPOSITES  IN  TERMS  OF  COMPLIANCE 


<r; 

< 

s'„ 

si 

Si 

si 

S'M 

c' 

< 

si 

Si 

s' 

Again,  the  primes  can  all  be  eliminated  without  affecting  the  validity  of  the 
relationship.  As  stated  previously,  die  choice  of  primed  and  unprimed  coor- 
dinates is  arbitrary  and  can  vary  from  unidirectional  to  multidirectional  com- 
ponents. The  symmetry  relations  for  the  transformed  compliance  can  be 
shown  by  including  interaction  terms  (T^CT^  and  <£,0^  in  addition  to  (T^O^  in  the 
stored  energy  expression  in  Eq.  17.  We  can  then  show  that 


S’l6’S'61'  S'26*S,62 


(124) 


The  transformation  equations  for  the  compliance  are  taken  from 
matching  like  terms  in  the  last  two  steps  of  Eq.  123  and  will  be  shown  in 
Table  25.  The  missing  equations  can  be  derived  in  an  analogous  fashion  from 
Eq.  121  to  123.  This  table  is  analogous  to  the  transformed  modulus  in  Table  19. 


Table  25.  TRANSFORMATION  OF  COMPLIANCE  OF  ON- AXIS 

UNIDIRECTIONAL  COMPOSITES  IN  POWER  FUNCTIONS 


•S'/, 

•S',, 

Ss. 

c * 

m4 

n 4 

2m*  n* 

m*ns 

S' 

« 

n 4 

m4 

2m*n* 

m*n * 

S' 

m*n* 

m*n 3 

m4+n 4 

- m*n* 

S' 

°€t 

4m*n* 

4m*n 3 

-8m*n 3 

(m*-n*)* 

I m 

2m3n 

- 2 mn3 

2(mn3~msn  ) 

mn3-msn 

E9 

2mn3 

-2m3n 

2(m3n -mn3) 

m3n  -mn3 

Note  that  the  difference  between  this  table  and  Table  19  for  the  transformation 

of  modulus  can  be  traced  to  the  use  of  engineering  shear  strain.  For  each 

component  with  single  subscript  6,  the  coefficients  on  each  row  shall  be 

multiplied  by  2.  For  components  with  double  subscript  6,  such  SI , the 

o6| 

coefficients  shall  be  multiplied  by  4.  In  the  last  column  of  Table  25,  the 
effect  of  double  subscript  6 is  to  divide  each  coefficient  by  4.  All  the 
differences  between  Tables  25  and  29  can  be  accounted  for  with  these 
corrections. 


The  multiple- angle  formulation  of  the  transformation  of  compliance 
follows  precisely  the  same  pattern  as  that  for  the  transformed  modulus. 

The  multiple- angle  trigonometric  identities  in  Eq.  95  can  be  substituted  into 
the  coefficients  in  Table  25.  By  following  the  process  from  Eq.  96-98  we 
can  derive  the  multiple -angle  representation  of  the  transformed  compliance 
in  a matrix  multiplication  table  as  follows: 


Table  26.  TRANSFORMED  COMPLIANCE  FOR  ON-AXIS 
UNIDIRECTIONAL  COMPOSITES  IN  MULTIPLE 
ANGLE  FUNCTIONS 


/ 

u. 

U3 

u, 

cos20 

cos40 

U, 

-cos 20 

cos40 

u4 

-cos40 

u9 

’-4cos40 

sin20 

2 sin  40 

si  n20 

-2  sin  40 

R1  * 4 (S'll  ' S,22)2  + 4 (S'l6  + S'26)2  ’ U2 


or 


R1  * iU2 


(127) 


R2  - h <S'll  '+  S,22  - 2S’l2  - + 16  ,S,!6  - S'26>2 

■ u2 


or 


R2  = ±U3 


(128) 


From  the  relationship  above,  we  can  derive  the  transformation  equations  in 
terms  of  the  invariants.  As  mentioned  in  the  discussion  of  transformed 
modulus,  the  invariant  formulation  has  one  drawback  which  is  associated 
with  phase  angles.  For  this  reason,  the  multiple -angle  formulation  is  pre- 
ferred because  the  signs  are  self-contained.  For  the  compliance  of  an  on- 
axis  unidirectional  composite,  and  are  negative  if  the  longi- 
tudinal (the  1-axis)  stiffness  is  higher  than  the  transverse  (the  2-axis) 
stiffness.  See  data  in  Table  28.  But  for  completeness  of  all  three  alternative 
formulations  of  the  transformation  of  compliance,  that  of  the  invariant 
functions  is  listed  in  a matrix  multiplication  table  as  follows: 


Table  27.  TRANSFORMED  COMPLIANCE  OF  ON- AXIS 

UNIDIRECTIONAL  COMPOSITES  IN  INVARIANT 
FUNCTIONS 


/ 

R, 

RZ 

S',, 

u, 

-cos20 

-cos 40 

si. 

U, 

cos20 

- cos40 

s'„ 

U4 

cos40 

Sgg 

<JS 

4cos40 

S’,e 

~ si  n20 

-2  sin  40 

Sie 

- si n20 

2 sin  40 

9 3 

“X 

*“  ?;rv,v  ■;/  " • 

£ 

Note  that  the  signs  of  the  trigonometric  functions  are  changed  from  those  in 
Table  26  because  U_  and  U have  negative  values.  So  negative  signs  must 

£ j 

be  used  in  Eq.  127  and  128,  or 


R1  = -U2 


(129) 


r2  = -u3 


This  choice  of  negative  valueB  is  different  from  the  invariant  functions  of 
the  modulus  transformation  in  Table  21  where  the  positive  values  of  and 
were  picked.  This  was  bo  because  ^ is  greater  than  for  most 
unidirectional  composites. 


Finally,  compliance  components  also  provide  the  link  to  engineering 
constants.  Modulus  can  provide  this  link  only  when  the  material  is  ortiio- 
tropic,  or  in  its  symmetry  axes.  But  compliance  does  not  have  this  re- 
striction. This  can  be  easily  demonstrated  by  applying  a uniaxial  stress  in 
the  1' -direction  of  an  off-axis  coupon.  Then  from  the  stress-strain  rela- 
tion of  an  off-axis  material  in  Table  24,  we  have 

= S'  0"  ' , 6 ' =S'  <T'  6'  -S'  ,<T'  (130) 


1 ~ in  ’ “2  ~ 21  1 ' ~ 6 61  1 

We  can  immediately  establish  the  off-axis  Young's  modulus  and  Poisson's 
ratio  as 

(131) 


E'l  = 1/ S'n 


V * -s'2i/s'n 


(132) 


The  sign  of  the  induced  shear  strain  in  Eq.  130  depends  on  the  sign  of  the 
shear  coupling  S'^,  assuming  the  applied  stress  is  tsnsile. 

Fig.  40  shows  the  two  possible  shear  strains  induced  by  an  off-axis  uniaxial 


94 


*2*.^  /- 


■ ' ... 


1 

1 


tensile  test 


f 


I 


tr 


*/' 


(a) 


FiS*  40.  Off-avi  o 

axis  uniaxial  -i 

addiuon  to  n°rmai  8;ai;r;  .™8  te8t — m^ce,  in 

:°;hpling  — - i8  posz;^  8train  if  •— 

e c°efficient  is  negat.ve  < >>  °r  dative  ahear, 

Positive  ply  orientation  in  (a)  ind  C ' ^ m°8t  comPosit( 

beca-u2---flne:a^ 

4*  EXAMPLES  of  OFF-AXIS  C 

We  ahow  in  this  OMpLlANCE 

**•  orthotropic  com  8eCti°n  **  *«•*>«**  Com  .. 
compoaite.  When  AP°nents  of  complia„Ce  are  ii8t  ! .lSnCe  ^ T3°0/52< 

in  Table  25,  **'  8Uba«tuted  into  thecal  f^  ‘ ^ ^ 

Alternatively  w " ^ ^nee. 

we  uae  th.  r „ /'  e can  arrive  at  the  8am  . 

compliance  t ValU68  °*  Ui  in  the  multipi  Fan8f°rmed  COI*Ponenta  if 

given  in  i tl"0"™'0”  ..  Table  *,  „ ' of  th. 

tnaneformed  „ T ^ ^ ^ 

m®d  compliance  can  then  be  '“-"PoaHe..  Tbe 

The  numerical  reaulf  computed  from  the  , 

8 arC  li8ted  in  Table  29  and  ^ ati°n8  *“  Table  26‘ 

"d  eU""*  in  Fig.  41. 


95 


Table  28.  TYPICAL  VALUES  OF  LINEAR  COMBINATIONS  OF 
COMPLIANCE  FOR  ON-AXIS  UNIDIRECTIONAL 
COMPOSITES  (TPa)-l 


T300/5208 

B{4)/5505 

AS/3501 

Scotchply 

/1002 

Kevlar  49 
/Epoxy 


U1 

U2 

U3 

U4 

U5 

55.  53 

-45.78 

-4.22 

-5.77 

122.6 

43.42 

-24. 55 

-13. 94 

-15.06 

117.0 

61.62 

-52.18 

-2.20 

-4.38 

132.0 

83.5 

-47.5 

-10.2 

-16.9 

200.8 

126.4 

-84.33 

-28.86 

-33. 33 

319.4 

Table  29.  TRANSFORMED  COMPLIANCE  FOR  T300/5208 
UNIDIRECTIONAL  COMPOSITES  (TPa)'1 


Sil 

S’ 

22 

Si2 

S66 

S1 6 

S26 

5.52 

97.09 

-1.55 

139.4 

0 

0 

13.  77 

93.06 

-3.66 

131.0 

-30.20 

-15.58 

34.  75 

80.53 

-7.88 

114.1 

-46. 96 

-32. 34 

59.75 

59.75 

-9.99 

105.7 

-45.78 

-45.78 

80.53 

34.75 

00 

00 

• 

r- 

1 

114.1 

-32.34 

-46.96 

93.  06 

13.77 

-3.66 

131.0 

-15.58 

-30. 20 

97.09 

5.52 

-1.  54 

139.4 

0 

0 

The  general  remarks  on  the  transformed  compliance  are  very  similar 
to  those  on  die  modulus.  We  will  simply  repeat  the  relevant  features  without 
further  detailed  discussions. 

a.  Mirror  image  exist  between  S'.  and  S'  , and  S'  , and  S'  , . This 

11  CC  lo  £0 

can  be  seen  from  Table  26  and  Fig.  41. 

b.  The  amplitude  of  S'j2  is  now  one  quarter  that  of  S'^.  The  fre- 
quency of  49  remains  the  same  for  both  transformed  components. 

See  Table  26  and  Fig.  41. 

c.  Again,  only  the  shear  coupling  components  are  affected  by  the  sign 
of  the  angle  of  rotation. 

d.  Because  of  the  symmetry,  only  three  transformed  components 
Q'u#  Q'12  (or  Q.  ).  and  Q'^  need  to  be  drawn. 

e.  In  Fig.  42  the  dashed  lines  show  the  approximation  of  the  trans- 

1 

formed  components  using  only  the  first  one  or  two  columns  of 
the  complete  transformation  equations  in  Table  26.  Because  of 
the  particular  values  for  T300/5208,  the  approximations  give 
excellent  results.  This  is  analogous  to  the  transformed  modulus 
shown  in  Fig.  36. 

f.  In  Fig.  43  , the  dashed  lines  show  the  contributions  of  the  two  pre- 
dominant components  S 22  an<*  S^,  assuming  and  S12  are  zero. 

These  compliance  curves  are  expected  to  be  controlled  more  by  the 
matrix  than  by  the  fiber.  That  is  why  the  two-term  approximations 


in  Fig.  43  agree  well  with  the  exact  compliance.  The  opposite  was 

true  in  the  modulus  where  the  dominant  role  was  played  by  Qjj  or 

the  longitudinal  stiffness  E . This  was  shown  in  Fig.  38. 

L 


U2Sin20  Sm20 

-60  I* 

Fig.  42.  Comparison  of  exact  and  approximate  transformed  compliance 
in  terms  of  the  multiple -angle  functions  for  T300/5208.  The 
dashed  lines  are  approximations  without  the  last  column  in 
Table  26  or  U,»0. 
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Fig.  43.  Contributions  of  S22  and  S&6  to  the  transformed  compliance 

for  T 300/5208.  The  solid  lines  are  the  exact;  the  dashed  lines 
the  contribution  of  S22  and  terms.  The  Sjj  and  S^2  are 
taken  to  be  zero. 
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Relationships  exist  between  the  derivatives  of  the  S'  and  S'  , 

X X bfa 

and  the  shear  coupling  components  as  follows: 

&S' 


11 


= -2  S' 


16 


*S'22 

SiT*  - 2S' 


26 


(133) 


(134) 


h. 


From  these  relations,  we  can  easily  obtain  the  angles  where 
maxima,  minima,  points  of  inflection  exist  in  the  transformed 
components  of  compliance,  similar  to  Fig.  37  for  the  modulus. 
Engineering  constants  will  also  vary  with  ply  orientation.  These 
constants,  however,  are  not  governed  by  any  transformation 
equations  like  modulus  and  compliance,  although  they  can  be 
related  to  transformed  components  of  compliance  as  indicated 
in  Eq.  131  and  132.  We  can  similarly  define  the  transformed 
engineering  shear  modulus  as: 


G'l2  * 1/S'66 


(135) 


These  transformed  engineering  constants  can  be  computed  from 

the  transformed  compliance  data  in  Table  29,  and  are  plotted  in 

Fig.  44  as  functions  of  ply  orientations.  Also  shown  in  this  figure, 

as  dashed  lines,  are  the  transformed  components  of  Q'  and  Q'  , 

1 1 66 

Note  that  these  modulus  components  have  identical  or  nearly 
identical  values  when  ply  orientation  is  0 or  90.  For  example; 
when  § = 0 


From  Table  9#  Qjj  = 181.8  GPa 

• 

Q22  = 10.34 

If 

(136) 

[ 

Q66  = 7*17 

ri 

From  Table  7,  E_  * 181 

Li 

it 

E * 10.3 

T 

n 

(137) 

- *4 

GLT  = 7-17 

it 

101 

* » * . 

Poisson's  Ratio 


Young's  Modulus,  GPa 


Shear  Modulus,  GPa 


Fig.  44.  Transformed  engineering  constants  for  T300/5208. 

T ransformed  components  of  modulus  taken  from  Table  22  and 
Fig.  35  are  also  shown.  Note  the  significant  difference  of  the 
transformed  quantities  between  the  modulus  components  and  the 
engineering  constants.  In  order  to  avoid  making  big  mistakes, 
we  recommend  that  engineering  constants  be  limited  to  those 
applied  to  the  symmetry  axis  only,  like  those  in  Table  7. 
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In  off-axis  orientations,  the  transformed  modulus  and  the 
transformed  engineering  constants  are  no  longer  close.  The 
former  can  be  severed,  folds  higher  than  the  latter.  In 
laminated  and  built-up  structures  the  most  direct  material 
properties  for  stiffness  are  the  modulus  components,  not  the 
engineering  constants.  The  basic  issue  lies  in  the  fact  that 
engineering  constants  are  essentially  one-dimensional  properties. 

It  is  with  great  difficulty  and  with  much  unnecessary  complexity 
if  we  insist  on  using  engineering  constants  to  describe  the 
elastic  behavior  of  laminated  and  built-up  structures.  Unfortunately, 
many  design  and  experimental  procedures  are  built  on  one- 
dimensional concepts  where  engineering  constants  abound.  This 
is  a carryover  by  those  who  have  gotten  accustomed  to  the  use 
of  ordinary  materials. 
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INVERSE  RELATIONSHIP  BETWEEN  MODULUS  AND  COMPLIANCE 
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The  off-axis  stress-strain  relations  as  listed  in  Tables  18  and  24, 
are  based  on  modulus  and  compliance,  respectively,  are  repeated  here 
as  Tables  30  and  31  esccept  all  primes  have  been  removed. 


Table  30.  OFF-AXIS  STRESS-STRAIN  RELATION  FOR 
UNIDIRECTIONAL  COMPOSITES  IN  TERMS 
OF  MODULUS 


0„ 

o,* 

0« 

0*/ 

Ogg 

0*e 

0*, 

06* 

Oss 

Table  31.  OFF-AXIS  STRESS-STRAIN  RELATION  FOR 
UNIDIRECTIONAL  COMPOSITES  IN  TERMS 
OF  COMPLIANCE 


s„ 

S,* 

S,s 

s*, 

Sgg 

S*s 

Ss, 

Ss* 

S$s 

The  difference  between  these  stress- strain  relations  is  that  the  role  of  stress 
and  strain  are  the  inverse  of  each  other.  In  Table  30,  the  strain  is  the  inde- 
pendent variable;  in  Table  31,  the  stress  is  the  independent.  We  inverted  the 
on-axis  stress-strain  relations  in  Section  1 when  we  went  from  Eq.  8 to  11  by 
simply  solving  the  simultaneous  equations.  We  need  only  to  repeat  the  same 
process  for  the  off-axis  case,  where  shear  coupling  terms  are  no  longer  aero. 


We  can  proceed  with  the  inversion  or  solution  of  these  equations  by  the 
method  of  determinant  as  follows: 

We  will  assume  that  we  are  given  the  equations  in  Table  30.  We  will 
first  obtain  the  determinant  of  the  modulus  components: 


Determinant  of  Modulus  = det  Q. . - A 

Q11Q22Q66  + 2Q12Q26Q61  " Q22Q16  " Q66Q12  “ Qll°62 


(138) 


^1  = <Q22Q66  - °26>  ' A 


22 


(Q11Q66  " Q16)  / A 


12 


66 


= (Q16Q26  “ Q12Q66)  / A 


^11^22  ” ^12^  ^ 


(139) 


’16 


(Q12Q26  " Q22Q16)  / A 


26 


<Q12Q16  “ Q11Q26)  ^ 


We  have  gotten  the  components  of  compliance  from  those  of  modulus.  If  we 

are  given  the  compliance  and  want  to  know  the  modulus,  we  simply  interchange 

the  Q..  for  S. . in  Eq.  138  and  139. 
ij  iJ 
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Thus  there  are  two  ways  that  we  can  compute  the  off-axis  modulus  or 
compliance.  This  is  diagrammed  in  Fig.  45.  Our  starting  point  consists 
of  a set  of  engineering  constants  as  those  listed  in  Table  7.  There  are  2 ways 
of  getting  the  transformed  compliance  and  modulus: 

a.  From  engineering  constants,  compute  on-axis  compliance  and 
modulus,  shown  in  Tables  8 and  9,  respectively.  Compute 
transformed  modulus  using  its  transformation  equations  in 
Table  20  and  compute  transformed  compliance  using  those  in 
Table  26. 

b.  Alternatively,  we  can  go  directly  from  the  transformed  modulus 
to  the  transformed  compliance  by  inversion  in  Eq.  138  and  139 
or  the  transformed  compliance  to  the  transformed  modulus  also 
by  inversion. 


Transformation  t |^7j 

f 

Inversion 


Transformation  \~T7 

•*  Sj) 


Fig.  45. Relation  of  off-axis  compliance  and  modulus.  Can  be 
obtained  in  two  ways:  (a)  separate  transformations  of 
modulus  and  compliance;  (b)  transformation  of  modulus 
then  followed  by  inversion  to  obtain  transformed  compliance; 
or  transformation  of  compliance  then  followed  by  inversion 
to  obtain  transformed  modulus. 

From  the  operational  standpoint,  the  inversion  following  transformation,  the 
2nd  method,  is  a little  easier  to  perform  than  the  1st  method,  the  two  trans- 
formations. But  the  difference  is  small.  We  should  be  well  versed  in  both 
methods  and  choose  the  one  that  is  best  for  a given  situation. 
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SECTION  IV 

IN-PLANE  STIFFNESS  OF  SYMMETRIC  LAMINATES 

SCOPE 

The  stiffness  of  multidirectional  laminates  consisting  of  plies  and  ply- 
assemblies  with  arbitrary  ply  orientations  will  be  described.  The  com- 
posite laminates  are  limited  to  those  stacking  sequences  which  have  mid  - 
plane  symmetry;  i.  e. , ply  orientations  in  the  lower  half  of  the  laminates 
are  exactly  the  same  as  those  in  the  upper  half.  Such  laminates  will 
behave  like  a homogeneous  anisotropic  plate.  We  will  see  that  the  effective 
modulus  of  the  composite  laminate  is  simply  the  arithmetic  average  of  the 

modulus  of  the  constituent  plies.  Simple  formulas  and  charts  can  be  made 

/ 

so  laminate  modulus,  in  terms  of  the  constituent  plies,  can  be  quickly 
established.  The  key  relation  for  the  in-plane  modulus  of  any  laminate  is: 


PRINCIPAL  NOMENCLATURE  AND  DEFINITIONS 


In-plane  modulus  of  multidirectional  laminates,  in  Nm  * f 
or  Pam. 

In-plane  compliance  of  multidirectional  lamina it  is  the 

-1  -1 

inverse  of  A..,  in  N m,  or  (Pam)  . V 

An  in-plane  engineering  constant:  the  Young's  modulus 
along  the  1-axis. 

An  in-plane  engineering  constant:  the  Young's  modulus 
along  the  2 -axis. 

An  in-plane  engineering  constant:  the  longitudinal  sljear 
modulus. 

Total  thickness  of  a laminate. 


n.  h = Thickness  of  i-th  ply  assembly  (one  or  more  plies 
i o 

having  the  same  ply  orientation  0 );  i = l,2,...m. 
h/n  = Unit  ply  thickness. 


n^h^  = Thickness  of  ply  assembly  having  n plies  at  0 ply 

orientation, 
cos  0,  or 

Total  number  of  ply  assemblies  in  a laminate. 

Stress  resultant,  in  N i = 1,2,6. 


sin  0 or 

Total  number  of  plies  in  a laminate 


h /h 
i o 

Number  of  plies  in  i-th  ply  assembly;  i = 1,  2, . . . m. 


Number  of  plies  having  0 ply  orientation. 

Modulus  of  ply  assembly  with  9 ply  orientation. 

Modulus  of  i-th  ply  assembly;  i = 1,  2,  ...m. 

Linear  combinations  of  modulus  for  the  multiple -angle 
transformation;  i = 1,2,3,  4,5. 

Integrals  of  trignometric  functions  for  the  evaluation 
of  the  in-plane  modulus;  i = 1,  2,  3,4. 

nQ/n 

Volume  fraction  of  ply  assembly  with  0 orientation. 

Nj/h  = Average  stress  across  thickness  of  a laminate; 
i - 1,  2,  6. 

In-plane  strain  components;  i = 1,2,6. 
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LAMINATE  CODE 

A multidirectional  composite  laminate  is  defined  by  the  following 
code  to  designate  the  stacking  sequence  of  ply  assemblies: 

p)3/902/45/-453J  5 (140) 

This  code  is  represented  diagrammatically  in  Fig.  46  and  47,  and  contains  the 
following  features: 

• Starting  from  the  bottom  of  the  plate,  at  z = h/2,  the  first 
ply  assembly  has  three  plies  of  0-degree  orientation; 
followed  by  the  next  assembly  with  two  90-degree  plies; 
followed  by  one  45-degree  ply;  and  finally  the  last  assembly 
with  three  -45-degree  plies. 

• The  subscript  S denotes  that  the  laminate  is  symmetric  with 
respect  to  the  midplane  or  the  z = 0 plane.  The  upper 
half  of  the  laminate  is  the  same  as  the  lower  half  except  the 
stacking  sequence  is  reversed  in  order  to  maintain  the  mid- 
plane symmetry. 

• A subscript  T is  used  to  designate  the  total  laminate,  with- 


I 


2.  IN-PLANE  STRESS-STRAIN  RELATION  FOR  LAMINATES 
In  the  derivation  of  the  stress -strain  relation  of  a multi- 
directional laminate,  we  need  to  make  the  following  simplifying  assump- 


tions: 


The  laminate  is  symmetric;  i.  e. , 
0(«)  = 0(-z) 


(143) 


Q..(z)  = Q..(-z) 
ij  lJ 


(144) 


Thus,  both  the  ply  orientation  and  the  ply  material  modulus 
are  symmetric  with  respect  to  the  midplane  of  the  laminate. 


The  strain  remains  constant  across  the  laminate  thickness. 
We  will  use  superscript  zero  to  signify  the  assumed  constant 
in-plane  strain  components  as  follows: 

€j(z)  = e° 


= 6 2 


(145) 


€6(z)  = 66 


Note  that  there  is  no  z-dependency  or  z-variation  across  the 
laminate  thickness.  This  assumption  is  reasonable  when  the 
thickness  of  the  laminate  is  small  in  comparison  with  the 
length  and  width  and  the  laminate  is  symmetric. 


■ I 


s » ■■  . • 


Since  the  stress  distribution  across  the  multidirectional  laminate  is 

not  constant  because  the  modulus  varies  from  ply  to  ply,  it  is  much  easier 
to  define  an  average  stress  than  an  actual  stress  across  the  laminate.  This 

average  stress  can  be  used  to  define  the  stress-strain  relation  of  the 
laminate.  The  stress,  in  this  case,  will  be  the  average  stress,  and  the 
strain,  the  in-plane  strain  in  Eq.  145.  We  can  then  calculate  the  ply 
stress,  or  stress  at  any  ply  within  the  laminate  from  the  in-plane  strain. 

We  will  show  this  later  after  the  stress -strain  relation  for  the  laminate  is 
established.  The  average  stress  is  defined  as  follows: 


(Tdz 


In  Fig.  48,  we  show  the  relationship  between  the  actual  stress  from  ply 
to  ply  and  the  average  stress  across  the  laminate  by  the  averaging  pro- 
cess of  Eq.  146. 


Fig.  48.  Definition  of  average  stress.  Comparison  between  the 
corresponding  components  of  the  actual  ply  stress  and 
the  average  laminate  stress  is  shown. 


Up  to  this  point,  we  have  used  primed  coordinates  (or  the  applied 
stress  or  strain,  and  unprimed  coordinates  for  the  material  symmetry 
axes.  This  was  shown  in  Fig.  34(a).  When  we  go  beyond  the  properties 
of  unidirectional  composites  into  laminated  composites  and  structures, 
it  is  better  to  change  our  notation  to  the  unprimed  as  the  coordinates  for 
our  laminate  or  structures.  A multidirectional  laminate  and  its  reference 
coordinate  1-2  are  shown  in  Fig.  34(b)  and  repeated  here  in  Fig.  49. 

Since  there  are  many  ply  orientations  in  this  laminate,  we  will  have  to 
identify  the  specific  orientation  for  stress  or  strain  transformation.  The 
use  of  primed  coordinates  will  not  be  adequate  if  more  than  one  set  of 
material  symmetry  axes  exist. 


/ 


Fig.  49.  Reference  coordinate  system  1-2  for  typical  multidirectional 
laminates.  Differing  from  Fig.  34(a)  the  primed  coordinates 
are  not  used  here. 

Substituting  the  stress -strain  relation  for  any  ply  orientation  into 
Eq.  146,  we  have 

■ h*/[Que.  * Q12e2  + °16e6]  dZ  (147) 

Substituting  the  assumed  constant  strain  in  Eq.  145,  we  have 

5 ’hfori  + Q12S2  + Q16e6]d*  <U81 

Since  the  in -plane  strain  components  are  independent  of  z,  we  can  factor 
them  out  of  the  integral  signs.  Only  the  modulus  is  left  Inside  the 
integral  because  the  modulus  varies  from  ply  to  ply  depending  on  each 
ply  orientation. 


1 


i 


■ t/v*  6°  ♦ * A**  e6 

(149) 

= h [AU®1  + A|2€2  * A1666] 

(150) 

Similarly, 

*2  = ~ [*12®  i + A22S2  + A26®6] 

(151) 

\ * T [A16S°  + A26®2  * A66®6] 

(152) 

where 

Aji  = /Qud..  A22  = /Q22d.  , 

A12=  /Q  12dz’  A66  = /Q 66d2 

A16=  /°16d2-  A26  * /Q 26dl 

► 

(153) 

where  A„  is  the  equivalent  modulus  for  a multidirectional  laminate.  This 
modulus  is  simply  the  average  of  the  modulus  of  the  constituent  plies. 

There  is  a difference  of  a length  in  the  physical  dimension  of  modulus 

-2  -1  ' 

Q^.  in  Pa  or  Nm  ; and  that  of  A„  in  Nm  . v 

We  can  further  define  stress  resultants  as 


Ni  * hal 


N,  = h(T 
2 2 


(154) 


N, 


hir6 


-1 


Note  the  unit  of  stress  resultant  is  Nm  , or  force  per  unit  width  of  a 
laminate  with  thickness  h.  The  in-plane  stress-strain  relation  for 
laminates  is  actually  the  stress  resultant  versus  in-plane  strain  relation. 
The  latter  is  derived  by  combining  Eq.  154  with  150  et  al.  We  have: 
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1 


A1161  + ^2 


+ A16€6 


A12G1  + A22€2  + A26e6° 


(155) 


N = A ,€°  + A,,€°  + A,, 6° 

6 16  1 26  2 66  6 

This  set  of  simultaneous  equations  can  be  inverted  to  yield  the  in-plane 
strain  in  terms  of  the  stress  resultant.  This  process  is  exactly  the  same 
as  that  described  in  Section  3.  5,  and  the  resulting  determinant  in  Eq.  138, 
and  the  inverted  modulus  components  or  the  compliance  components  in 
Eq.  139.  In  other  words,  Eq.  155  is  based  on  modulus  A.,  of  the  laminate, 
and  v/e  wish  to  find  the  corresponding  compliance  a„  by  inversion  such  that 


- anNl  + a12N2  + a!6N6 


= a,„N  + a__N_  + a_ . N 
12  1 22  2 26  6 


(156) 


‘ a16Nl  + a26N2  + a66N6 

where  a.,  take  the  place  of  S..  in  Table  29,  and  A.,  take  the  place  of  Q 
13  lJ  lJ 

in  Table  22.  We  will  show  both  stress -strain  relations  in  Eq.  155  and  156 

dn  matrix  multiplication  tables  as  follows: 

Table  32.  IN- PLANE  STRESS-STRAIN  RELATION  OF  SYMMETRIC 
LAMINATES  IN  TERMS  OF  MODULUS 


< 

«? 

N, 

A„ 

A/, 

A/, 

", 

A,i 

A„ 

A„ 

", 

A «/ 

Af? 

A$s 

Table  33.  IN- PLANE  STRESS-STRAIN  RELATION  OF  SYMMETRIC 
LAMINATES  IN  TERMS  OF  COMPLIANCE 


N, 

N* 

Ns 

< 

an 

°,6 

< 

o 

a» 

°*s 

°6, 

as* 

These  stress -strain  relations  are  valid  for  the  in-plane  deformation  of  sym- 
metric laminates.  If  stress  resultants  are  given,  we  can  find  the  induced 
in-plane  strain  immediately  from  Table  33.  Then  the  on-axis  ply  strain 
can  be  obtained  by  strain  transformation  from  the  initial  1-2  axes  to  the 
orientation  of  a specific  ply  or  ply  assembly.  The  ply  stress  is  nothing 
more  than  the  ply  strain  multiplied  by  the  on-axis  modulus.  The  complete 
process  going  from  stress  resultants  to  on-axis  ply  strain  and  ply  stress 
is  illustrated  in  Fig.  50. 


(a) — 
Stress 
Resultants 


• (b)  — 

In -Plane 
Strain 


+ (c)  - 

On  - Axis 
Ply 
Strain 


- (d) 

On -Axis 
Ply 
Stress 


Fig.  50.  On-Axis  ply  strain  and  stress  calculations.  From  given  stress 

resultants  applied  to  a multidirectional  laminate;  we  can  go 

From  (a)  to  (b):  use  in-plane  stress-strain  relation  in 

laminate  compliance  from  Table  33. 

From  (b)  to  (c):  use  positive  strain  transformation  in  Tables 

14  et  al.  The  primed  strain  in  Table  14 

o 

shall  be  replaced  by  in -plane  strain  6 ^ 

From  (c)  to  (d):  use  the  on -axis  stress- strain  relation  of 
unidirectional  composites  in  terms  of 
modulus,  in  Table  6. 


os  r*rv,vv 


From  the  compliance  a_  in  Table  33,  we  can  calculate  the  effective 
engineering  constants,  following  the  process  used  in  the  off-axis  unidirec- 
tional composites  in  Eq.  130,  et  al.  We  will  have: 


In-plane  longitudinal  modulus  = E°  = ■■  * 

1 auh 


In-plane  transverse  modulus  = E°  = — * , ■ 

2 


In-plane  Poisson's  ratio  = v 


12 


a 


12_ 

11 


In -plane  shear  modulus  = G 


12  ' 


(157) 


Again  we  want  to  emphasize  that  engineering  constants  are  the  constants 
associated  with  simple  tests  such  as  uniaxial  tensile  and  compressive  tests 
and  simple  shear  tests.  They  are  the  results  of  1-dimensional  tests  and 
represent  1-dimensional  characteristics  of  laminates.  But  composites  are 
rarely  used  in  1-dimensional  configuration.  The  2-dimensional  properties 
of  composites  are  much  different  from  the  1-dimensional  properties  of 
ordinary  materials.  Engineering  constants  are  not  easy  to  work  with  and 
can  lead  to  big  errors.  We  discussed  this  issue  on  off-axis  unidirectional 
composites  in  Fig.  44,  and  will  do  it  again  on  multidirectional  laminates 
later  in  this  section. 

3.  EVALUATION  OF  IN -PLANE  MODULUS 

We  have  only  mentioned  that  the  in-plane  modulus  of  a multi- 
directional laminate  is  the  arithmetic  average  of  the  off-axis  modulus  of 
the  individual  plies  or  ply  assemblies.  The  averaging  process  is  shown 
as  integrals  in  Eq.  153.  We  will  now  describe  the  steps  needed  to  per- 
form the  integrations  so  that  the  contribution  of  the  ply  modulus  to  the 
laminate  modulus  can  be  defined. 


cos49  dz 


We  can  repeat  the  process  for  other  in-plane  components 


U.cos49 


cos46  dz 


A16 

= /Vs  = 

J^Uj/2  sin26  + U3»in4el  dz 

- iU2V3 

+ U3V4 

(166) 

A26 

= 4 u,vo 

2 2 3 

- °3V4 

(167) 

where 

/•  h/2 

v3 

= I sin29  dz 

(168) 

J -h/2 

/.  h/2 

v4 

= / sin48  dz 

(169) 

./  -h/2 

We  have  seen  that  the  evaluation  of  the  in-plane  modulus  A.,  is  reduced  to 

U 

the  evaluation  of  four  integrals,  defined  by  to  V , These  relations  can 
be  summarized  in  a matrix  multiplication  table  as  follows: 


Table  34.  FORMULAS  FOR  IN-PLANE  MODULUS  OF  LAMINATES 


/•h/2 

r 

[l,2,3,4]7  _h/2 

cos26,  cos  46,  sin26«  8ia46| 

(170) 
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This  is  a condensed  definition  of  the  V.  where  the  numerals  in  the  bracket 
on  the  left-hand  side  applies  to  the  corresponding  term  on  the  right-hand 
side  of  Eq.  170.  The  value  of  V.  are  dependent  on  the  variation  of  ply 
orientations  in  the  multidirectional  laminate.  It  is  implicitly  assumed 


that  the  laminate  consists  of  plies  of  the  same  unidirectional  composite. 


Because  sine  and  cosine  functions  are  bounded  between  -1  and  + 1 , the  V's 
are  bounded  by  the  same  limits,  as  we  soon  shall  see.  The  similarity 
between  Tables  34  and  20  is  the  result  of  the  same  transformation  equations. 
In  the  limit  when  the  laminate  has  only  a ply  orientation,  we  recover  Table 
20  from  34  because  the  integrants  in  Eq.  170  are  constant.  The  V.  be- 
come simply  the  trigonometric  functions  times  the  laminate  thickness: 

Vj  = hcos29 


- hcos40 

V = hsin29 
3 


hsin40 


Before  we  perform  explicit  integrations  for  the  evaluation  of  V. 
from  Eq.  170,  we  must  define  the  geometric  quantities  of  a symmetric 
laminate  in  Fig.  51. 


Fig.  51.  Definitions  of  terms  in  a symmetric  laminate. 

The  index  for  ply  assembly  goes  from  0 to  m/2  when 
m is  the  total  number  of  ply  assemblies.  The  i-th 
ply  assembly  is  located  between  z^^and  z^,  and  has  a 
thickness  h:. 


I 


$ 


Let  us  now  evaluate  the  first  V.  in  Eq.  170.  We  first  normalize  it 


with  respect  to  the  laminate  thickness: 


v,  *T 


1 I / 

r ' V. 


h/2 


cos26  dz 


(172) 


h/2 


If  the  laminate  is  symmetric,  we  only  need  to  evaluate  on  half  of  the  thick- 
ness, say,  from  z = 0 to  z = h/2.  The  new  limits  of  integration  call  for  a 
new  interpretation  of  the  laminate  code  as  defined  in  Eq.  140.  The  starting 
point  of  the  ascending  ply  sequence  has  been  reversed  from  the  z = 0 to 
z = h/2.  Only  the  upper  half  of  a symmetric  laminate  needs  to  be  evaluated. 
Thus, 


1 -hio 


h/2 


cos28  dz 


(173) 


Since  each  ply  assembly  is  assumed  to  have  the  same  ply  orientation  and 
material,  this  integration  can  now  be  replaced  by  summation  as  we  move 
from  ply  assembly  to  ply  assembly. 


* 2 

\ ’ h L co.ae,  [z.  - Vl] 


(174) 


_2_ 

h 


cos20.  h 
i i 


(175) 


where  h.  = thickness  of  i-th  ply  assembly;  where  i begins  from  the  midplanex 
Let  v.  = volume  fraction  of  plies  with  0.  orientation 

(176) 


2h./h 
i o 


J 


If  each  index  i in  Eq.  175  represents  a unique  ply  orientation,  we  can  now 
substitute  Eq.  176  into  175. 


(177) 


* m/2 


i 


cos26  v. 
i i 


where 


= VjCos29^  + v^cos202  + v^cos2fi^  + . . . 


vl  + V2  + v3  + 


= l 


(178) 

(179) 


Thus,  is  simply  the  rule  of  mixtures  equation,  or  the  weighted  average 
of  the  cos26  functions.  Since  cosine  functions  can  never  be  greater  than 
unity  (or  less  than  minus  unity),  each  term  in  Eq.  178  is  always  equal  to 
or  less  than  the  corresponding  term  in  179.  We  can  therefore  conclude 


that  is  bounded  as  follows: 


>!< 

-1  < V < 1 


(180) 


By  applying  the  identical  process  to  the  remaining  V 
get  the  following: 

* 


* 

V 

3 


VjCos49^  + 


v sin20  + 
1 1 


V CO840 
2 2 


v_sin20. 
2 2 


+ v^cos40 

+ v sin20 
3 3 


in  Eq.  170,  we  can 

3 + . . . (181) 

+ . . . (182) 


V*  = v^sin40^  + v^sin40^  + v^sin40^  + . . . (183)  . 

With  these  simple  equations,  we  can  easily  compute  the  in-plane  modulus 
of  multidirectional  laminates  with  any  ply  orientations.  The  information 
needed  is:  each  ply  orientation  and  the  volume  fraction  of  that  ply  orienta- 
tion; i.  e. , 8^  and  v^,  respectively.  Then  from  Eq.  178,  181-183,  we  can 
calculate  the  V\*.  From  Table  34  we  can  compute  the  A.^  for  any  multi- 
directional laminates. 
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When  normalized  V\  or  V.*  are  used.  Table  35  can  be  rewritten  in  a 
matrix  multiplication  table  as  follows: 


Table  35.  FORMULAS  FOR  NORMALIZED  IN-PLANF 
MODULUS 


/ 

4 

“s 

> 

U, 

V2 

Agg/h 

U, 

V* 

v2 

A,2/h 

A^/h 

-V? 

A,e/h 

V*/2 

•4* 

A 26 

V*/2 

~v* 

The  units  of  normalized  in-plane  modulus  will  be  the  same  as  the  modulus 

4 

of  unidirectional  composites.  Direct  comparison  between  unidirectional 
and  multidirectional  composites  is  possible  with  the  normalized  in-plane 
modulus.  If  the  laminate  is  unidirectional,  the  V.  are  those  listed  in 
Eq.  171.  Using  the  definition  of  the  normalized  V\  in  Eq.  172,  the  values 
in  Eq.  171  can  be  substituted  into  Table  35.  We  will  recover  the  trans- 
formation relation  of  unidirectional  composites  in  Table  20. 
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4.  CROSS -PLY  LAMINATES 

We  will  now  examine  some  commonly  encountered  symmetric  lamin- 
ates and  determine  the  values  of  their  in-plane  modulus  and  compliance. 

First,  we  will  study  cross -ply  composites.  The  ply  orientations 
are  limited  to  0 and  90  degrees.  In  Table  36,  all  the  values  of 
the  trigonometric  functions  which  will  be  needed  are  listed. 

Table  36.  VALUES  OF  TRIGONOMETRIC  FUNCTIONS  FOR 
CROSS- PLY  LAMINATES 


ei 

cos20. 

cos40. 

l 

sin28. 

sin  46. 

l 

0 

1 

i 

0 

0 

90 

-1 

i 

0 

0 

Substituting  these  trigonometric  functions  into  Eq. 

* 


178,  et  al,  we  have 


(184) 


where  v = the  volume  fraction  of  the  0 ply  orientation.  By  taking  these 
8 

values  and  putting  them  into  Table  35,  we  will  have  the  normalized  in- 
plane modulus  for  cross -ply  composites  as  functions  of  volume  ratios. 
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This  is  done  in  Table  37  where  matrix  multiplication  is  Implied. 


Table  37.  FORMULAS  FOR  IN-PLANE  MODULUS  FOR 
CROSS- PLY  COMPOSITES 


/ 

4 

A„/h 

u, 

v0  ~ v90 

/ 

A„/h 

U, 

v9o  ~vo 

/ 

A,g/h 

U4 

-/ 

u5 

-/ 

A/6  “ Ag€  - 0 


Note  that  only  the  first  two  components  are  affected  by  the  volume  frac- 
tions of  the  constituent  plies  in  the  laminate.  The  remaining  four  components 
are  constant  or  zero.  Cross -ply  laminates  are  orthotropic  because  the 
shear  coupling  terms  are  zero.  If  we  substitute  the  definitions  of  the  U. 
from  Eq.  99  into  the  last  two  equations  in  Table  37,  we  will  have: 


Vh  = U4 

’ U3 

= °12 

A66/h  * U5 

- U3 

(185) 


Following  are  some  numerical  examples  of  the  calculation  of  the 
In-plane  modulus  and  compliance  of  cross-ply  laminates.  T300/5208  will 
bo  used  as  our  sample  material.  The  elastic  modulus  in  terms  of  the 
U is  listed  in  Table  22.  Combining  the  modulus  data  with  the  formulas 


In  Table  37,  we  arrive  at  the  following  expressions: 


An/h  = 76.37  + (vQ  -v9Q)  85.73  + 19.  71 


A22/h 

= 76.37  - 

(vo  - 

v ) I 
90 

Al2/h 

= 22.61  - 

19.71 

= 2.90 

A66/h 

= 26.88  - 

19.71 

= 7.17 

> 

o 

II 

A26  = ° 

(186) 


The  results  from  these  equations  are  plotted  in  Fig.  52,  using  the  volume 
fraction  of  90 -degree  plies  as  the  abscissa.  Note  that  the  rule  of  mixtures 
relations  apply  in  the  first  two  components.  Both  are  linear.  The  other 
two  components  in  Eq.  186  are  constant  for  all  volume  fractions. 


100  Percent  90°  Plies 


Fig.  52.  In-plane  modulus  of  cross-ply  composites.  Note  all  lines  are 
straight  indicating  linear  relationship  as  functions  of  90-degree 
plies.  The  Poisson  and  shear  components  are  constant  and 
independent  of  volume  fractions. 


| 128 


We  can  also  calculate  the  compliance  components  by  inverting  the 
modulus  at  a given  volume  fraction.  The  inversion  process  though  must 
be  repeated  for  other  fractions  as  well.  Let  us  assume  that  the  volume 
fractions  are  50  percent.  The  modulus  components  are: 

Vh  * A22  /h  = 96.08,  Al2/h  = 2.90,  A^/h  = 7.17  (187) 

Using  the  matrix  inversion  method  described  in  Eq.  138  and  139,  we  have 

the  following  solutions  where  we  substitute  A../h  for  Q. .,  and  the  answer  is 

U U 

a,.h  in  place  of  S. 
ij  U 

detA..  = 66.126  x 103°h3(Fh)3 
ij  ' 


allh  = a22h  = 10‘41<TPa) 

a12h  =”*3141 

a66h=  139.44  » 

a16  = a26  = 0 


(188) 

, 


These  data,  and  similar  data  for  other  volume  fractions  of  90 -degree 
plies,  are  plotted  in  Fig.  53.  Note  that  the  rule  of  mixtures  or  linear 
relation  no  longer  exists.  From  the  compliance  components,  we  can 
get  the  following  engineering  constants  using  Eq.  157: 

E°  = 1/10.41  = 96.0  GPa 

^°2=  . 3141/10.41  = .0301  (189) 

G°  = 1/139  = 7. 17  GPa 
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Percent  90°  Plies 

Fig.  53,  In-plane  compliance  of  cross-ply  composites.  Shear  compliance 
is  constant.  Highly  nonlinear  relationship  exists  in  the  a^  and 
80  t*le  ru*e  °*  mixtures  equation  is  not  applicable.  The 
Poisson  compliance  is  also  nonlinear  but  very  small. 


Since  cross -ply  laminates  are  orthotropic,  we  could  have  calculated  and 
obtained  the  same  engineering  constants  in  Eq.  189  from  those  in-plane 
modulus  components  in  Eq.  187  directly  using  the  relations  in  Eq.  13. 
Q,.  A, , /h 


where 


m m 

l A11  A22j  [ 96.08ZJ 


96.08 

1.001 


96.  0 GPa 


This  agrees  with  the  longitudinal  modulus  in  Eq.  189. 


(190) 


(191) 


The  constraining  effect  of  the  90  degree  ply  is  responsible  for  the  low 
Poisson's  ratio  in  Eq.  189.  Since  the  value  of  m in  Eq.  190  is  almost  unity, 
therefore 

^/h  * E°  (192) 

This  is  true  only  for  cross-ply  composites.  Since  A^/h  follows  the  rule 
of  mixtures  relation,  we  can  say  that  E^  will  follow  approximately  the 
same  relation. 


5.  ANGLE -PLY  LAMINATES 

Angle -ply  laminates  form  another  very  common  group  that  deserves 
special  attention.  In  this  group,  there  are  only  two  ply  orientations  which 
have  the  same  magnitude  but  opposite  signs.  The  laminate  is  balanced, 
which  moans  that  there  are  an  equal  number  of  plies  in  both  the  positive 
and  negative  ply  orientations  or,  there  exists  a bisector. 


Thus  for  angle -ply  laminates  we  have: 

*1  * t0-  *2  * -« 


and  V1  = v2  = 

Substituting  these  values  into  Eq.  178  et  al. 


* 1 

Vj  = ^(cos  20  + cos20)  = cos20 

V = cos40 
2 


* 


V 


* 

4 


(193) 


(194) 
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**  » 
. v* 


* 


0 


The  formulas  for  the  in-plane  modulus  for  angle -ply  laminates  are  listed 
in  Table  38,  where  matrix  multiplication  is  implied.  They  are  obtained 
by  substituting  the  values  from  Eq.  194  into  Table  35. 


Note  that  the  first  four  rows  of  this  table  are  identical  to  those 
in  Table  20  for  the  unidirectional  modulus  Q„  except  where  the  ply  orien- 
tation 0 in  Table  20  must  be  replaced  by  the  angle  0 in  the  angle- 


ply  laminate.  The  shear  coupling  terms  vanish  for  the  angle -ply  laminate 
because  of  the  last  line  in  Eq.  192,  or  V^=V4=0  . 

Table  38.  FORMULAS  FOR  IN-PLANE  MODULUS  FOR 
ANGLE-PLY  LAMINATES 


/ 

A„/h 

u, 

cos2</> 

cos4<fi 

AM/h 

u, 

-cos2(fi 

cos4<fi 

Al3!/h 

U4 

-cos  4 

4 

-cos  4cfi 

A~  -A 

„=o 

Because  of  the  similarity  between  this  table  and  Table  20,  we  can  imme- 
diately convert  the  numerical  results  for  T300/5208  in  Fig.  35  by  simply 
replacing  9 by  ±0  and  deleting  the  shear  coupling  terms.  The  comparable 
components  of  in-plane  modulus  for  T300/5208  angle-ply  laminates  as 
functions  of  ±0  are  equal  to  those  in  the  first  four  columns  in  Table  22. 
The  last  two  columns  are  zero.  Like  cross-ply  laminates,  angle-ply 
laminates  are  also  orthotropic. 

For  [45/-45]  angle  ply,  we  have 

cos20  •*  0,  cos40  = -1 
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Using  the  data  for  T300/5208  from  Table  22,  we  have  from  the  formulas 
for  in-plane  modulus  in  Table  38 


A1 1 / h = A 22/h=  76.37  - 19.71  = 56.66  GPa 


A12/h  = 22.61  + 19.71  = 42.32 
A66/h  = 26.  88  + 19.  71  = 46.  59 


(195) 


A16  ■ A26  “ ° 

With  the  exception  of  these  shear  coupling  terms,  these  values  are  identical 
to  those  for  0 = 45  degrees  in  Table  22.  The  in-plane  modulus  of  angle-ply 
laminates  are  listed  in  Table  39  and  plotted  in  Fig.  54.  Using  the  inversion 
method  applied  in  Eq.  188,  we  have  for  the  [45/-45]  in  Eq.  195: 

*0 

det  A..  = 66.  126  x 10  h (Pa) 

U 


ailh  = a22h  = 39,91  (TPa) 


a!2h  = "29*81 


-1 


(196) 


*66h’21-46 
a16  ’ *26  ' 0 

These  compliance  values  together  with  those  for  other  values  of  0 are  shown  as 
solid  lines  in  Fig.  55.  The  dashed  lines  are  the  transformed  components 
of  the  compliance  of  unidirectional  composites  S. . taken  from  Fig.  40. 
Components  other  than  the  shear  coupling  components  between  A../h  and 
the  transformed  Q..  are  identical  as  shown  in  Fig.  54.  Complete  differences, 
however,  exist  between  a. .h  and  the  transformed  S...  The  solid  versus 

lj  ij 

dashed  lines  in  Fig.  55  show  these  differences. 


133 


• £* S h&M  ' ' - 


±0 

All/h 

A22/h 

A12/h 

A66/h 

A16  A26- 

0 

181.8 

10.  3 

2.90 

7.  17 

0 

15 

160.4 

11.9 

12.75 

17.  05 

0 

30 

109.3 

23.6 

32.46 

36.78 

0 ' 

45 

56.6 

56.6 

42.32 

46.  59 

o 

60 

23.6 

109.  3 

32.46 

36.78 

0 

75 

11.9 

160.4 

12.75 

17.05 

0 

90 

10.  3 

181.8 

2.90 

7.  17 

0 

Fig.  54.  In-plane  modulus  of  angle-ply  laminate  of  T300/5208  composite. 

With  the  exception  of  the  non-zero  shear  coupling  components, 
the  curves  above  are  identical  to  the  transformed  modulus  of 
unidirectional  T300/5108  shown  in  Fig.  35  where  the  coordinates 
are  defined  in  the  parenthesis.  These  curves  here  are  the  laminate 
modulus  as  a function  of  lamination  angle  0. 


135 


Fig.  55.  In-plane  compliance  of  T300/5208  angle-ply  laminates.  The 
solid  lines  are  based  on  the  data  in  Table  39,  and  are  not 
transformation  curves.  The  dashed  lines,  taken  from  Fig.  40, 
are  the  transformed  components  of  compliance  of  the  same 
unidirectional  composite,  with  the  coordinates  defined  in 
parenthesis;  i.  e. , 9 vs.  S^..  Close  similarity  between  A„/h 
and  Q'  does  not  exist  between  a. .h  and  S'... 


The  engineering  constants  for  off-axis  unidirectional  composites 
will  also  be  completely  different  from  those  for  angle- ply  laminates.  The 
engineering  constants  for  0 = 45  degrees,  which  is  simply  the  ±45-degree 
angle  ply,  can  be  obtained  directly  from  the  results  in  Eq.  196  using  the 
relations  in  Eq.  157* 

E°  = E°  = !/auh  = 39.91’1  = 25.05  GPa 

= — ai2 / au  = 29.81/39.91  = .746  (197) 

G12  = l/a66h  = 21. 46_1  = 46. 59  GPa 

The  corresponding  engineering  constants  for  an  off-axis  unidirectional 
T300/5208  were  calculated  from  the  data  in  Table  29: 

E2'  = 59.75"1  = 16.73  GPa 
9.99/59.75  = .167  (198) 

105.  7“  1 = 9.46  GPa 

Compare  like  constants  in  Eq.  197  and  198.  We  see  that  the  values  for 
angle -ply  laminates  in  Eq.  197  are  much  higher  than  the  off-axis  uni- 
directional in  Eq.  198*  The  Poisson's  ratios  in  Eq.  197  exceed  the  upper 
limit  for  ordinary  materials,  which  is  1/2.  This  is  theoretically  admissible 
for  non-isotropic  materials.  The  comparison  of  these  engineering  constants 
between  off-axis  unidirectional  and  angle-ply  laminates  as  functions  of  ply 
orientation  0 and  lamination  angle  ±0  are  shown  in  Fig.  56.  The  significant 
increase  in  the  angle-ply  laminates  over  that  of  the  off-axis  unidirectional 
over  the  entire  range  of  angles  is  very  apparent.  The  increase  is  caused 
by  the  constraining  influence  of  a laminate.  The  plies  are  coupled  and  are 
not  free  to  deform  independently. 


E1 


12 


G ' = 

12 


Poisson  s Ratio 


Youngs  Modulus,  GPa 


i \x-o;,(0) 


Shear  Modulus,  GPa 


Fig.  56.  Comparison  between  elastic  constants  of  angle  ply  and 
unidirectional  composites.  The  variations  of  elastic 
constants  are  shown  as  dashed  and  solid  lines.  The 


dashed  lines  are  identical  to  those  in  Fig.  44 


6. 


QUASI- ISOTROPIC  LAMINATES 

If  the  following  conditions  for  the  in- plane  modulus  of  a laminate 
are  satisfied,  the  laminate  is  quasi-isotropic: 

A11  = A22’  2A66  = A11  " A12  (199) 

The  last  relation  is  analogous  to  that  for  the  modulus  Q in  Eq.  23.  This 

ij 

will  become  evident  later.  The  two  relations  in  Eq.  199  reduce  the  number 
of  independent  components  of  the  in-plane  modulus  from  four  to  two.  This 
is  a necessary  condition  for  isotropy. 

Intuitively,  if  ply  orientations  are  random,  we  would  expect  an 
isotropic  laminate.  Directionality  would  disappear.  For  example,  chopped 
fiber  composites  are  quasi- isotropic.  If  we  examine  Eq.  170  which  defines 
the  Vj,  it  is  reasonable  to  expect  the  integrated  trigonometric  functions  to 
vanish.  Physically,  when  there  is  equal  probability  of  fibers  oriented  in 
any  direction,  or  there  is  a continuous  variation  in  fiber  orientation,  the 
cyclic  terms  in  Table  34  as  defined  by  the  V.  will  vanish.  The  in-plane 
modulus  components  will  converge  toward  the  invariants  in  the  first 
column  of  the  table.  The  in-plane  modulus  becomes: 


AU/h-  A22/h-  U1 


A12/h-  U4 


(200) 


/h  - U. 
5 


The  conditions  for  isotropy  in  Eq.  199  are  satisfied  because  of  the 
relations  between  the  invariants  as  described  in  Eq.  102.  We  have  a 
quasi-isotropic  material. 
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Since  this  is  an  isotropic  material,  we  can  find  the  quasi-isotropic 
engineering  constants,  as  follows: 


For  the  second  laminate,  we  have: 


v = V = v = V = — 

1 2 3 4 4 


(204) 


V*  = ^ £cos  0 + cos  180  + cos  90  + cos  (-90)J  = 0 
* 1 

^2/  = “fcos  0 + cos  360  + cos  180  + cos(-180)J  = 0 

V3  = ^ £sin  0 + sin  180  + sin  90  + sin  (-90)J  = 0 

V4  = sin  0 + sin  360  + sin  180  + sin(-180)J  = 0 


Since  all  the  V.  are  zero,  the  laminates  are  quasi-isotropic.  In  fact, 
we  can  generalize  that  any  laminate  with  "m"  ply  assemblies  spaced  at 
ply  orientations  of  "180/m"  degrees  will  be  quasi-isotropic.  In  the  first 
case  we  had  m = 3;  in  the  second  case,  m fc  4.  With  symmetric  laminates 
moreover,  we  must  double  the  number  of  ply  assemblies.  The  minimum 
number  of  plies  are  6 and  8,  respectively.  The  first  laminate  is  also 
called  Pi/3;  and  the  second  Pi/4. 

There  is  a very  practical  reason  for  quasi-isotropic  laminates  be- 
yond being  isotropic  like  ordinary  materials.  This  configuration 
represents  the  minimum  performance  that  we  can  expect.  Since  we  may 
be  uncomfortable  in  dealing  with  directionally  varying  properties,  we  can 
just  use  the  quasi-isotropic  laminate.  A direct  substitution  of  this  laminate 
for  the  ordinary  material  can  be  done  without  hesitation,  because  this 
substitution  is  no  different  from  substitution  of  any  other  ordinary  materials. 

The  qua  si -isotropic  Young's  modulus  of  T300/5208,  as  listed  in 
Eq.  202,  is  equal  to  the  Young's  modulus  of  aluminum.  But  there  is  a 
minimum  of  40  percent  savings  in  weight.  When  directionality  is  judiciously 
added,  the  advantages  of  composites  are  overwhelming. 
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The  quasi-isotropic  laminates  are  used  as  the  starting  point  of 
optimization  of  ply  orientation.  If  minimum  weight  is  a criterion,  the 
quasi- isotropic  laminate  should  be  the  upper  bound  of  the  weight.  An 
optimized  material  taking  full  advantage  of  the  directionality  of  properties 
should  only  have  lower  weight  than  the  qua  si- isotropic  configuration. 

7.  GENERALIZED  PI/4  LAMINATES 

This  is  a family  of  laminates  having  four  ply  orientations  spaced  at 
45-degree  intervals.  The  normal  Pi/4  laminates  have  four  ply  assemblies 
with  equal  thickness  and  are  therefore  quasi-isotropic.  Generalized  Pi/4 
laminates  refer  to  those  with  arbitrary  thicknesses  in  ply  assemblies, 
including  zero  thickness  for  one  or  more  ply  assemblies.  We  will  now 
list  all  the  trigonometric  functions  and  their  values  for  our  ply  orientations 
in  Table  40. 

Table  40.  VALUES  OF  TRIGONOMETRIC  FUNCTIONS  FOR  IN- PLANE 
MODULUS  OF  GENERALIZED  PI/4  LAMINATES 


•i 

cos 26^ 

co»40. 

sin26. 

sin40^ 

o 

1 

1 

0 

0 

90 

-1 

1 

0 

0 

45 

0 

-1 

1 

0 

-45 

0 

-1 

-1 

0 

* 

V4  = 0 


With  these  values,  the  formulas  for  in-plane  modulus  in  Table  35  can  be 
specialized  for  our  generalized  Pi/4  laminates.  This  is  done  in  a matrix 
multiplication  table  as  follows: 

Table  41.  FORMULAS  FOR  IN-PLANE  MODULUS  OF  GENERALIZED  i 


PI/4  LAMINATES 


Note  when  all  the  v^  are  equal,  we  recover  the  quasi-isotropic  laminates. 
When  the  45 -degree  plies  are  zero,  we  recover  the  formulas  for  cross - 
ply  laminates  listed  in  Table  37.  When  we  have  a special  angle -ply  laminate 
with  the  lamination  angle  0 equal  to  45  degrees,  we  recover  from  Table  41 
the  special  ±45  laminate  from  Table  38.  Finally,  Pi/4  laminates  are 
orthotropic  when  the  45 -degree  plies  are  balanced,  or  when 


(206) 


When  this  is  true,  the  shear  coupling  components  become  zero. 


The  formulas  in  Table  41  can  be  represented  by  a series  of  diagrams 

or  plots.  First  of  all,  the  components  of  modulus  are  all  linear  functions 

of  the  ply  fractions.  The  components  are  proportional  to  four  linear 

combinations  of  the  ply  fractions;  viz.,  v^,  v^q*  v^g  + v 45  ant* 

v^g  - v . In  Fig.  56  we  show  the  in-plane  modulus  of  generalized 

Pi/4  laminates  for  T300/5208  composite.  The  first  chart  shows  component 

A^/h.  This  chart  is  valid  for  balanced  as  well  as  unbalanced  laminates; 

i.  e. , independent  of  the  values  of  v and  v . The  same  can  be  applied 

45  ”45 

to  the  second  chart  on  components  A /h  and  A,  ,/h.  Component  A /h  is 

12  60  22 

not  shown  because  it  can  be  found  by  interchanging  v^  with  v from  the 
A^/h  chart.  Only  in  unbalanced  laminate,  that  is  the  +45  has  different 
number  of  plies  from  -45,  will  the  last  chart  in  Fig.  57  become  necessary. 

The  open  hexagon  in  each  diagram  represents  the  properties  of  the 
quasi -isotropic  laminate.  Note  all  relationships  in  Table  41  are  described 
by  straight  lines.  This  did  not  happen  by  accident.  In  fact,  it  is  important 
to  choose  the  correct  parameters  so  linear  relationships  exist  between  the 
ply  and  the  laminate  properties.  The  correct  property  set  for  the  case 
of  stiffness  of  the  laminate  is  the  components  of  modulus  of  the 
unidirectional  composites.  If  another  property  set  is  chosen,  nonlinear 
relationship  would  result.  When  the  property  set  consisting  of  engineering 
constants  are  chosen,  the  straight  lines  in  Fig.  57  will  be  replaced  by 
curved  lines.  In  fact  the  curved  lines  have  been  referred  to  by  a trade 


name:  the  carpet  plot.  The  moral  of  the  story  is  that  modulus  is 
the  simplest  property  for  the  description  of  the  stiffness  of  laminated 
composites.  The  carpet  plot  is  an  unnecessarily  complicated  way  of 
showing  properties  of  composites. 


Fig.  57.  In-plane  modulus  of  generalized  Pi/4  laminates  of  T 300/  5208 

composite.  Quasi-isotropic  points  are  shown  as  open  hexagons 


SCOPE 

The  flexural  stiffness  of  symmetric  sandwich  laminates  with  honeycomb 
core  and  multidirectional  composite  facing  will  be  covered.  The  strain  in 
the  laminate  is  defined  as  linearly  proportional  to  the  curvature.  Flexural 
rigidity  can  then  be  defined  as  the  modulus  for  flexural  stiffness  and  can  be 
related  to  the  bending -curvature  relation.  The  contribution  of  the  core  and 
the  effect  of  stacking  sequence  on  the  flexural  rigidity  will  also  be  discussed. 
The  key  relation  for  the  flexural  modulus  of  any  laminate  with  or  without 
core  is: 


I 


PRINCIPAL  NOMENCLATURE 


The  width  of  beams 

Flexural  modulus  of  multidirectional  laminates,  with  or 
without  core,  in  Nm;  i,  j- ■ 1,2,6. 

D../h*3  or  12D../h3  when  z = 0. 
ij  »J  c 

Flexural  compliance  of  multidirectional  laminates,  and  is 
the  inverse  of  D„,  in  (kNm)  S i,j  = 1,2,6. 

A flexural  engineering  constant:  the  effective  Young's  modulus 
along  the  1-axis,  in  GPa. 

A flexural  engineering  constant:  the  effective  Young's  modulus 
along  the  2-axis,  in  GPa. 

A flexural  engineering  constant:  the  effective  longitudinal 
shear  modulus,  in  GPa. 

Total  thickness  of  laminate,  with  or  without  core,  in  m. 
h3[l  - z*3J/ 12 

n.h  » Total  thickness  of  the  i-th  ply  assembly;  i = 1,  2, . . .m 
i o 

(h  - 2z  )/n  = Unit  ply  thickness 
c 

Momen  of  inertia 
Curvature,  in  m Si*  1,2,6. 

Total  number  of  ply  assemblies  in  a laminate. 

Moment,  in  Nm  Si*  1,  2,  6. 

Total  number  of  plies  in  a laminate. 

Total  number  of  plies  in  the  i-th  ply  assembly;  i = 1,  2, ...  m. 

Modulus  of  the  i-th  ply  assembly,  in  GPa;  i,  j • 1,  2,  6. 

Modulus  of  the  ply  assembly  with  6 orientation;  i,  j * 1,  2,6. 

Linear  combinations  of  modulus  for  the  multiple -angle 
transformation  of  modulus,  in  GPa;  i = 1,  2,  3,4,  5. 


1.  LAMINATE  CODE 

The  same  laminate  code  convention  as  that  used  for  the  in-plane 
modulus  in  Eq.  140  will  be  followed  for  the  flexural  modulus.  For  symmetric 
laminates  we  can  add  the  half  depth  of  the  core  in  the  code;  for  example: 

|03/902/45/-453/zcJ  g (207) 

The  orientations,  ply  assemblies  and  the  core  for  this  laminate  are  shown  in 
Fig.  58.  The  plies  are  arranged  in  an  ascending  order  from  the  bottom  or 
the  z = -h/2  face.  This  again  can  be  a source  of  confusion.  The  code  in 
Eq.  207  applies  to  the  lower  half  of  a symmetric  laminate  starting  from 
the  bottom  face.  The  stacking  sequence  in  the  upper  half  of  the  laminate 
is  in  reverse  order  of  the  code  shown  in  Eq.  207.  The  actual  integration 
for  the  calculation  of  the  flexural  modulus  is  applied  over  the  upper  half 
of  the  laminate  which  extends  from  z * 0 to  z = h/2. 


In  the  case  of  the  in-plane  modulus,  only  the  volume  fractions  of  the 
ply  assemblies  are  important.  This  is  clearly  shown  in  Eqs.  153  and  170.  The 
actual  stacking  sequence  does  not  affect  the  in-plane  modulus.  Whether  the 
laminate  code  is  intended  to  follow  an  ascending  or  descending  order  is  of 
no  consequence  to  the  in-plane  modulus.  This,  however,  is  no  longer  true 
for  the  flexural  modulus  that  we  will  discuss  in  the  section.  The  positions 
of  ply  assemblies  in  a laminate  have  direct  effect  on  the  flexural  modulus. 

That  is  why  we  are  discussing  the  laminate  code  again. 

2.  MOMENT-CURVATURE  RELATIONS 

In  the  flexural  behavior  of  laminates,  moment  and  curvature  are  the 
key  variables,  similar  to  the  role  of  stress  resultant  and  in-plane  strain 
in  the  in-plane  behavior  of  the  last  section.  The  counterpart  of  the  stress - 
strain  relation  for  the  in-plane  behavior  is  the  moment-curvature  relation 
for  the  flexural  behavior.  The  elastic  constants  for  the  latter  relation  will 
be  called  the  flexural  modulus  and  flexural  compliance.  It  is  the  purpose 
of  this  section  to  develop  definitions  of  moment  and  curvature,  and  their 
relationship  to  each  other. 

The  distribution  of  ply  stresses  can  be  symmetric  and  antisymmetric 
with  respect  to  the  midplane.  In  Section  4,  the  stress  distribution  was  sym- 
metric and  this  was  shown  in  Fig.  48.  In  Fig.  59  we  will  repeat  the  sym- 
metric distribution  of  Fig.  48,  and  will  also  show  ths  case  of  antisymmetric 
distribution. 

As  the  result  of  symmetric  stress  distribution  in  Fig.  59(a),  we 
can  represent  the  variable  stress  by  an  average  stress  (J^  or  stress 
resultant  N^.  The  in-plane  behavior  of  symmetric  laminates  can  be 
completely  characterized  using  the  average  stress  or  stress  resultant. 

When  the  stress  distribution  is  anti- symmetric,  as  shown  in  Fig.  59(b), 
the  average  stress  across  the  entire  laminate  thickness  is  zero.  One 
approach  of  dealing  with  the  anti- symmetric  stress  distribution  is  to 
define  a new  quantity:  the  moment,  to  take  the  place  of  stress  resultant. 

The  simplest  or  first  moment  has  three  components: 
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The  unit  of  moment  is  N,  or  Nm/m;  i.  e 


a moment  per  unit  width  of  a 


laminate  with  thickness  h, 


The  sign  of  the  components  of  moment  is  also  critical.  The  bending 
components  of  moment,  like  the  normal  components  of  stress  and  strain, 
are  easy  to  rationalize  and  readily  defined.  A bending  moment  is  positive 
if  the  average  induced  stress  in  the  upper  half  of  the  laminate  is  positive. 

In  Fig.  60(a)  we  define  the  positive  component  for  in  Fig.  60(b),  the 
positive  M^.  When  or  M^  is  negative,  the  average  induced  stress  in  the 
upper  half  of  the  laminate  will  be  negative.  We  use  average  stress  here 
because  in  a laminated  material  it  is  possible  to  have  both  positive  and  nega- 
tive stresses  in  each  half  of  the  laminate.  Fig.  59(b)  shows  this  possibility. 


(a)  (b)  (c) 

Fig.  60.  The  positive  directions  of  components  of  moment. 

Bending  moments  are  shown  in  (a)  and  (b). 

In  (c),  the  positive  twisting  moment  is  shown  and  the  induced 
shear  stresses  are  idealized  as  linear  in  the  z coordinate. 

Note  positive  twisting  moment  appears  as  clockwise  torque  on 
the  positive  1-axis  face;  counterclockwise  on  the  positive  2 -axis 
face.  The  right-hand  rule  is  not  invoked  here.  The  effect  of  the 


positive  twisting  moment  can  be  duplicated  by  four  self-equili- 
brating  forces  acting  at  the  corners  as  shown.  These  arrows  are 
useful  in  determining  the  direction  of  twisting  as  illustrated  in 
Fig.  25  by  Timoshenko  and  Woinowsky-Krieger*. 


*S.  Timoshenko  and  S.  Woinowsky-Krieger,  Theory  of  Plates 
and  Shells.  McGraw-Hill,  1959.  pp.  80-81 


The  sign  convention  for  twisting  moment  follows  the  same  rule; 
viz. , a positive  shear  stress  on  the  upper  half  of  the  laminate  is  associated 
with  the  positive  twisting  moment.  The  positive  shear  stress  component  is 
defined  in  Fig.  6.  Fig.  60(c)  shows  the  result  of  positive  twisting  moment 
and  the  induced  shear  stress  distribution.  All  the  arrows  will  reverse  their 
directions  if  the  twisting  moment  is  negative.  We  are  not  imposing  the  right- 
hand  rule  for  the  sign  convention  except  that  the  coordinates  and  angle  of  rota- 
tion follow  this  rule.  If  the  right-hand  rule  is  followed,  as  shown  by  Timoshenko, 
we  must  distinguish  the  twisting  moment  on  the  1-axis  face  as  M^,  an<*  *^at  on 

the  2 -axis  as  M,,;  or  M and  M , respectively. 

21  xy  yx 

Then  we  have  the  following  relations: 


Therefore 


M._  = 

M 

12 

xy 

M,,  = 

M 

21 

yx 

M.„  = 

12 

21 

(209) 


(210) 


The  important  issue  here  is  not  what  sign  convention  we  use.  We  must  under- 
stand the  rationale  and  be  consistent.  Again  we  would  like  to  mention  how 
critical  signs  are  when  we  work  with  composite  materials.  A wrong  guess 
is  often  inconsequential  for  ordinary  materials,  but  can  be  disastrous  for 
composites.  The  signs  for  shear  stress,  shear  strain,  twisting  moment  here 
and  twisting  curvature,  which  we  will  introduce  presently,  are  all  sources  of 
uncertainty  and  error. 

We  will  now  derive  the  strain -displacement  relation  for  the  bending  of 
a plate  similar  to  that  for  the  in -plane  stretching  of  a plate  in  Eq.  1 and  4. 


We  will  assume  that  the  plate  is  initially  flat  as  shown  in  Fig.  61(a).  After 
bending,  the  plate  can  be  described  by  a function  w where: 

w = w(x,  y)  (211) 

It  is  implied  that  the  vertical  displacement  of  each  point  does  not  vary  in 
the  z -direction.  The  normal  to  the  plate  does  not  stretch  or  deform.  It 
only  rotates  as  the  plate  is  bent  or  twisted.  Fig.  61(b)  is  an  illustration 
of  a bent  plate. 


Fig.  61.  Definition  of  a plate  or  laminate  before  and  after  bending. 

The  deformed  midplane  is  described  by  a function  w(x,  y). 

The  rotation  of  the  normal  to  the  midplane  can  be  directly  related  to 
the  first  derivative  at  the  same  point  in  the  plate.  This  is  shown  in  Fig.  62, 
where  two  cases  of  the  bent  plane  is  shown  for  the  purpose  of  establishing 
the  correct  sign  convention.  Consistency  between  Fig.  60  and  61  is 
maintained  if  we  recognize  that  a change  in  sign  in  the  derivation  is  necessary. 


i 


If 


(a)  (b) 


Tension  Compression 

at  z>o  at  z>o 


Fig.  62.  Sign  convention  of  midplane  displacements.  For  a concave 
downward  deformation  in  (a),  the  derivative  of  w is  negative, 
and  a negative  sign  must  be  added  to  the  displacement-derivative 
relation  in  Eq.  212.  When  the  curvature  is  reversed  in  (b), 
the  derivative  of  w is  now  positive. 


The  definition  of  displacement-derivative  relation  is: 


u * -z8 


Similarly,  we  can  derive  the  displacement  along  the  y-axis  as: 


v * 


(212) 


(213) 


" •>.-  , . 
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From  the  last  two  equations,  and  the  strain-displacement  relations  of 
Eq.  1 and  4,  we  can  show  that: 


We  can  relate  the  second  derivatives  to  curvatures  k,  as  follows 


The  bending  curvatures  are  the  reciprocals  of  the  radii  of  curvature.  The 
relations  in  Eq.  215  can  be  found  in  elementary  calculus,  and  can  be  related 
to  simple  graphical  illustrations.  The  sign  convention  used  here  is  based 
on  that  established  in  Fig.  60.  The  twisting  curvature  is  difficult  to  illus- 
trate and  is  not  normally  covered  in  elementary  text.  We  derived  our  rela- 
tion through  the  use  of  the  strain -displacement  relation  in  Eq.  4.  Substituting 
the  definitions  in  Eq.  215  into  214,  we  have: 


This  assumed  linear  strain  distribution  is  shown  in  Fig.  6 3.  A more  general 
assumed  state  of  strain  than  both  Eq.  145  and  216  would  be  the  sum  of  the  two. 
This  combined  strain  will  be  used  as  the  basis  of  general,  nonsymmetrical 
laminates  which  we  will  cover  in  the  future. 


Z Z 


Fig.  63.  Assumed  linear  strain  distribution  across  laminate  thickness. 

Maximum  strain  values  reached  at  the  upper  and  lower  faces. 

They  are  equal  but  opposite  in  signs  when  the  laminate  is 
symmetric. 

We  can  now  derive  the  moment -curvature  relations  by  substituting 
the  assumed  strain  into  the  definition  of  moment  in  Eq.  208.  We  must  first, 
however,  substitute  the  off-axis  stress -strain  relations  listed  in  Table  18  or  30 
into  Eq.  208.  This  will  express  the  stress  components  in  terms  of  the 
strain  components.  As  we  did  in  the  derivation  of  the  in-plane  stress - 
strain  relations,  we  deleted  the  primes  from  ail  the  relations  in  Table  18 
for  the  sake  of  simplicity; that  was  also  done  in  Table  30. 
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From  Eq.  208 


From  Table  30 
From  Eq.  216 


M1  = J Ojzdz 

= /[Quei  + Q12€2  * Q1666]ldZ 


■/[V.  + Q12k2  + Q16k6]  A* 


(217) 

(218) 


Since  k are  constant,  not  dependent  on  z,  they  can  be  factored  out, 
i 


mi  {/v2<i']ki  +[/Qi2*2dz]k2 


Ml=  Vl  + D12k2  + D!6k6 


Similarly: 


M, 


D12kl  * °22k2  * °26k6 


(219) 


M, 


D16kl  + D26k2  + °66k6 


where 


Dn  = /v2dz*  D22  -f 


Q22Z  d2 


D12  =/°12Z  dZ'  D66  =/*Q"z_dz 


66 


(220) 


D16  =/  °16 


2 

z dz. 


°26  -fQ 


26 


2, 

z dz 


We  have  thus  derived  the  moment -curvature  relations  in  Eq.  219  and  defined 
the  flexural  modulus  in  Eq.  220. 


Now  by  inverting  the  moment-curvature  relations  we  can  obtain  the 
following  relations  in  terms  of  flexural  compliance,  duplicating  the  same 
steps  used  in  the  inversion  of  Eq.  138  and  139. 


k.  = dUMl  ♦ d12M2  * d16M6 


d12Ml  + d22M2  + d26M6 


(221) 


k6  ’ d16Ml  f * d"M' 


' - •*  * v 


The  relationships  above  can  be  presented  in  matrix  multiplication  tables  as 
follow  8: 


Table  42.  MOMENT -CURVATURE  RELATION  OF  SYMMETRIC 
LAMINATES  IN  TERMS  OF  MODULUS 


*, 

kg 

ke 

M, 

On 

0,2 

D,e 

Mg 

0,2 

022 

0^6 

Me 

D,e 

02S 

Dee 

Table  43.  MOMENT-CURVATURE  RELATION  OF  SYMMETRIC 
LAMINATES  IN  TERMS  OF  COMPLIANCE 


M, 

M2 

Me 

K 

d„ 

d/2 

die 

k2 

d/2 

d22 

dee 

ke 

die 

d2e 

dee 

I 


We  can  now  define  the  effective  flexural  engineering  constants.  From 
the  compliance  relations  in  Table  43,  we  know  that  under  simple  bending  of 
only,  the  resulting  curvature  along  the  1-axis  is: 


kl  = duMi  (222) 

The  rigidity  of  this  material  is: 

M1  1 

IT  = T"  <223> 

i u 

From  elementary  theory,  we  know  that  rigidity  of  a homogeneous  beam  is 

Rigidity  = El  (224) 

where  E is  the  homogeneous  Young's  modulus;  and  I is  the  moment  of  inertia. 
By  equating  the  two  relationships,  we  have: 


1 1 
k,  “ d 


(223) 


Similarly,  we  can  show: 


El  = 

l/dll 

(225) 

II 

1/Idn  = 12/h3du 

(226) 

= 1. 

I = h3/12 

Ef  _ 
2 

12/h3d22 

1 

f 

v = 

12 

-d  /d„ 

12  11 

> (227) 

Gf2  = 

12/h3d., 

66 

The  superscript  f denotes  effective  flexural  engineering  constants.  These 
are  the  constants  if  the  beam  or  plate  of  our  multidirectional  laminates  is 
treated  like  a homogeneous  material. 


■ : ..  - V.  - 

- .. 


i-  - w-i  ■ 'unwn 


3.  EVALUATION  OF  FLEXURAL  MODULUS 

We  will  now  evaluate  the  components  of  flexural  modulus  by  performing 
the  integration  of  the  components  in  Eq.  220.  Similar  to  the  case  of  in-plane 
modulus  in  Section  4,  we  will  first  substitute  the  off-axis  modulus  of  the  uni- 
directional composites  using  the  multiple -angle  transformation  relations 
listed  in  Table  20.  The  modulus  components  Q..  in  the  integrands  of  Eq.  220 
are  the  off-axis  components  except  that  primes  have  been  deleted.  Using 
the  unprimed  quantities,  we  have 


2, 

z dz 


Dn  ' fQ u 

-f[vi  * u 


From  Fq.  220 

From  Table  20  = / |U,  + U^cos20  + U2Cos4oJ  z^dz 

Since  U.  are  independent  of  z for  a laminate  with  the  same  composite. 


= U^fz^dz  + U ^jcoelQ  z^dz  + 


u3/‘ 


- u, h*3  f + u3v2 


cos46  z dz  (228) 

(229) 


where:  h 


*3 


/*h/2  2 | 

fh/2  2 

/•h/2 

= / z dz  = 2/ 

z dz 

= 2/ 

./-h/2  7 

0 

y z 

c 

- h3 17i  M 

3 - i3r> 

l - Z*31 

12  [\1-h/2/. 

- 12  L1 

L * 1 

c J 

2, 

z dz 


2 r o h/2 

■ T 1 

J z 


(230) 


* Volume  fraction  of  core  * 2z  /h 


, -.r 

1 ./-h/2 

’•■f. 


cos26  z dz 


h/2 

/2 


cos 28  z dz 


cos46  z dz 


(231) 
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It  is  assumed  that  the  honeycomb  core  has  no  stiffness  in  the  1-2  coordinate 
system.  That  is  the  reason  why  the  lower  limit  of  integration  is  set  at  the  half 


depth  of  the  core,  z . 

c 


Similarly, 


°22  * U1  h ' U2V1  + U3V2 


(232) 


°12  U4  h ‘ U3V2 


(233) 


°66  = U5  h - U3V2 


(234) 


°16  = ~2  U2V3  + U3V4 


(235) 


Kh> 


where 


°26  * T U2V3  - U3V4 


< 


sin20  z dz 


sin40  z dz 


(236) 


(237) 


Here  again,  the  evaluation  of  the  flexural  modulus  reduces  to  the  evaluation  of 
the  four  weighted  trigonometric  integrals,  V to  V^.  Similar  to  Eq.  170, 
we  can  combine  the  definitions  of  V.  into  one  expression: 


1,2,3,  4]  = 2J 2 


[coa29,  cos46,  sin29,  sin49]z  dz  (238) 


We  can  also  put  all  the  formulas  for  the  components  of  the  flexural  modulus 
into  a matrix  multiplication  table  as  in  Table  44.  Note  the  similarity  between 
this  table  and  the  formulas  for  in-plane  modulus  in  Table  34.  The  definitions 
of  the  V\,  however,  are  different. 


■ 

^ •*  , % .V  » 

, **  * > - 4 


Table  44.  FORMULAS  FOR  FLEXURAL  MODULUS  OF 
SYMMETRIC  SANDWICH  LAMINATES 


h*3 

Ug 

us 

@n 

u, 

K 

Vg 

@gg 

U, 

-K 

Vg 

D/g 

u4 

-Vg 

@66 

U5 

-Vg 

@ 16 

V3/2 

v4 

@26 

V,/2 

-V4 

where  h*3 - (I~z*3)h3/I2 


Fig.  64.  Schematic  diagram  of  a symmetric  sandwich  laminate.  There 
are  m ply  assemblies  and  n plies  in  the  laminate  using  indices 
i and  t,  respectively.  Assuming  the  half  depth  of  the  core  is 
equal  to  a multiple  of  unit  plies,  the  half  depth  can  be  designated 
by  i ■ c or  n£  » 6 in  this  figure. 
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Let  us  try  to  evaluate  the  first  term  in  Eq..  238. 

■h/2 


V1  = 2 


/: 


cos20  z dz 


(239) 


If  we  assume  that  each  ply  assembly  would  have  the  same  unidirectional 
material,  the  integration  in  Eq.  239  can  be  replaced  by  a summation. 

See  Fig.  64  for  the  definitions  of  indices  of  summation. 

m/2 

— I s .5  1 

(240) 


v = -r-  cos20.  |z.  - z.  1 

1 3 ifei  • L 1 -'J 


Simplification  of  this  summation  in  terms  of  volume  fractions  such  as  that 
for  the  in-plane  modulus  in  Eq.  178  is  not  possible  because  of  the  cubic 
terms  here.  Some  simplification  is  possible  if  the  half  thickness  of  the 


Substituting  Eq.  243: 


where  n equals  the  total  number  of  plies  including  the  core  thickness 

expressed  in  equivalent  number  of  plies.  The  variables  in  the  bracket 

can  be  expressed  in  terms  of  plies  using  Eq.  242.  The  formulas  for 

the  other  three  V will  take  the  same  form.  Only  the  trigonometric 
i 

function  changes;  i.  e. , cos29^  in  Eq.  245  is  replaced  by  sin29.,  and  others. 

The  bracketed  quantity  in  Eq.  245  is  therefore  a weighting  factor. 

In  the  case  of  the  in-plane  modulus,  the  weighting  factor  was  the  volume 
fraction  of  each  ply  orientation.  We  had  the  rule  of  mixtures  relation. 

In  the  case  of  flexural  modulus,  this  weighting  factor  put  heavier  em- 
phasis on  the  outer  plies.  Again,  if  we  assume  that  all  plies  have  the 
same  thickness,  and  the  same  core  half-depth  is  a multiple  of  the  unit 
ply  thickness,  we  can  establish  the  numerical  values  of  this  weighting 
factor  starting  with  the  midplane  as  zero  and  move  upward  toward  the 
top  surface  where  n/2-th  ply  is  located.  The  evaluation  of  this 
weighting  factor  is.  listed  in  Table  4 5.  Eq.  245  can  be  rewritten  as  follows: 

n/2 


v * = 4 j cos29t  [t3  - (t-1)3]  (246) 


Ply  Order 
Number,  t 

t 

t-1 

t3  - (t-1)3 

l 

l 

0 

1 

2 

2 

1 

7 

3 

3 

2 

19 

4 

4 

3 

37 

5 

5 

4 

61 

6 

6 

5 

91 

7 

7 

6 
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8 

8 

7 

169 
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9 

8 

217 

10 

10 

9 

271 

11 

11 

10 

331 

12 

12 

11 

397 

The  index  t is  used  here  to  distinguish  from  the  index  i in  Eq.  245.  The 
latter  index  is  intended  for  the  number  ply  assemblies;  and  the  former  index 
the  total  number  of  all  plies.  The  two  indices  will  be  equal  if  each  ply 
assembly  has  only  one  ply. 
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Using  the  numerical  values  listed  in  Table  45,  Eq.  246  can  now  be 
written  as: 


* 8 
V1  = T 

n 


Jcos201  + 7cos262  + 18cos203  + 37cos204  + . . .J  (247) 


If  adjacent  plies  have  the  same  ply  orientation,  we  have  for  example,  ply 
assemblies  with  two  plies  each. 


91  = 02  ' °3  = V ‘ * 


Then 


T = I"8cos20  + 56cos20  + 152cos20  + . . . 1 

1 3 L 1 3 5 J 


If  we  have  a honeycomb  core  with  a half-depth  of  4-ply  thickness,  or 


z = 4h 
c o 


(248) 


(249) 


the  first  ply  or  ply  assembly  for  the  facing  will  start  with  t = 5 in  Table  45. 

V*  = —z  |61cos201  + 91cos202  + 127cos203  + I69cos204  + . . .J 
n (250) 

where  the  total  number  of  plies  n must  include  the  half-depth  of  the  core 
which  is  equal  to  4 plies.  If  we  have  a 3 -ply  facing,  the  value  of  n/2  is  7. 


4.  FLEXURAL  BEHAVIOR  OF  UNIDIRECTIONAL  LAMINATES 

If  our  laminate  is  unidirectional,  the  ply  orientation  0 is  a fixed  value, 
independent  of  the  z coordinate.  The  trigonometric  functions  in  Eq.  238 
can  be  taken  outside  of  the  integrals.  The  resulting  V are: 

jjg  J 

Vr,  _ _ -i = fcos20,  cos40,  sin20,  sin40]  h (251) 

11*  2,  3,  4J 

where  h ^ is  defined  in  Eq.  230.  For  this  simplified  case,  the  formulas 
for  the  flexural  modulus  in  place  of  Table  44  can  be  restructured  as  follows: 


Table  46.  FORMULAS  FOR  THE  FLEXURAL  MODULUS 
OF  UNIDIRECTIONAL  COMPOSITES 


Note  that  the  constants  in  this  table  are  identical  to  those  of  the  trans- 
formed in-plane  modulus  of  unidirectional  composites  in  Table  20.  The  only 

$ 

difference  is  the  normalizing  factor  needed  for  the  flexural  modulus.  The  h 
factor  is  defined  in  Eq.  230.  Thus,  we  can  obtain  the  off-axis  flexural 
modulus  by  multiplying  the  off-axis  Q„  by  this  factor.  The  transformed 
flexural  modulus  will  have  exactly  the  same  as  the  transformed  in-plane 
modulus  in  Fig.  35.  This  is  shown  in  Fig.  65.  All  the  remarks  about  the 
transformed  in -plane  modulus  following  Fig.  35  are  equally  applicable  to 
the  flexural  modulus. 


Fig.  65.  Transformed  flexural  modulus  of  unidirectional  T300/5208.  These 
are  the  same  curves  as  those  in  Fig.  35  with  the  exception  of  the 
normalizing  factor  in  Table  46.  The  normalized  modulus  are 

3 

defined:  D*  = D /h*  , where  h*  is  defined  in  Eq.  230. 


By  comparing  Tables  46  and  20,  we  have: 


(252) 


Identical  factor  shall  be  applied  to  the  other  components  of  the  flexural 
modulus  D... 

From  this  simple  relation  in  Eq.  252,  we  can  immediately  write  down  the 
flexural  compliance  S_  divided  by  the  normalizing  factor.  We  have  now: 


(253) 


where  the  transformed  in-plane  compliance  can  be  found  from  Table  29. 
and  Fig.  41.  The  latter  is  repeated  in  Fig.  65  where  the  normalizing  factor 
has  been  added.  Thus,  the  off-axis  flexural  compliance  is  the  same  as  that 
of  the  in-plane  compliance  with  the  exception  of  the  h ^ factor.  So  long  as 
a beam  and  plate  consist  of  symmetric,  homogeneous  factings,  their 
flexural  modulus  and  compliance  can  be  obtained  directly  from  the  in-plane 
modulus  and  compliance,  respectively.  We  only  need  to  know  the  values 
of  the  normalizing  factor,  as  shown  in  Eq.  252  and  253.  The  flexural 
rigidity  of  unidirectional  composites  is,  therefore,  as  simple  to  calculate 
as  that  of  ordinary  materials. 
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As  we  have  defined  the  flexural  rigidity  in  Eq.  223  et  al,  we  know  that  a 
symmetric  sandwich  beam  with  unidirectional  facings  will  have  the  following 
rigidity: 


h3  l-z*3| 


(EDi  = — = 

11 


(254) 


where  (EI^  = bending  rigidity  of  a beam  cut  along  the  1-axis  of  our  laminate 
b = width  of  beam 

We  can  make  the  following  remarks  about  the  flexural  rigidity  of  beams  and 
plates  relying  on  the  expressions  in  Eq.  252  and  253,  respectively: 

• Rigidity  is  highly  dependent  on  the  thickness  h.  If  we  double  the 
thickness,  we  will  get  a cubic  increase  in  return,  or  8 times  the 
rigidity. 

• Removal  of  materials  near  the  midplane  is  a very  effective  way 
of  reducing  the  weight  without  much  sacrifice  in  the  rigidity. 

If  one-third  of  the  material  at  the  center  is  removed;  i.  e. , 
zc  = 1/3,  the  loss  in  rigidity  as  measured  by  z^  is  only  1/27 
of  the  solid  beam  or  plate. 

Both  remarks  are  valid  for  a composite  and  ordinary  materials  so  long  as 
the  facing  material  is  homogeneous.  If  multidirectional  composites  are 
used  for  the  facing,  the  remarks  above  are  only  true  qualitatively.  We 
will  discuss  this  later  in  this  section. 

For  an  off-axis  unidirectional  composite  facing,  the  beam  will  twist 
under  pure  bending.  This  is  the  equivalent  of  the  shear  coupling  in  the  in- 
pl&ne  behavior  of  off-axis  materials.  The  nonzero  shear  coupling  compo- 
nents and  Q26  will  result  in  nonzero  D^,  D^,  d^,  and  d^.  From 
Table  43,  we  can  relate  the  induced  curvatures  by  bending  moments.  For 
exampl«»  if  we  apply  bending  to  our  off-axis  beam  as  shown  in  Fig.  67, 


from  Table  43: 


dllMl 


d M 
12 


d16M 


1 

1 


(255) 


The  first  curvature  is  due  to  normal  bending;  the  second,  the  Poisson's 
coupling;  and  the  third,  the  twisting  coupling.  The  question  now  is  how 
will  the  twist  occur:  how  much,  and  in  what  direction.  This  is  the  recurring 
question  associated  with  shear  stress  and  shear  strain.  Again,  we  must  pay 
attention  to  the  sign  convention. 


/ 

4*  M,  'W'M, 

Fig.  67.  Pure  bending  of  an  off-axis  beam.  Positive  ply  orientation 

and  positive  moment  are  shown.  Heavy  arrows  show  the  direc- 
tion of  movements  of  the  four  corners,  similar  to  Fig.  60(c). 

We  know  from  Fig.  66  that  the  shear  coupling  terms  for  T300/5208  and  for 
most  practical  composites  are  negative  for  positive  ply  angles.  Since  the 
moment  in  Fig.  67  is  also  positive,  we  know  from  Eq.  255  that  the  twisting 

curvature  must  be  negative.  Now  refer  to  Fig.  60(c)  where  we  showed  the 
effect  of  a positive  twisting  moment  on  the  stress  distribution  and  the  possi- 
ble directions  of  displacements  indicated  by  heavy  arrows.  Hence  a posi- 
tive curvature  will  be  a clockwise  rotation  about  the  1-axis.  For  our  beam 
in  Fig.  67,  we  have  negative  curvature.  Therefore,  the  twisting  curvature 
caused  by  the  bending  moment  will  be  a counterclockwise  rotation  along 
the  l-axis.  This  rotation  is  represented  by  the  heavy  arrows  shown  at  the 
corners. 
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5. 


FLEXURAL  MODULUS  OF  CROSS-PLY  LAMINATES 
Cross  -ply  laminates  are  the  simplest  multidirectional  laminates. 
Repeating  the  values  of  the  trigonometric  functions  in  Table  36,  we  have 
the  following: 


Table  47.  VALUES  OF  TRIGONOMETRIC  FUNCTIONS 
FOR  CROSS-PLY  LAMINATES 


9. 

cos20. 

cos40. 

sin20. 

sin40. 

l 

l 

l 

i 

l 

0 

1 

1 

0 

0 

90 

-1 

1 

0 

0 

Let  us  study  the  effect  of  stacking  sequence  on  the  flexural  modulus 
of  symmetric  laminates.  We  will  use  a 16-ply  laminate  with  three  different 
stacking  sequences  as  shown  in  Fig.  68. 


Fig.  68.  Cross -ply  laminates  with  16  plies  but  different  number  of  ply 
assemblies;  viz.,  m=4,  8 and  16. 
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From  the  second  column  of  Table  47,  we  know  that 


cos4d.  = cos40  = 1 (256) 

A Ct 

$ 

Following  the  pattern  of  Eq.  247  for  V we  can  immediately  write  down  the 
analogous  relation  for  V^. 

8 512 

V*  = — (1  + 7 + 19  + 37  + 61  + 91  + 127  + 169)  = = 1 (257) 

2 n 512 

where  n = 16  was  used.  Because  of  the  special  relation  in  Eq.  256,  the 
V will  remain  constant,  independent  of  the  stacking  sequences  shown  in 

Ct 

Fig.  68. 


Knowing  the  values  of  cos26.  from  the  first  column  of  Table  47,  we 
can  substitute  the  values  into  Eq.  247  for  the  case  of  m=4  or  Fig.  68(a). 


V*  = 5J2  " 7 - 19  " 37  + 61  + 91  + 127  + 169)  = cr^=  T (258> 


512 


Note  that  the  first  ply  from  the  midplane  upward  is  a 90-degree  ply.  We 
have  mentioned  before  that  there  is  a difference  between  the  laminate  code 
as  defined  by  Eq.  207  which  follows  an  ascending  order  from  the  bottom  sur- 
face of  the  laminate.  The  stacking  sequence  starting  from  the  midplane  is  the 
opposite  of  the  laminate  code  for  all  symmetric  laminates,  as  shown  in  Fig.  68, 
Note  also  that  when  the  ply  orientation  for  the  entire  laminate  is  the  same, 
will  be  the  same  as  the  relationship  shown  in  Eq.  257,  or  will  be 
unity.  This  is  the  case  of  a unidirectional  composite  shown  in  Table  46. 

$ 

The  values  of  for  the  cases  of  m= 8 shown  in  Fig.  68(b),  and  that 
of  m=16  in  Fig.  68(c)  are  listed  below,  respectively: 


For  m=8 


V1  = 5l2  (_1  - 7 + 19  + 37  -61-91  + 127  + 169)  =||!  = -|-  (259) 
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For  m=16. 


V*  =—  (-1  + 7 - 19  + 37  - 61  + 91  - 127  + 169)  = — = — 
512  512  16 


(260) 


It  appears  that  a pattern  has  been  established  that  for  cross -ply  symmetric 
laminates  with  increase  ply  assemblies  as  shown  in  Fig.  68. 


V.  — “*  i m - 4,  6,  12,  . . . 
i m 


(261) 


where  m = total  number  of  ply  assemblies.  As  m increases,  approaches  zero; 


i.e. , the  cross -ply  laminate  becomes  square -symmetric,  or  D = D . 

11  £»£• 


Summarizing  the  results  for  this  family  of  cross -ply  laminates  in  which 
the  total  number  of  ply  assemblies  is  a variable,  we  can  enter  the  values  of 


V.  into  Table  44  and  arrive  at  Table  48.  Care  must  be  exercised  in  the 

l 


proper  use  of  normalizing  factors. 


Table  48.  FORMULAS  FOR  FLEXURAL  MODULUS  OF  A SOLID 
SYMMETRIC  [0/90]  CROSS-PLY  LAMINATE 


1 


/ 

02 

03 

D„* 

u, 

3/m 

/ 

u, 

-3/m 

/ 

0,2 * 

04 

-/ 

Dee* 

05 

-/ 

D,e=  D*e=  O 

Djf  = 12 Djj/h3 , z*  =0 
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Note  that  only  V is  affected  by  the  stacking  sequence.  We  only  showed  the 
case  of  changing  the  number  of  ply  assemblies.  Other  stacking  sequences 
than  those  shown  in  Fig.  68  are,  of  course,  possible;  an  example  of  which 
maybe  [°2^9°4^02]s*  The  value  for  vt  wiU  be  different  from  that  shown  in 

Eq.  258  and  Table  48.  The  effect  of  V on  the  flexural  modulus  is  the  degree 
of  anisotropy,  or  the  difference  between  and  D In  the  limit  when 
we  have  infinite  numbers  of  alternating  plies,  our  laminate  will  become 
qua  si -homogeneous.  The  property  of  the  laminate  will  be  square  symmetric, 
not  isotropic.  This  difference  between  square  symmetric  and  isotropy  was 
illustrated  in  Eq.  22  and  23. 

Let  us  calculate  the  flexural  modulus  of  cross-ply  laminates  as 

shown  in  Fig.  68.  The  formulas  for  this  class  of  special  laminates 

are  listed  in  Table  48.  Using  the  data  for  T300/5208,  we  have  for  16 -ply 

laminates  with  no  core,  or  z =0. 

c 

h = 16h  = 16  x 125  X lo"6  = 2xlo'3m 


+ 


3U2/m  + U 


>] 


(262) 


or 


or 


= 666  x 10 ~12J^76.  37  + 3 x 85.  73/m  + I9.7lj 
= 666  x 10"12  x 160.  3 x 109 


= 106.  8 Nm  for 
= 85.  4 Nm  " 

= 74.  7 Nm  " 

= 64.  0 Nm  for 


4-ply  assembly  laminate  or  m=4 

8 " " " m=8 

16  " " " m=16 


quasi-homogeneous  laminate  or  m=00 


(263) 


Similarly,  3 

D22  * IT  [U1  ' 3V“  + U3]  <264> 

= 21.16  Nm  for  4-ply  assembly  laminate  or  m=4 
or  = 42.  58  Nm"  8 " " " m=8 

or  = 53.28  Nm"  16  " " » m=16 

or  = 63.99  Nm  " quasi -homogeneous  laminate 

Components  of  an<i  are  not  dependent  on  or  m,  the  number 
of  ply  assemblies.  So  their  constant  values  are  1.93  and  4.77  Nm, 
respectively. 


Du,  Dm>  in  Nm 


Fig.  69.  Flexural  modulus  components  as  functions  of  ply  assemblies 
for  a T300/5208  laminate.  Note  that  as  ply  assemblies  m 
increases,  the  modulus  components  approach  the  modulus  of 
the  quasi -homogeneous  laminate,  although  many  assemblies 
are  needed  for  good  convergence. 


Fig.  69  shows  the  sensitivity  of  the  two  normal  components  of  the 
flexural  modulus  of  T300/5208  as  functions  of  the  number  of  ply  assemblies. 
The  convergence  for  this  case  does  not  appear  to  be  very  rapid. 


If  we  introduce  a honeycomb  core  into  our  cross -ply  laminate,  we 
want  to  show  how  the  flexural  modulus  can  be  calculated.  Let  us  examine 
three  cross -ply  laminates  in  Fig.  70.  These  laminates  are  sandwich  con- 
structions with  facing  materials  identical  to  those  solid  laminates  shown 
in  Fig.  68.  The  number  of  ply  assemblies  are  different  among  these  lami- 
nates. The  core- half-depth  is  equal  to  four  plies. 


(a)  (b)  (c) 

Fig.  70.  Cross-ply  sandwich  laminates.  This  symmetric  laminate  has 
2-8  ply  facings  and  4-ply  thick  half- depth  of  core.  Total 
thickness  of  laminate  is  24  equivalent  plies.  Three  different 
numbers  of  ply  assemblies  are  shown;  m=4,  8 and  16.  This 
figure  shows  the  same  facing  laminates  as  those  solid  laminates 
in  Fig.  68. 
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The  flexural  modulus  of  these  sandwich  laminates  can  be  readily  cal 
culated  by  substituting  the  nonzero  trignometric  functions  into  Eq.  251. 
For  the  case  of  2 -ply  assembly  laminate  in  Fig.  70(a),  or  m = 4. 


We  notice  the  trend  again  as  the  number  of  ply  assemblies  m increase 
* 

the  value  of  V decreases  by  the  following  relation: 


The  sandwich  laminates  approaches  square -symmetric  as  m increases 

We  need  the  following  values  before  we  can  use  the  formulas  for 
the  flexural  modulus  in  Table  44. 


Note  this  is  the  same  as  that  in  Eq.  270.  This  is  expected  for  homogeneous 
material  shown  in  Table  46.  The  same  correction  factor  for  the  sandwich 
core  is  applied  to  the  first  column  of  Table  44.  Summarizing  the  results  in 
Eq.  269  to  272,  we  can  show  the  formulas  for  the  sandwich  laminates  in 
Fig.  70  in  a matrix  multiplication  table  as  follows: 


Table  49.  FORMULAS  FOR  FLEXURAL  MODULUS  OF 
A SYMMETRIC  SANDWICH  LAMINATE 
WITH  [0/90]  FACINGS 


/6/9m 


!6/9m 


We  are  now  ready  to  calculate  the  flexural  modulus  of  our  sandwich 
laminates  assuming  the  facing  material  is  T300/5208. 


From  Table  43  and  Table  49  for  m*4  in  Fig.  67(a); 


26  x 76.  37/27  + 16  x 85.  73/9  x 4 + 26  x 19.  71/27 


Similarly,  from  Table  49 


This  laat  value  is  the  modulus  for  a quasi -homogeneous  laminate.  This  is 
the  same  value  to  which  will  also  converge.  Thus,  we  will  have  a 
square -symmetric  laminate  in  flexure.  The  results  of  these  components  of 
the  flexural  modulus  will  be  plotted  in  Fig.  71. 


Note  the  substantial  increase  in  the  modulus  components  of  the  sand 
wich  construction  here  over  the  solid  laminates  shown  in  Fig.  69.  First 
of  all,  there  is  a thickness  increase  from  16  to  24  plies.  If  our  laminate 
were  homogeneous  and  solid  or  without  a core,  the  increase  in  the 
flexural  modulus  components  will  be  the  cube  of  the  thickness  ratio.  In 
our  particular  case,  it  will  be: 


(24/16)  = 3.375 


208(fb  m=ao 


Flexural  modulus  for  a sandwich  laminate  of  T300/5208  as 
functions  of  the  number  of  ply  assemblies.  The  trend  is  very 
similar  to  that  shown  in  Fig.  69  for  solid  laminates.  The 
dashed  lines  are  bound  on  the  modulus  in  Fig.  69  corrected 
for  thickness  and  core  depth.  Note  the  approach  by  direct 
correction,  although  valid  for  homogeneous  facing  materials, 
is  not  so  for  discrete  laminates.  Calculation  for  specific 
stacking  sequence  of  the  facing  is  necessary. 


• '■>  I 


On  the  other  hand,  if  we  have  a sandwich  construction,  there  should  be  a 
reduction  proportional  to 


which  represents  the  effect  of  core  if  the  facings  were  homogeneous. 
Assuming  the  core  thickness  for  our  laminate  is  the  same  as  those  in 
Fig.  70,  the  value  of  the  core  correction  in  Eq.  281  will  be  26/27  as  in 
Eq.  272.  Thus,  the  net  effect  of  thickness  increases  and  the  presence  of 
core  is  simply  the  product  of  Eq.  280  and  281: 


3.  375  x 26/27  = 3.  25 


The  results  are  shown  as  dashed  lines  in  Fig.  71.  Note  the  difference 
between  the  dashed  and  solid  lines.  The  point  is  that  the  flexural  modulus 
of  multidirectional  laminates  must  be  calculated  for  each  stacking  sequence 
and  core  depth.  Simple  ratios  such  as  that  in  Eq.  282  are  applicable  only 
to  sandwich  construction  with  homogeneous  facing  material.  Eq.  252 
for  example,  provide  the  ratios  for  thickness  and  core  depth  corrections. 
But  this  is  limited  to  unidirectional  material  only.  Laminated  composites 
must  be  assessed  on  individual  basis.  Attempts  to  simplify  the  calculation 
°f  flexural  modulus  by  smearing  the  laminated  facing  can  lead  to  erroneous 
results.  Since  precise  calculation  is  straight  forward,  such  shortcuts  are 
not  warranted. 


6.  FLEXURAL  MODULUS  OF  SYMMETRIC  LAMINATES 

In  the  last  section  we  saw  that  cross-ply  laminates  are  orthotropic, 
or  square  symmetric  when  the  number  of  ply  assemblies  approach  infinity. 

We  will  see  in  this  section  that  symmetric  laminates  with  or  without  core 

are  generally  anisotropic.  Balanced  laminates  such  as  an  angle-ply  laminate  are 

orthotropic  in  their  in-plane  modulus  but  are  generally  anisotropic  in  their 


flexural  modulus.  The  reason  for  this  is  the  fact  that  the  position  of  each 
ply  is  unique.  Each  ply  occupies  unique  values  along  the  z-axis.  The 
shear  coupling  terms  and  of  a +0  ply  cannot  be  cancelled  by  those 
of  a -0  ply  unless  the  positions  of  these  plies  are  judiciously  selected. 

We  will  show  later  that  the  shear  coupling  terms  can  be  cancelled  if 
antisymmetric  laminates  are  used.  So  there  are  two  methods  of 
obtaining  orthotropic  flexural  modulus: 

• Use  of  on-axis  plies  only.  This  is  the  case  of  cross-ply 
laminates. 

• Use  of  antisymmetric  laminates.  This  will  be  discussed 
later. 

Th^notivation  to  make  laminates  orthotropic  (and  symmetric)  is  often 
driven  by  the  availability  of  the  in-plane  stress  analysis  and  flexural  analysis 
of  laminates.  Most  current  analytical  tools  are  limited  to  orthotropic  and 
homogeneous  plates.  It  seems  ridiculous  to  limit  the  use  of  composite 
materials  by  the  availability  of  analytical  tools.  Although  it  is  not  within 
the  scope  of  the  introductory  text  to  cover  the  stress  analysis  of  anisotropic 
bodies,  it  is  important  to  understand  how  anisotropy,  and  nonhomogeneity 
emerge  in  composite  laminates  and  to  what  degree  they  can  be  manipulated. 

To  rule  out  anisotropy,  for  example,  forces  us  to  impose  arbitrary 
limits  on  the  use  of  composite  materials.  In  this  case,  composite  materials 
demand  more  advanced  analytical  tools  than  ordinary  materials.  In  many 
finite-element  techniques,  these  features  are,  or  can  be,  incorporated.  But 
much  more  progress  is  needed  before  we  can  feel  comfortable  with  anisotropy. 


7.  PLY  STRESS  AND  PLY  STRAIN  ANALYSIS 

The  ply  stress  and  ply  strain  in  a symmetric  laminate  with  or  without 
core  due  to  flexure  can  be  determined  following  the  procedure  in  Section  IV 
for  the  in-plane  stretching.  This  is  shown  in  Fig.  72  which  is  analogous  to 
Fig.  50  for  the  in-plane  behavior. 


Moment  Curvature  In  ~ Plane  On -Axis  On-Axii 

Strain  pty  piy 

Strain  Stress 

Fig.  72.  Ply  stress  and  strain  in  a symmetric  laminate  under  flexure: 

From  (a)  to  ( b):  Use  moment- curvature  relations  in  Table  42. 

From  (b)  to  (c):  Use  curvature- strain  equation  in  Eq.  216. 

Use  top  surface  of  each  ply  assembly  for  the 
z- value;  i.  e. , z^  for  the  i-th  ply  assembly 
shown  in  Fig.  64. 

From  (c)  to  (d):  Use  strain  transformation  to  transform  the 
laminate  strain  to  the  on-axis  strain  on  the 
top  surface  of  each  ply  assembly. 

From  (d)  to  (e):  Use  the  on-axis  stress- strain  relation  to 
determine  the  corresponding  on-axis  ply 
stress. 


The  process  of  determining  the  ply  stress  and  ply  strain  is  straight 
forward.  The  motivation  is  to  assess  the  strength  of  each  ply  within  the 
laminate.  The  strength  calculation  in  terms  of  strength  ratios  will  be  covered 
in  a later  section.  We  have  shown  in  Fig.  72  that  the  highest  value  of  z 
be  used  for  the  ply  strain  determination  as  we  go  from  (b)  to  (c).  The  highest 
z,  or  z.  for  the  i-th  ply  assembly  in  Fig.  64,  will  have  the  highest  strain 
and  stress  in  each  ply  assembly  by  virtue  of  Eq.  216.  It  is  the  highest  strain 
components  that  will  govern  the  strength  of  that  ply  assembly.  This  is  obvious 
if  all  components  of  the  curvature  have  the  same  sign,  positive  or  negative. 

But  is  it  not  so  obvious  if  the  signs  are  mixed.  We  can  show  that  regardless 
of  the  signs  of  the  curvature  components,  the  highest  z will  govern  their 
strength.  This  will  be  done  in  the  section  on  strength,  to  be  presented  later. 

Future  volumes  shall  include  the  following  topics: 

General  Laminates 
Thermal  Stresses 
Strength  Analysis 
Mic  r omec  hanic  s 
Fracture 
Fatigue 
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